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Abstract 

In this paper, we first introduce the concept of Laurent differentiaUy essential systems 
and give a criterion for Laurent differentiahy essential systems in terms of their supports. 
Then the sparse differential resultant for a Laurent differentially essential system is de- 
fined and its basic properties are proved. In particular, order and degree bounds for the 
sparse differential resultant are given. Based on these bounds, an algorithm to compute 
the sparse differential resultant is proposed, which is single exponential in terms of the 
number of indeterminates, the Jacobi number of the system, and the size of the system. 
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1 Introduction 

The multivariate resultant, which gives conditions for an over-determined system of polyno- 
mial equations to have common solutions, is a basic concept in algebraic geometry [12\ [T9| 
[23| [26} [27l HU [49]. In recent years, the multivariate resultant is emerged as one of the most 
powerful computational tools in elimination theory due to its ability to eliminate several 
variables simultaneously without introducing much extraneous solutions. Many algorithms 
with best complexity bounds for problems such as polynomial equation solving and first 
order quantifier elimination, are based on the multivariate resultant [U O [HI [151 US]- 

In the theory of multivariate resultants, polynomials are assumed to contain all the 
monomials with degrees up to a given bound. In practical problems, most polynomials are 
sparse in that they only contain certain fixed monomials. For such sparse polynomials, 
the multivariate resultant often becomes identically zero and cannot provide any useful 
information. 

As a major advance in algebraic geometry and elimination theory, the concept of sparse 
resultant was introduced by Gelfand, Kapranov, Sturmfels, and Zelevinsky jl9|, I49j . The 
degree of the sparse resultant is the Bernstein-Kushnirenko-Khovanskii (BKK) bound [2] 
instead of the Bezout bound [El [40l [50], which makes the computation of the sparse resul- 
tant more efficient. The concept of sparse resultant is originated from the work of Gelfand, 
Kapranov, and Zelevinsky on generalized hyper geometric functions, where the central con- 
cept of ^-discriminant is studied |17j . Kapranov, Sturmfels, and Zelevinsky introduced the 
concept of ^-resultant |28j . Sturmfels further introduced the general mixed sparse resultant 
and gave a single exponential algorithm to compute the sparse resultant [l9l[50]. Canny and 
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Emiris showed that the sparse resultant is a factor of the determinant of a Macaulay style 
matrix and gave an efficient algorithm to compute the sparse resultant based on this matrix 
representation [13^ I14j. D'Andrea further proved that the sparse resultant is the quotient of 
two Macaulay style determinants similar to the multivariate resultant 

Using the analogue between ordinary differential operators and univariate polynomials, 
the differential resultant for two linear ordinary differential operators was implicitly given 
by Ore [39] and then studied by Berkovich and Tsirulik [l] using Sylvester style matrices. 
The subresultant theory was first studied by Chardin [7] for two differential operators and 
then by Li [38] and Hong [24] for the more general Ore polynomials. 

For nonlinear differential polynomials, the differential resultant is more difficult to define 
and study. The differential resultant for two nonlinear differential polynomials in one variable 
was defined by Ritt in [441 p. 47]. In [55^ p. 46], Zwillinger proposed to define the differential 
resultant of two differential polynomials as the determinant of a matrix following the idea 
of algebraic multivariate resultants, but did not give details. General differential resultants 
were defined by Carra-Ferro using Macaulay's definition of algebraic resultants [6]. But, 
the treatment in [6] is not complete. For instance, the differential resultant for two generic 
differential polynomials with positive orders and degrees greater than one is always identically 
zero if using the definition in [6]. In [54) . Yang, Zeng, and Zhang used the idea of algebraic 
Dixon resultant to compute the differential resultant. Although efficient, this approach is 
not complete, because it is not proved that the differential resultant can always be computed 
in this way. Differential resultants for linear ordinary differential polynomials were studied 
by Rueda-Sendra [Ml [47] . In [16] , a rigorous definition for the differential resultant of n + 1 
differential polynomials in n variables was first presented and its properties were proved. A 
generic differential polynomial with order o and degree d contains an exponential number 
of differential monomials in terms of o and d. Thus it is meaningful to study the sparse 
differential resultant which is the main focus of this paper. 

Our first observation is that the sparse differential resultant is related with the non- 
polynomial solutions of algebraic differential equations, that is, solutions with non-vanishing 
derivatives to any order. As a consequence, the sparse differential resultant should be more 
naturally defined for Laurent differential polynomials. This is similar to the algebraic sparse 
resultant [191 150j . where non-zero solutions of Laurent polynomials are considered. 

Consider n + 1 Laurent differential polynomials in n differential indeterminates Y = 
{yi, ■ ■ • ,yn}: 

k 

r, = ^UikM,k{i = 0,...,n), (1) 

fc=0 

where Uik S £ are differentially independent over Q and Afj^ are Laurent differential mono- 
mials in Y. As explained later in this paper, we can assume that Mik are monomials with non- 
negative exponent vectors Uik- Let Sj = ord(Pi, Y) and denote Mjfc/Mjo = Y[]=i YliLoivf^y'''^' 
A (Y[si] -^a,i:-aio ^ where yj''' is the l-th. derivative of yj and yI"*'] is the set {yj''' : 1 < j < 
n, < / < Si}. Let Uj = (ujo, un, . . . , uuj = 0, . . . , n) be the coefficient vector of Pj. 

The concept of Laurent differentially essential system is introduced, which is a necessary 
and sufficient condition for the existence of sparse differential resultant. Pq, . . . , Pn are called 
Laurent differentially essential if [Pq, . . . , Pn] n Q{uo . . . , u„} is a prime differential ideal of 
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codimension one, where [Pq, . . . , P™] is a differential ideal generated in Q{Y, Y^^; uq, . . . , u„}. 
This concept is similar to (but weaker than) the concept of essential supports introduced by 
Sturmfels in [50], but its properties are more complicated. Precisely, we have 

h Ij 

Theorem 1.1 For Pj given in (OP, let qj = max"^Qord(Pi, y^) and dij = '^ik Yl tikji^j 

k=0 1=0 

(i = 0, . . . ,n; j = 1, . . . ,n) where xj are algebraic indeterminates. Denote 

( C^Ol dQ2 • • • d^n \ 
d\\ d\2 • • • d\n 

\ dnl dn2 ■ ■ ■ dnn / 

to be the symbolic support matrix of Then the following assertions hold. 

1) The differential transcendence degree of Q{uq . . . , Un){-^^, ■ ■ ■ , -^f^) over Q(uo . . . , u„) 
equals rk(Mp). 

2) [Pq, . . . , Pn] nQ{uo, . . . , u„} is a prime differential ideal of codimension n + 1 — rk(Mp). 
So Pq, • • ■ ,Pn form a Laurent differentially essential system if and only i/rk(Mp) = n. 

3) Pq, . . . , Pn form a Laurent differentially essential system if and only if there exist ki (1 < 
ki < li) such that rk(Mfc(j^,,,^fc^) = n where is the symbolic support matrix for 
the Laurent differential monomials Mq/^q/Mqq, . . . , M^kn/Mno- 

With the above theorem, computing the differential transcendence degree of certain differ- 
ential polynomials is reduced to computing the rank of certain symbolic matrix. Similar to 
the case of linear equations, this result provides a useful tool to study generic differential 
polynomials. As an application of the above result, the differential dimension conjecture 
[451 p. 178] for a class of generic differential polynomials is proved. For the n + 1 Laurent 
differential monomials Mo^q/Moo, • • • ,M„fc^/M„o (1 < < U) mentioned in 3) of Theorem 
ll.H a more efficient algorithm to compute their differential transcendence degree over Q is 
given by reducing their symbolic support matrix to a standard form called T-shape. 

Before introducing the properties of the sparse differential resultant, the concept of 
Jacobi number is given below. Let G = {51,..., (7^} be n differential polynomials in 
Y = {2/1, . . . Let Sij = ord{gi,yj) be the order of gi in yj if yj occurs effectively in 

fi and Sij = —00 otherwise. Then the Jacobi bound, or the Jacobi number, of G, denoted as 
Jac(G), is the maximum number of the summations of all the diagonals of S. Or equivalently, 

n 

Jac(G) = max^Si(^(j), 

i=l 

where a is a permutation of {1, . . . ,n}. Jacobi's Problem conjectures that the order of the 
zero dimensional component of G is bounded by the Jacobi number of G |43j . 

The properties of the sparse differential resultant are summarized in the following theo- 
rem. 
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Theorem 1.2 The sparse differential resultant R(uo, . . . , u„) G Q{uo, . . . , u„} o/Pq, . . . , Pn 
has the following properties. 

1) R(uo, . . . , u„) is differentially homogenous in each Uj (i = 0, . . . , n). 

2) h, = ord(R, u,) < J, = Jac(P.) where P. = {P^, . . . ,P^}\{Pf }. 

3) Let ^o(IPo7 • • • ) IPn) ihe set of all specializations of the coefficients Uik ofV-i under 
which Pj = (i = 0, . . . , n) have a common non-polynomial solution and 2o(Po, • • • , Pn) 
the Kolchin differential closure of Zq{^q, • • • , Pn). Then Zo(Po, • • • , P„) = V(sat(R)) . 

4) Assume that Pj (z = 0, . . . , n) have the same set A of monomials. The differential toric 
variety Xj[ associated with A is defined and is shown to be an irreducible projective 
differential variety of dimension n. Furthermore, the differential Chow form US\. \37\l 
of is R. 

5) (Poison Type Product Formula) Let uq appear in R and to = deg(R, Uqq"^). Then 
there exist ^^-fc in certain differential field Jv (t = ^, ■ ■ ■ ^to) such that 

R = ^n(noo + 5^uofce.fc)(^°\ 

T = l k = l 

where A is a polynomial in Q(ui, . . . , u„)[uq''''\'Uqq"^]. Furthermore, if 1) any n of the 
Pj (i = 0, . . . , n) form a differentially independent set over Q(uo, . . . , u„) and 2) for 
each j = 1, . . . ,n, Bj £ Span^jajfe — Ojo : k = 1, . . . ,li;i = 0, . . . , n}, then there exist 
Tjrk £ -Tv {t = 1, . . . ,tQ; k = \, . . . ,n) such that 

to 

R=^n 

T = l 

where rjr = {r]ri, ■ ■ ■ ,r]rn) and e-i is the exponent vector of yi. Moreover, r]r {t = 
1, . . . , to) (^fc generic points of the prime differential ideal [Pi, ... , P„] : m C -7^(uo, . . . , 
u„){Y}, where m is the set of all differential monomials in Y. 

6) deg(R) < nr=o("^« + 1)^'^^ < (m + l)ELo(-^^+i) < (m + where rm = 
deg(Pi, Y), m = maxj{mj}, and J = Ya=q Ji- 

7) Let Si = ord(Pi, Y). Then R has a representation 

n n hi 

H M<J^+i)'^^^(^) . R = ^ £ G., (P.) 

i=0 i=0 j=0 

where dj € Q[uf . . . , uL^"^, YI'*]] with h = max{hi + Si} such that deg(Gij(Pi)(^')) < 
+ l + TJl=o{h, + l)deg(Mio)]deg(R). 



Po(??r) 

Moo(??r) 
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Although similar to the properties of algebraic sparse resultants, each property given 
above is an essential extension of its algebraic counterpart. For instance, it needs lots of 
efforts to obtain the Poison type product formula. Property 2) is unique for the differential 
case and reflects the sparseness of the system in certain sense. 

More properties for the sparse differential resultant are proved in this paper. For instance, 
the explicit condition for the equation system ([1]) to have a unique solution for Y is given. The 
sparse resultant for differential polynomials with non-vanishing degree terms are also defined, 
which gives conditions for the existence of solutions instead of non-polynomial solutions. 

Let Pj (i = 0, . . . , n) in ([T]) be generic differential polynomials containing all monomials 
with order < Si and degree < rrii and R(uo, . . . ,u„) the differential resultant of Pq, . . . ,Pn- 
Then a BKK style degree bound is given: 

Theorem 1.3 For each i € {0, 1, . . . , n}, 

deg(R, Ui) < M{{Qji)j^ifl<i<s-sp Qio, Qi,fc-i, Qi,k+i, Qi,s-s,) 

k=0 

where Qji is the Newton polytope of (Pj)^') as a polynomial in y^^\ . . . ,y^n o-nd A4{S) is the 
mixed volume for the polytopes in S. 

In principle, the sparse differential resultant can be computed with the characteristic 
set method for differential polynomials via symbolic computation [l5l EJ EJ SSI \52\- But 
in general, differential elimination procedures based on characteristic sets do not have an 
elementary complexity bound [20]. 

Based on the order and degree bounds given in 2) and 6) of Theorem 11.21 & single 
exponential algorithm to compute the sparse differential resultant R is proposed. The idea 
of the algorithm is to compute R with its order and degree increasing incrementally and to 
use linear algebra to find the coefficients of R with the given order and degree. The order 
and degree bounds serve as the termination condition. Precisely, we have 

Theorem 1.4 With notations introduced in Theorem \1.SX the sparse differential resultant 

o/Pq, . . . ,Pn can be computed with at most 0[- ) Q-arithmetic 

operations, where I = Yl^=oi^i + 1), m = max"^Qmj, and J = Y17=o 

From Theorem 11.41 the complexity of this algorithm is single exponential in terms of I, and 
J. The sparseness is reflected in the quantity / which is called the size of the system and the 
Jacobi number J. Note that even the complexity of computing the algebraic sparse resultant 
is single exponential |491 114j. The algorithm seems to be the first one to eliminate several 
variables from nonlinear differential polynomials with a single exponential complexity. 

The rest of the paper is organized as follows. In Section 2, preliminary results are 
introduced. In Section 3, the sparse differential resultant for Laurent differentially essential 
systems is defined. In Section 4, Theorem ll.il is proved. In Section 5, properties 1) - 5) of 
Theorem 11.21 are proved. In Section 6, properties 6) and 7) of Theorem 11.21 and Theorems 
11.31 and 11.41 are proved. In Section 7, the paper is concluded and several unsolved problems 
for differential sparse resultant are proposed. 
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2 Preliminaries 



In this section, some basic notations and preliminary results in differential algebra will be 
given. For more details about differential algebra, please refer to \A5 \ 1291 ^ I48 t [IE]. 

2.1 Differential polynomial algebra and Kolchin topology 

Let be a fixed ordinary differential field of characteristic zero, with a derivation operator 
6. An element c £ J- such that (5c = is called a constant of J-. In this paper, unless 
otherwise indicated, 6 is kept fixed during any discussion and we use primes and exponents 
(i) to indicate derivatives under 6. Let denote the free commutative semigroup with unit 
(written multiplicatively) generated by d. 

A typical example of differential field is which is the field of rational functions in a 
variable x with 5 = 

Let 5 be a subset of a differential field Q which contains F. We will denote respectively 
by F[S]^ ^(S), J^{S}, and J^{S) the smallest subring, the smallest subfield, the smallest 
differential subring, and the smallest differential subfield of Q containing T and S. If we 
denote Q{S) to be the smallest subset of Q containing S and stable under 5, we have T{S} = 
T[Q{S)] and T{S) = T{Q{S)). A differential extension field ^ of is said to be finitely 
generated if Q has a finite subset S such that Q = J^{S). 

A subset E of a differential extension field G of J- is said to be differentially dependent 
over if the set (6'a)e(=e,oes is algebraically dependent over and is said to be differentially 
independent over J", or to be a family of differential indeterminates over in the contrary 
case. In the case S consists of one element a, we say that a is differentially algebraic or 
differentially transcendental over respectively. The maximal subset of Q which are 
differentially independent over T is said to be a differential transcendence basis of G over 
T. We use d.tr.degt?/J^ (see [291 P- 105-109]) to denote the differential transcendence degree 
of Q over J^, which is the cardinal number of Vt. Considering J-" and Q as ordinary algebraic 
fields, we denote the algebraic transcendence degree of Q over F by iiAegQ/F. 

A homomorphism (/? from a differential ring {TZ,5) to a differential ring {S,5i) is a dif- 
ferential homomorphism if ip o 5 = 5i o ip. If TZq is a common differential subring of TZ and 
S and the homomorphism if leaves every element of TZq invariant, it is said to be over TZq. 
If, in addition TZ is an integral domain and 5 is a differential field, (p is called a differential 
specialization of TZ into S over Rq. The following property about differential specialization 
will be needed in this paper, which can be proved similarly to Theorem 2.16 in [16]. 

Lemma 2.1 Let Pi(U,Y) G 7'(Y){U} {i = l,...,m) where U = (ni,...,n^) and Y = 

{yi, . . . , yn) are sets of differential indeterminates. If the set (Pi(U, Y))^'^'^^ (i = 1, . . . , m; j = 
1, . . . , Hi) are algebraically dependent over F{V), then for any differential specialization U to 
U'' C overF, (Pj(U'^, Y))('^»j)(i = l,...,m;j = l,...,ni) are algebraically dependent over 
F. In particular, i/Pj(U, Y) (z = 1, . . . ,m) are differentially dependent over F{I]), then for 
any differential specialization V to V C F over F, Pi(V,Y) (i = 1, . . . ,m) are differentially 
dependent over F . 

A differential extension field £" of is called a universal differential extension fields if 
for any finitely generated differential extension field F\ of in £" and any finitely generated 
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differential extension field T2 of T\ not necessarily in f , can be embedded in E over J^i, 
i.e. there exists a differential extension field J-3 in E that is differentially isomorphic to T2 
over T\. Such a differential universal extension field of T always exists ([Ml Theorem 2, p. 
134]). By definition, any finitely generated differential extension field of J- can be embedded 
over T into and £^ is a universal differential extension field of every finitely generated 
differential extension field of J- . In particular, for any natural number n, we can find in E 
a subset of cardinality n whose elements are differentially independent over T . Throughout 
the present paper, E stands for a fixed universal differential extension field of T . 

Now suppose Y = {yi,y2, ■ ■ ■ , Un} is a set of differential indeterminates over E. For any 
y G Y, denote by y*^^-*. The elements of -7^{Y} = F[y^l^^ : j = 1, . . . , n; /c G N] are called 
differential polynomials over J-" in Y, and J-"{Y} itself is called the differential polynomial ring 
over T in Y. A differential polynomial ideal X in J-'{Y} is an ordinary algebraic ideal which 
is closed under derivation, i.e. 5{X) C X. And a prime (resp. radical) differential ideal is 
a differential ideal which is prime (resp. radical) as an ordinary algebraic polynomial ideal. 
For convenience, a prime differential ideal is assumed not to be the unit ideal in this paper. 

By a differential affine space we mean any one of the sets E^ (n G N). An element 
r/ = (?7i, . . . ,r/„,) of will be called a point. Let S be a subset of differential polynomials 
in J^{Y}. A point r/ = (r/i, . . . , rjn) G is called a differential zero of S if f{r]) = for any 
/ G S. The set of differential zeros of S is denoted by V($]), which is called a differential 
variety defined over J^. The differential varieties in E^ (resp. the differential varieties in E^ 
that are defined over T) are the closed sets in a topology called the Kolchin topology (resp. 
the Kolchin J^-topology). 

For a differential variety V which is defined over J^, we denote ^.{V) to be the set of all 
differential polynomials in J-'{Y} that vanish at every point of V . Clearly, I(V') is a radical 
differential ideal in J^{Y}. And there exists a bijective correspondence between Kolchin J-"- 
closed sets and radical differential ideals in J^{Y}. That is, for any differential variety V 
defined over J-", V(I(y)) = V and for any radical differential ideal X in J^{Y}, I(V(I)) = X. 

Similarly as in algebraic geometry, an J^-irreducible differential variety can be defined. 
And there is a bijective correspondence between J^-irreducible differential varieties and prime 
differential ideals in J-'{Y}. A point G V(X) is called a generic point of a prime ideal 
X C J-'{Y}, or of the irreducible variety V(X), if for any polynomial P G J-'{Y} we have 
P{rj) = <^=> P G X. It is well known that p. 27] a non-unit differential ideal is prime if 
and only if it has a generic point. 

Let Z be a prime differential ideal in J^{Y} and ^ = (Ci,---,^™) a generic point of 
X |291 p. 19]. The dimension of X or of V(X) is defined to be the differential transcen- 
dence degree of the differential extension field ^(Cii • • • 1 ■^n) over F, that is, dim(X) = 
d.tr.degJ'(6,---,^n)/-7^- 

We will conclude this section by introducing some basic concepts in projective differential 
algebraic geometry which will be used in Section [5^31 For more details, please refer to [33|.i36j. 
And unless otherwise stated, in the whole paper, we only consider the affine differential case. 

For each / G N, consider a projective space P(/) over E. By a differential projective 
space we mean any one of the sets P(Z) (/ G N). Denote zq, zi, . . . , to be the homogenous 
coordinates. Let X be a differential ideal of J^{z} where z = {zq, zi, . . . , zi}. Denote X : z = 
{/ G Flzjl Zjf G Xfor each j = 0, . . . , /}. 
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Definition 2.2 LetX he a differential ideal ofJ-{z}. I is called a differentially homogenous 
differential ideal of J^{z} if X : z, = X and for every P G X and a differential indeterminate 
A over -P(Az) G F{X\X in the differential ring J^{A, z}. 

Consider a differential polynomial P G £{z} and a point a G P(Z)- Say that P vanishes 
at a, and that a is a zero of P, if P vanishes at \a for every A in £. For a subset ^ of P(Z)) 
let I(^) denote the set of differential polynomials in J-{x} that vanishes on ^ . Let N{S) 
denote the set of points of P(i) that are zeros of the subset S of £^{z}. And a subset V of 
P(Z) is called a projective differential T -variety if there exists S C .F{z} such that V = Y{S). 
There exists a one-to-one correspondence between projective differential varieties and perfect 
differentially homogenous differential ideals. And a projective differential J^-variety V is ^- 
irreducible if and only if I{V) is prime. 

Let X be a prime differentially homogenous ideal and ^ = (^o, ^i, • • • , 6) be a generic point 
of X with ^0 7^ 0- Then the differential dimension of V(X) is defined to be the differential 
transcendence degree of J-'{{^Q^^k)i<k<l) over T. 

2.2 Characteristic sets of a differential polynomial system 

Let / be a differential polynomial in J-'{Y}. We define the order of / w.r.t. yi to be the 
greatest number k such that appears effectively in /, which is denoted by ord(/, j/j). 
And if yi does not appear in /, then we set ord(/, yi) = — oo. The order of / is defined to 
be maxj ord(/, yj), that is, ord(/) = maxj ord(/, yj). 

A ranking ^ is a total order over 0(Y), which is compatible with the derivations over 
the alphabet: 

1) 69yj > 9yj for all derivatives 9yj G 6(Y). 

2) em > e2yj =^ sem > se^y, for f?iy„92%- e e(Y). 

By convention, 1 < 9yj for all 9yj G 0(Y). 

Two important kinds of rankings are the following: 

1) Elimination ranking: yi > yj =^ S'^yi > 5^yj for any k,l > 0. 

2) Orderly ranking: k > I =^ S^yi > ^'y^, for any i,j G {1, 2, . . . . r/}. 

Let be a differential polynomial in J-'{Y} and ^ a ranking endowed on it. The greatest 
derivative w.r.t. ^ which appears effectively in p is called the leader of p, which will be 
denoted by Up or ld(p). The two conditions mentioned above imply that the leader of 9p is 
9up for G 0. Let the degree of p in Up be d. As a univariate polynomial in Up, p can be 
rewritten as 

p = hu'l + h-iu'^-^ + • • • + lo. 

Id is called the initial of p and is denoted by Ip. The partial derivative of p w.r.t. Up is 
called the separant of p, which will be denoted by Sp. Clearly, Sp is the initial of any proper 
derivative of p. The rank of p is Up, and is denoted by rk(p). 

Let p and q be two differential polynomials and Up the rank of p. q is said to be partially 
reduced w.r.t. p if no proper derivatives of Up appear in q. q is said to be reduced w.r.t. p if q 
is partially reduced w.r.t. p and deg{q, Up) < d. Let .4, be a set of differential polynomials. A 
is said to be an auto-reduced set if each polynomial of A is reduced w.r.t. any other element 
of A. Every auto-reduced set is finite. 
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Let A = Ai, A2, ■ ■ ■ , At be an auto-reduced set with Si and Ij as the separant and initial 
of Ai, and / be any differential polynomial. Then there exists an algorithm, called Ritt's 
algorithm of reduction, which reduces / w.r.t. ^ to a polynomial r that is reduced w.r.t. A, 
satisfying the relation 

t 

nsfl^^-/ = r,mod[^], 

i=l 

where di, Cj (i = 1, 2, . . . , i) are nonnegative integers. The differential polynomial r is called 
the differential remainder of / w.r.t. A. 

Let A be an auto-reduced set. Denote H_4 to be the set of all the initials and separants 
of A and to be the minimal multiplicative set containing H_4. The saturation ideal of A 
is defined to be 

sat(^) = [A] ■.H^ = {p:3h£ H^, s.t. hp G [A]}. 

An auto-reduced set C contained in a differential polynomial set S is said to be a charac- 
teristic set of 5, if 5 does not contain any nonzero element reduced w.r.t. C. A characteristic 
set C of an ideal J' reduces to zero all elements of J'. If the ideal is prime, C reduces to zero 
only the elements of and = sat(C) ( \29\ Lemma 2, p. 167]) is valid. 

In terms of the characteristic set, the cardinal number of the characteristic set of I is 
equal to the codimension of I, that is n — dim(Z). When I is of codimension one, it has the 
following property. 

Lemma 2.3 14^' Let X he a prime differential ideal of codimension one in J-'{Y}. 

Then there exists an irreducible differential polynomial A such that X = sat(^) and {A} is 
the characteristic set ofX w.r.t. any ranking. 

3 Sparse differential resultant for Laurent differential poly- 
nomials 

In this section, the concepts of Laurent differential polynomials and Laurent differentially 
essential systems are first introduced, and then the sparse differential resultant for Laurent 
differentially essential systems is defined. 

3.1 Laurent differential polynomial 

Let T be an ordinary differential field with a derivation operator 5 and J-'{Y} the ring 
of differential polynomials in the differential indeterminates Y = {yi, . . . , Let £ be 
a universal differential field of J-". For any element e € e^'^^ is used to denote the set 
{eW : i = 0,...,k}. 

The sparse differential resultant is closely related with Laurent differential polynomials, 
which will be defined below. 

Definition 3.1 A Laurent differential monomial of order s is a Laurent monomial in vari- 
ables yW = {yl''^)i<i<n;0<k<s- More precisely, it has the form OLi nfc=o(yi''V'''' where 
dik are integers which can be negative. A Laurent differential polynomial is a finite linear 
combination of Laurent differential monomials with coefficients from £. 
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Clearly, the collections of all Laurent differential polynomials form a commutative dif- 
ferential ring under the obvious sum, product operations and the usual derivation operator 
S, where all Laurent differential monomials are invertible. We denote the differential ring of 
Laurent differential polynomials with coefficients in T by J-^{yi,yi^, . . . ,yn,yn^}, or simply 
by^{Y,Y-i}. 

Remark 3.2 J'{Y,Y-i} = T{yi ,yi^ , ■ ■ ■ ,ynjyn^} is only a notation for Laurent differen- 
tial polynomial ring. It is not equal to T[yj'^\ (y^"^)^^^ : k > 0]. 

Denote <S to be the set of all differential ideals in J^{Y, Y~^}, which are finitely generated. 
Let m be the set of all differential monomials in Y and T the set of all differential ideals in 
J^{Y}, each of which has the form 

[/i, . . . , M : m = {/ G T{Y}\ 3M G m, s.t. M ■ / G [/i, . . . , /.]} 

for arbitrary /j G J"{Y}. Now we give a one-to-one correspondence between <S and T. 
The maps : S — >■ T and ijj : T — > S are defined as follows: 

• Given any I = [Fi, . . . , Fg] G 5. Since each Fi G -F{Y, Y~^}, we can choose a vector 

(Ml, ...,Ms) G such that MjFj G T'jY} (i = 1, . . . , s). We then define (p(I) = 
[MiFi,...,MsFs] : inc7-{Y}. 

• For any J = [h, ...,/,]: m G T, define iP{J) = [/i, . . . , /,] in ^{Y, Y"!}. 

Lemma 3.3 The above maps (j) and are well defined. Moreover, (poip = idj- and ip o (j) = 
ids. 

Proof: tp is obviously well-defined. To show that is well-defined, it suffices to show that 
given another (iVi, . . . , A^^) G with NiFi G /"{Y} (i = 0, . . . , n), [MiFi, . . . , MsFg] : m = 
[NiFi, NsFg] : m follows. It follows directly from the fact that N^Fi G [MiFi, M^Fs] : 
m and MiFi G [NiFi, NsFs] : m. 

For any X = [Fi, . . . , F,] G S, V'o</.(X) = V([MiFi, . . . , M,F,] : m) = [M^Fi, M^F,] = 
X C J-"{Y, Y~^} where MiFi G J-'{Y}, since Laurent differential monomials are invertible. So 
we have o (f) = ids- And for any J = [/i, . . . , /^J : m G T, ^ o 'tp{J) = (l){[fi, . . . , fr]) = J. 
Thus, (f)oil) = idj- follows. □ 

From the above, for a finitely generated Laurent differential ideal X = [Fi, . . . ,^5], al- 
though (t){X) is unique, different vectors (Mi, . . . , Mg) G can be chosen to give different 
representations for (piX). Now the norm form for a Laurent differential polynomial is intro- 
duced to fix the choice of (Mi, . . . , Ms) G when we consider (/>(X). 

Definition 3.4 For every Laurent differential polynomial F G (J{Y, Y~^}, there exists a 
unique laurent differential monomial M such that 1) M ■ F £ £{Y} and 2) for any Laurent 
differential monomial T with T ■ F E ^^{Y}, T ■ F is divisible by M ■ F as differential 
polynomials. This M ■ F is defined to be the norm form of F, denoted by F^. The order of 
F^ is defined to be the effective order of F, denoted by Eord(F). Clearly, Eord(F) < ord(F). 
And the degree of F is defined to be the degree of F^, denoted by deg(F). 
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In the following, we consider zeros for Laurent differential polynomials. 

Definition 3.5 Let £^ = £\{a € £"[ B/c € N, s.t. a^'^) = 0}. Let F be a Laurent differential 
polynomial in J-{Y, Y^^}. A point (ai, . . . , a„) E (f"^)" is called a non-polynomial differential 
zero of F if F{ai, . . . ,a„) = 0. 

It becomes apparent why non-polynomial elements in £^ are considered as zeros of 
Laurent differential polynomials when defining the zero set of an ideal. If F € I, then 
{y^^^)~^F G I for any positive integer k, and in order for {y'^^^)~^F to be meaningful, we 

need to assume yf^^ ^ 0. We will see later in Example 13. 2 H how non-polynomial solutions 
are naturally related with the sparse differential resultant. 

3.2 Definition of sparse differential resultant 

In this section, the definition of the sparse differential resultant will be given. Since the study 
of sparse differential resultants becomes more transparent if we consider not individual dif- 
ferential polynomials but differential polynomials with indeterminate coefficients, the sparse 
differential resultant for Laurent differential polynomials with differential indeterminate co- 
efficients will be defined first. Then the sparse differential resultant for a given Laurent 
differential polynomial system with concrete coefficients is the value which the resultant in 
the generic case assumes for the given case. 

Suppose A = {Mio, Mil, . . . , J (i = 0, 1, . . . , n) where M^^ = ni=i Il/toCyf )'''^^' = 
is a Laurent differential monomial of order Sj with exponent vector aik € Z"^*'"*" ' 
and for ki ^ k2, Oiki 7^ aiA:2- Consider n + 1 generic Laurent differential polynomials defined 
over Ai {i = 0,1, ... ,n): 

h 

F, = Y,UikM,k{i = 0,...,n), (2) 

fc=0 

where all the Uik are differentially independent over Q. The set of exponent vectors Sj = 
{oik ■ k = 0, ... ,li} is called the support of Pj. The number |Sj| = /j + 1 is called the size of 
Pj. Note that Si is the order of Pj and an exponent vector of Pj contains n(sj + 1) elements. 
Denote 

Uj = {uio,Uii, . . . ,Uin) {i = 0,... ,n) and u = {uik : i = 0, . . . ,n;k = 1, . . . (3) 

To avoid the triviality, li > 1 (i = 0, . . . , n) are always assumed in this paper. 

Definition 3.6 A set of Laurent differential polynomials of form ^ is called a Laurent 
differentially essential system if there exist ki{i = 0, . . . , n) with 1 < ki < li such that 

d.tr.degQ(-^, . . . , ^^^^)/Q = n. In this case, we also say that AQ,...,An or So, • • • , 
Sn form a Laurent differentially essential system. 

Although Mjo are used as denominators to define differentially essential system, the 
following lemma shows that the definition does not depend on the choices of MjQ. 

Lemma 3.7 The following two conditions are equivalent. 
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1. There exist ko, . . . ,kn with 1 < ki < li such that d.tr.degQ( ^^^° , . . . , ^''^ ) / Q = n. 

2. There exist pairs {ki,ji) (i = 0, . . . , n) with ki 7^ € {0, . . . , k} such that 
d.tr.degQ(^, ...,§^)/Q = n. 

Proof: 1)^^ 2) is trivial. 

Now suppose 2) holds. Fix the n + 1 pairs {ki,ji), and without loss of generality, suppose 
-j^^, . . . , are differentially independent over Q. We need to show 1) holds. Suppose 

the contrary. Then we know that for any rrii € {1, . . . , -jj^, • • • , are differentially 

dependent over Q. Since = / -^^1 it follows that (i = 1, . . . , n) are differentially 
dependent over Q, which is a contradiction. □ 

Let [Po, . . . , Pn] be the differential ideal in Q{Y, Y"^; uq, . . . , u„} generated by Pj. By 
Lemma [331 [IPoj ■ ■ ■ > IPn] correspondents to [Pg'^, . . . , P^] : m C Q{Y; uq, . . . , u„} in a unique 
way. Moreover, we have the following lemma. 

Lemma 3.8 [Pq, . . . , P„] n Q{uo, . . . , u„} = ([P^,...,P^] : m) n Q{uo, . . . , u„}. 

Proof: Denote Pf = M^Pj (i = 0, . . . , n) where Mi are Laurent differential monomials. It 
is obvious that the right elimination ideal is contained in the left one. For the other direc- 
tion, let G be any element in the left ideal. Then there exist Hij € Q{Y,Y~^;uo,. . . ,u„} 

such that G = E^,H,,¥'f\ So G = Y.i,H,,{%-f^ = T.^,M^ff'' with G 
Q{Y, Y-i; Uq, . . . , u„}. Thus, there exists an M G m such that MG G [P^, . . . , P^] and 
G G ([P^, . . . ,P^] : m) n Q{uo, . . . ,u„} follows. □ 

In the whole paper, when talking about prime differential ideals, it is assumed that they 
are distinct from the unit differential ideal. The following result is the foundation for defining 
the sparse differential resultant. 

Theorem 3.9 Let Pq, . . . ,P„ he Laurent differential polynomials defined in Then the 
following assertions hold. 

• ([Pg^, . . . , P^] : m) is a prime differential ideal in Q{Y, uq, . . . , u„}. 

• ([Pq^, . . . , P^] : m) n Q{uo, . . . , u„} is of codimension 1 if and only if Pq, . . . , Pn form 
a Laurent differentially essential system. 

Proof: Let rj = (r/i, . . . , rjn) be a generic point of [0] over Q(u), where u is defined in ([3]). Let 

Then we claim that 6* = (r/i, ...,?]„; Co, ^oi, ••• , uqi^', ■ • • ; Cn, ""ni, • • • , UniJ is a generic point 
of ([P^, . . . , P^] : m), which follows that ([P^, Pf , . . . , P^] : m) is a prime differential ideal. 

Denote P^ = MjPj {i = 0, . . . , n) where where Mj are Laurent differential monomials. 
Clearly, Pf = M^Pi vanishes at 6 {i = 0,...,n). For any / G ([P^, Pf , . . . , P^] : m). 
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there exists an M G m such that Mf e [P^,Pf , . . . ,P^]. It follows that f{6) = 0. Con- 
versely, let / be any differential polynomial in Q{Y, uq, . . . , u„} satisfying f(6) = 0. Clearly, 
P^,PC . . . ,P^ constitute an autoreduced set with Uio as leaders. Let fi be the differential 
remainder of / w.r.t. this autoreduced set. Then /i is free from Ujo (i = 0, . . . , n) and there 
exist ki>0 such that niLo(^i^«o)*'" • / = /i, mod [P^, Pf , . . . , P^^]. Clearly, fi{9) = 0. 
Since /i G Q{u,Y}, /i = 0. Thus, / G [P^, Pf , . . . , P^] : m. So [P^, Pf , . . . , P^] : m is a 
prime differential ideal with 9 as its generic point. 

Consequently, ([P^,P^, . . . ,P^] : m) n Q{uo, . . . ,u„} is a prime differential ideal with 
a generic point C = (Co, ^^oi, • • • , uqi^^; ■ ■ ■ Cn,Uni, ■ ■ ■ ,UniJ. From it is clear that 
d.tr.degQ(C)/Q < Y17=o + ^- there exist pairs {ik,jk) {k = 1, . . . ,n) with 1 < jk < kj, 
and ifci 7^ ik2 (^i 7^ ^2) such that -j^^, • • • , -ff^ differentially independent over Q, 
then by Lemma 12. H Cin---;Ci„ a.re differentially independent over Q(u). It follows that 
d.tr.degQ(C)/Q = Er=o + Thus, ([P^, Pf , . . . , P^] : m) n Q{uo, . . . , u„} is of codimen- 
sion 1. 

Conversely, assume that ([P^,Pf,...,P^] : m) n Q{uo,...,u„} is of codimension 1. 
That is, d.tr.degQ(C)/Q = Y17=o^^ + ^- want to show that there exist pairs {ik,jk) 
{k = l,...,n) with 1 < jfc < and ik^ / ik2 {h / k2) such that ^wr^, • • • , are 

differentially independent over Q. Suppose the contrary, i.e., -jj^^^^, ■ ■ ■ , ^"''pff 
ferentially dependent for any n different ik and jk G {1, . . . ,lif.}. Since each Cj^. is a linear 
combination of {jk = ^,---ikk)^ it follows that Ciij---Xi„ are differentially de- 

pendent over Q(u). Thus, we have d.tr.degQ(C)/Q < Y17=o^i + a contradiction to the 
hypothesis. □ 

Combining Lemma 13.81 and Theorem 13.91 we have 

Corollary 3.10 [Po,Pi, . . . ,Pn] nQ{uo, . . . , u„} is a prime differential ideal of codimension 
one if and only if {Pj : i = 0, . . . ,n} is a Laurent differentially essential system. 

Now suppose {Pq, . . . , Pn} is a Laurent differentially essential system. Denote the differen- 
tial ideal [Pq, Pi, ... , P„] nQ{uo, . . . , u„} by I. Since I is of codimension one, by Lemma [2^ 
there exists an irreducible differential polynomial Il{u; uqq , . . . , Uno) = R(uo, . . . , u„) G 
Q{uo, . . . , u„} such that 

[Po, Pi, . . . , Pn] n Q{uo, . . . , u„} = sat(R) (5) 

where sat(R) is the saturation ideal of R. More explicitly, sat(R) is the whole set of differ- 
ential polynomials having zero differential remainders w.r.t. R under any ranking endowed 
on uo, . . . ,Un. 

Now the definition of sparse differential resultant is given as follows: 

Definition 3.11 R(uo, . . . , u„) G Q{uo, . . . , u„} in (5\) is defined to be the sparse differen- 
tial resultant of the Laurent differentially essential system Pq, . . . , Pn, denoted by Res^Q^...^^„ 
or RespQ^...^p„. And when all the Ai are equal to the same A, we simply denote it by Res^. 

From the proof of Theorem 13.91 and equation ([5]), R has the following useful property. 
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Corollary 3.12 Let R(uo, . . . , u„) he the sparse differential resultant of¥o, Pi, ... , P„. Then 
sat(R) C Q{uo, . . . , u„} is a prime differential ideal with a generic zero (u; (q, . . . , where 
Qi are defined in 

We give five examples which will be used throughout the paper. 
Example 3.13 Let n = 2 and Pj has the form 

Pj = UioVi + Uiiy'i + Ui2y2 [i = 0, 1, 2). 

It is easy to show that y'l' /y'l and y'^ jy'i clts- differentially independent over Q. Thus, 
Po,Pi,P2 form a Laurent differentially essential system. The sparse differential resultant 
is 

UOO Uoi U02 



R = Respg^; 



■uio un ui2 

U20 U2l U22 



Pay attention to the fact that R does not belong to the differential ideal generated by Pj in 
Q{Y; uq, . . . , u„} because each Pj is homogenous in yi,yi,y2 and R does not involve Y. 
That is why we use the ideal ([Po,Pi,P2] : m) C Q{Y; uq, . . . , u„} rather than [Po,Pi,P2] C 
Q{Y; uq, . . . , u^} in Theorem \3.9l Of course, R does belong to [Pq, . . . ,Pn] when regarded as 
a differential ideal of the Laurent differential polynomial ring Q{Y, Y^-*^; uq, . . . , u„}. 

The following example shows that for a Laurent differentially essential system, its sparse 
differential resultant may not involve the coefficients of some Pj. 

Example 3.14 Let n = 2 and Pj has the form 

Fo = ""00 + uoWiy'i, IPi = lilO + UlWl, ^2 = uio + uiiy'2. 

Clearly, Po,Pi,P2 form a Laurent differentially essential system. And the sparse differential 
resultant o/Po,Pi,P2 is 

R = uqiUio{uiiu\q - -uion'ii) + ^00^^11, 
which is free from the coefficients of P2 . 

Example 3.15 Let Aq = {i, 2/12/2}, Ai = {i, 2/1^/2} o^n-d A2 = {i, 2/12/2}- ^■^ ^'^^V 
verify that Ao,Ai,A2 form a Laurent differentially essential system. And ^es_Ao,Ai,A2 = 

U10U01U21U11UQQ - UioUooUiiU21Uq^ - uIiU21uIq - ■U01^^00«ll'"20- 

Example 3.16 Let n = 1 and Aq = Ai = {yf, (2/1)^, 2/1 2/1} • Clearly, Ao,Ai form a Laurent 
differentially essential system and Kes^ = ^ii^^oo ~ 2tioitiio'Uiinoo + ^oi'^io ~ 'Wi2^i02^'ii'Woo — 
U12U02U01U10 + UI2U01U00 + ■uioniino2- 

Example 3.17 Let n = 1 and Aq = Ai = {2/1 , 2/1 ) 2/i } • Clearly, Ao,Ai form a Laurent 
differentially essential system and Kes^ = —U12U01U00U10 — ui2Uq^u[q + 'Ui2^oi^i'Uoo + 

Ui2UoiUuUqq-UuUo2UooUio + UiiUq2u[qUoi+Uo2UoiuIq-uI^Uo2Uqq + UiiUo2UqiUio + UuuIqUi^ 
+^ll'"02'"00 - UuUq2UoiUio - niinoi-u'i2^00 + ^ioi^'l2'"lO - 'Wll^toi'"l2'UOO - u[iU02U0lUl0- 
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Remark 3.18 When all the Ai{i = 0, . . . , n) are sets of differential monomials, unless 
explicitly mentioned, we always consider Pj as Laurent differential polynomials. But when we 
regardVi as differential polynomials, Res^o_,,,^^„ is also called the sparse differential resultant 
of the differential polynomials Pj. In this paper, sometimes we regard Pj as differential 
polynomials where we will highlight it. 

We now define the sparse differential resultant for any set of specific Laurent differential 
polynomials over a Laurent differentially essential system. For any finite set A of Laurent 
differential monomials, denote by C{A) the set of Laurent differential polynomials of the 
form Y1iM£A'^^^^ where aM G £. Then C{A) can be considered as the affine space £^ or 
the projective space P(/ — 1) over £ where / = 

Definition 3.19 Let Ai = {Mjo, Mji, . . . , Mj;.} (z = 0,1,..., n) he finite sets of Laurent 
differential monomials which form a Laurent differentially essential system. Consider n + 1 
Laurent differential polynomials (Fq, Fi, . . . , Fn) S JliLo ^(•^i)- sparse differential resul- 
tant of Fq, Fi, . . . , Fn, denoted as ResFg^...^i?„, is obtained by replacing Uj by the corresponding 
coefficient vector of Fi in Res(uo, . . . , u„) which is the sparse differential resultant of the n+1 
generic Laurent differential polynomials in 

We will show in the next Section [3.31 that the sparse differential resultant ^esFo,...,F„ = 
will approximately measure whether or not the the over-determined equation system Fi = 
= 0, . . . , re) have a common non-polynomial solution. 

3.3 Necessary and sufficient condition for existence of non-polynomial so- 
lutions 

In the algebraic case, the resultant gives a necessary and sufficient condition for a system of 
homogenous polynomials to have common solutions. We will show that this is also true for 
sparse differential resultants in certain sense. 

To be more precise, we first introduce some notations. Let Aq,... ,An be a Laurent 
differentially essential system of monomial sets. Each element {Fq, . . . , Fn) € J~-{Ao) x • • • x 
C{An) can be represented by one and only one point (vo,...,v„) G £^0+^ x ••• x 
where v, = {vio,Vii, . . . ,Vii-) is the coefficient vector of F^. Let Zo{Ao, . . . , An) be the 
subset of f'o"*"^ X • • • X £^'"+1 consisting of points (vq, . . . , v„) such that the corresponding 
Fi = {i = 0, . . . ,n) have non-polynomial common solutions. That is, 

Zo{Ao, ...,An) = {(vo, . . . , v„) G X • • • X ^'"+1 : Fo = • • • = = have 

a common non-polynomial solution m{£^)^}. (6) 

The following result shows that the vanishing of sparse differential resultant gives a necessary 
condition for the existence of non-polynomial solutions. 

Lemma 3.20 ZoiAo, . . . , An) C V(sat(Res^o,...,^„)) . 

^Here, we can also consider the differential projective space P(li) over £ 
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Proof: Let Pq, . . . ,Pra be a generic Laurent differentially essential system corresponding to 
^0; • ■ ■ ) -^n with coefficient vectors uq, . . . , u„. By ([5]), [Fq, Pi, ... , P„] n Q{uo, . . . , u„} = 
sat(Res^o,...,^„). For any point (vq, . . . , v„) € Zo{Ao, ■ ■ ■ ,An), let (Pq, . . . ,Pn) G C{Ao) x 
• • • X C{An) be the differential polynomial system represented by (vq, . . . , v„). Let G be any 
differential polynomial in sat(Res^o^...^^„). Then G(vo, . . . , v„) € [Pq, • • • ,Pn] C £^{Y, Y~-^}. 
Since Pq, • • • , Pn have a non-polynomial common zero, G(vo, . . . , v„) should be zero. Thus, 
sat(Res^(,,...,^„) vanishes at (vq, . . . , v.„). □ 

Example 3.21 Continue from Example \3.13\. Suppose T = Q(x) and 5 = In this 
example, we have Resp(,^p^^P2 ^ 0. But yi = cux + cio, y2 = 0222;^ + C21X + C20 consist of a 
non-zero solution 0/ Pq = Pi = P2 = where Cij are distinct arbitrary constants. This shows 
that Lemma \3.2(A is not correct if we do not consider non-polynomial solutions. This example 
also shows why we need to consider non-polynomial differential solutions, or equivalently why 
we consider Laurent differential polynomials instead of usual differential polynomials. 

Let Zo{Ao, . . . , An) be the Kolchin differential closure of Zq^Aq, • • • , An) in £^'0+^ x • • • x 
£in+i_ Then we have the following theorem which gives another characterization for the 
sparse differential resultant. 

Theorem 3.22 Suppose the Laurent differential monomial sets Ai{i = 0, . . . , n) form a 
Laurent differentially essential system. Then Z{Aq, . . . .,An) = V(sat(ReSyiQ^...^^„)) . 



Proof: Firstly, by Lemma [3.201 Zq{Ao, . . . ,An) ^ V(sat(ReSy^o^...^^^)) . So Z{Ao, . . . ,An) = 
Zo{Ao, An) C V(sat(Res^o,...,^„)). 



For the other direction, follow the notations in the proof of Theorem 13.91 By Theo- 
rem [331 [1^0^ ) • • • ! IPn^ ] : m is a prime differential ideal with a generic point (rjX) where 
rj = {r]i, . . . ,r]n) is a generic point of [0] over Q{{uik)i=o,...,n;k^o) and C = (Co, lioi, • • • , 
uoig] . . . ; Cn, Uni, UniJ. Let {Fq, ...,Fn) G C{Ao) X • • • X C{An) be a set of Laurent 
differential polynomials represented by Clearly, r/ is a non-polynomial solution of Fi = 0. 
Thus, C £ ^o(-4o) • • • 5 -^n) C Z{Ao, . . . , An). By Corollary 13.121 ^ is a generic point of 
sat(Resyio^..._^„). It follows that V(sat(Res^Q^..._^„)) C Z{Ao, . . . , An). As a consequence. 



The above theorem shows that the sparse differential resultant gives a sufficient and 
necessary condition for a differentially essential system to have non-polynomial solutions 
over an open set of Y\7=o ^i-^i) in the sense of Kolchin topology. 

With Theorem 13.221 property 3) of Theorem 11.21 is proved. 

4 Criterion for Laurent differentially essential system in terms 
of supports 

Let Ai {i = 0, . . . ,n) he finite sets of Laurent differential monomials. According to Defini- 
tion l3.6l in order to check whether they form a Laurent differentially essential system, we need 
to check whether there exist Mj^. ,Mjj. € Ai{i = 0, . . . ,n) such that d.tr.degQ(MoA:(,/Mojo, 
. . . , Mnk„/Mnj„) /Q. = n. This can be done with the differential characteristic set method 




□ 
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via symbolic computation [45 j [3l H8| I16j . In this section, a criterion will be given to check 
whether a Laurent differential system is essential in terms of their supports, which is con- 
ceptually and computationally simpler than the naive approach based on the characteristic 
set method. 



4.1 Sets of Laurent differential monomials in reduced and T-shape forms 

In this section, two types of Laurent differential monomial sets are introduced, whose differ- 
ential transcendence degrees are easy to compute. 

Let Bi,B2, . . . , Bm be m Laurent differential monomials, where Bi = YYj=i Y[k>o(yj^^)'^^^'' ■ 
For each j € {1, . . . , n}, let qj = max^^oid{Bi, yj) . Let xi, . . . ,Xn be new algebraic indeter- 
minates and 

ij 

dij = '^dijkXj {i = 1,... ,mj = 1,... ,n) 

k=0 

univariate polynomials in Z[xj] respectively. If ord(i3j,yj) = — oo, then dij = and we 
denote deg{dij,Xj) = —oo. The vector {dii,di2, ■ ■ ■ ,din) is called the symbolic support vector 
of Bi. The following m x n matrix 



M 



( d\\ dyi ... d\n \ 
d21 d22 • • • d2n 

\ d'lril dm2 ■ ■ ■ dffin / 



is called the symbolic support matrix of i?i , . . . , B^ • 

Note that there is a one-to-one correspondence between Laurent differential monomials 
and their symbolic support vectors, so we will not distinguish these two concepts if there is 
no confusion. The same is true for a set of Laurent differential monomials and its symbolic 
support matrix. 

Definition 4.1 A set of Laurent differential monomials Bi,B2, . . . , i?m or its symbolic sup- 
port matrix M is called reduced if for each i < min(m, n), — oo ^ ord(i?j, yj) > ord{Bi^i:,yi), 
or equivalently — oo ^ deg{dii,Xi) > deg{di^k,i,Xi), holds for all k > 0. 

Note that a reduced symbolic support matrix is always of full rank since the term 
■|-|-min(m,n) ^ord{Bi,yi) ^.^j appear effectively in the determinant of the min(m, n)-th principal 
minor when expanded. 

Example 4.2 Let Bi = yly'ly'2,B2 = yl{y'2?yz{y'^? ,B^ = y[y'^. Then qi = 2,q2 = 1,^3 = 
1, and 

/ xl + 2 X2 
M = 3 2x2 2x3 -M 

\ Xi X3 

is reduced. 

Before giving the property of reduced symbolic support matrices, the following simple 
result about the differential transcendence degree will be proved. 
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Lemma 4.3 For r]i,r]2 in an extension field ofQ, d.tr.degQ(r/"^, ?7"^r72)/Q = d.tr.degQ(?7i, 

112) /Q, where 01,02 non-zero rational numbers. 

Proof: For any non-zero integer p, we have 

d.tr.degQ(r/i,r72)/Q = d.tr.degQ(?7i, ??2)/Q(???, r/2) + d.tr.degQ(r7?, r72)/Q 
= d.tr.degQ(?7f,r?2)/Q. 

So for each a G Q\{0}, d.tr.degQ(r?f , ?72)/Q = d.tr.degQ(r/i, r/2)/Q. Let ai=pi/qi {i = 1,2) 
wheie Pi, Qi are non-zero integers. Then, 

d.iT.degQ{vt\vTm)/Q = d.tr.degQ(r/l/''^^7?^/«^r72)/Q 

= d.tr.degQ(r/l'/^^^72)/Q (for Q(r/^/^^ ryf^''^) = Qivl^'^m)) 
= d.tr.degQ(r/i,r/2)/Q. 

□ 



Theorem 4.4 Let Bi, B2, ■ ■ ■ , Bj^ he m reduced Laurent differential monomials in Y. Then 
d.tr.degQ(Bi, B2,...,Bm)/Q = min(m, n). 



Proof: It suffices to prove the case m = n by the following two facts. In the case m > n, we 
need only to prove that Si, ... , B^ are differentially independent. And in the case m < ra, we 
can treat i/m+ii • • • j yn as parameters, then Bi,B2, ■ ■ ■ , B^ are still reduced Laurent differen- 
tial monomials. So if we have proved the result for m = n, d.tr.deg Q(i?i, i?2, • • • , Bm)/Q > 
d.tr.degQ(y^+i, . . .,yn){Bi,B2, . . . ,Bm)/Q{ym+i, ■■■,yn) =m follows. 

Since Bi,B2, ■ ■ ■ , Bn are reduced, wc have Oj = ord(i?j, yi) > for i < n. In this proof, a 
Laurent differential monomial will be treated as an algebraic Laurent monomial, or simply 
a monomial. Furthermore, the lex order between two monomials induced by the following 
variable order will be used. 



> 
> 
> 

> 
> 

> 



yi>y[>--->y{'"-'^ 



yn>yn>---> > yt^ > > ••• 



t/n-l ^ Vn-l ^ 



(01) (oi + l) 

y\ ^ >y\ > 



Under this ordering, we claim that the leading monomial of S^Bi{l < i < n,t e N) is 



LMit = 



Bi*y, 



(Of + t) 



7^ 



-. Here by leading monomial, we mean the monomial with the highest 



order appearing effectively in a polynomial. Let Bi 



N.ijj'f'^)^^ {I <i < n). If Ni 
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1, then the monomials of 6^Bi is of the form Y\k=oiyi^ y^^ where so,...,st are non- 
negative integers such that X]fc=o = Di and Yl\=i = i- Among these monomials, if 
Sk > for some I < k < t — 1, then sq is strictly less than Di — 1 and Yl^^oiVi"''^''^)^'' < 
iVi )^'~^yi = s — follows. Hence, in the case Ni = 1, the claim holds. Now 

Vi 

suppose Ni 7^ 1, then it is a product of variables with lex order larger than yl°^\ Then 
6*Bi = X;Lo {i)^''Ni6*~''{yi"''')^'- If A; = 0, then similar to the case Ni = 1, we can show 
that the highest monomial in NiS^y^"'^)^^ is A^i(2/5°'V'~^2/?'^*^- For each A; > 0, 6''Ni < Ni 

and 6''Ni6^-''{yl°'^)^' < Ni{yf'^)^^-^yf'^^^ = . Hence, the leading monomial of 

y^ " 

5^B, is iV.(y^V-^2/f = ^^4?^. 

y^ 

We claim that these leading monomials LMn = ' /p.) — {i = 1, . . . , m; t > 0) are alge- 

vi 

braically independent over Q. We prove this claim by showing that the algebraic transcen- 
dence degree of these monomials are the same as the number of monomials for any fixed t. Let 
= [y,,y[,...,y\''^-% Y* = [i/f^^'+^,...,^^^\ Bu = [B„ LMa, . . . , LMu] for 1 < i < 

n. We denote by BYi = (y^^^Y^ BYu = [{y^Y^ {yi"'Y'''yt^'\ {yl'^Y^^'yt'''^] 
for 1 < i < n. Then, by Lemma 14.31 we have 

n{t + l) > tT.degq{Bu,B2t,...,Bnt)/Q 
> tr.degQi(Sii,S2t,...,S„i)/Qi 
= tr.deg QiiBYu, BY2t, BYnt) /Qi 
= n{t + l) 

where Qi = Q(li, . . . ,Yn,Y* , . . . , Y*). Hence, this claim is proved. 

Now, we prove that Bi,...,Bn are differentially independent over Q. Suppose the 
contrary, then there exists a nonzero differential polynomial P € Q{zi,...,Zn} such that 
P{Bi, . . . , Bn) = 0. Let P = CfcPfc, where P^ is a monomial and G Q\{0}. Then, the 
leading monomial of Pk{Bi, . . . , i?„) is a product of LMa {i = 1, . . . ,n;t > 0). We denote 
this product by LMP^, then LMP^ ^ LMPj for k ^ j since these LMn are algebraically 
independent. But there exists one and only one product which has the highest order, which 
can not be eliminated by the others, which means that P{Bi, . . . , Bn) ^ 0, a contradiction. 

□ 

In general, we cannot reduce a symbolic support matrix to a reduced one. But, in the 
next section, we will show that any symbolic support matrix can be reduced to T-shape to 
be defined below. 

Definition 4.5 A set of Laurent differential monomials Bi, . . . , Bm or their symbolic sup- 
port matrix M is said to he in T-shape with index if there exist 1 < i < min(m, n), < 
j < min(m, n)—i such that all elements except those in the first i rows and the i+1, . . . , (i+j)- 
th columns of M are zeros and the sub-matrix consisting of the first i + j columns of M is 
reduced. The zero sub-matrix (^1,^2) in FigureUlis called the zero sub-matrix of M . 

In Figure [H an illustrative form of a matrix in T-shape is given, where the sub-matrices 
Ml and M2 of the matrix are reduced ones. It is easy to see that Mi must be an z x i square 
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Ml 






Zi 


M2 


Z2 



Figure 1: A T-shape Matrix 



matrix. Since the first i + j columns of a T-shape matrix M is a reduced sub-matrix, we 
have 

Lemma 4.6 The rank of a T-shape matrix with index {i,j) equals to i+ j. Furthermore, a 
T-shape matrix is reduced if and only if it is of full rank, that is, i + j = min (m, n). 

For a zero matrix S with k rows and / columns whose elements are zeros, we define its 
0-rank to be A; + /. A T-shape matrix M is not of full rank if and only if i + j < min(m, n). 
As a consequence, we have 

Lemma 4.7 A T-shape matrix of index {i,j) is not of full rank if and only if its zero sub- 
matrix is an (m — i) x (n — j) zero matrix with 0-rank m + n — i — j > max(m, n) + 1. 

The differential transcendence degree of m Laurent differential monomials in T-shape 
can be easily determined, as shown by the following result. 

Theorem 4.8 Let Bi,. . . , -Bm be m Laurent differential monomials and M their symbolic 
support matrix which is in T-shape with index {i,j)- Then d.tr.deg Q{Bi,B2, ■ ■ ■ ,Bm)/'^ = 
rk(M) =i+ j. 

Proof: Since M is a T-shape matrix with index («, j), by Lemma 14.61 the rank of M is i + j. 

Deleting the zero columns of the symbolic support matrix of -Bi+i, . . . ,Bm, we can get 
a reduced matrix. By Theorem 14.41 we have d.tr.deg Q(i?i+i, i?m)/Q = j- Since the 
symbolic support matrix of Bi,...,Bi is also a reduced one, by Theorem 14.41 we have 
d.tr.deg Q(Bi, ... ,Bi)/Q = z. Hence, 

d.tr.deg Q(5i,...,B^)/Q = d.tr.deg Q(5i, S^)/Q(5,+i, 5^) 

+d.tr.degQ(Si+i,...,5^)/Q 
< d.tr.deg Q(5i,...,5,)/Q + j 
= i + j- 

On the other hand, if we treat yi+i, . . . ,yi+j and their derivatives as parameters, the 
symbolic support matrix oi Bi, . . . , Bi is also a reduced one and the rank of this matrix 
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is i. By Theorem 1131 we have d.tT.degQ{yi+i, . . . ,yi+j){Bi, . . . , Bi) /Q{yi+i, . . . ,yi+j) = 
i. Since Sj+i, . . . ,Bm are monomials in yi+i, . . . , y-i+j (see Figure d]), Q(-Bj+i, . . . , Bm) C 
i+i, • • • -.Vi+j)- Hence, 



d.tr.degQ(Si,...,5™)/Q = d.tr.degQ(Bi, . . . , 5™)/Q(S,+i, . . . , 5™) 

+d.tr.degQ(Si+i,...,5^)/Q 
> d.tr.degQ(?/i+i, . . .,yi+j){Bi, . . . , Bi) /Q{yi+i, . . .,yi+j) +j 
= i + j- 

Thus, d.tr.degQ(5i, . . .,B^m)/Q = rk(M) = i + j. □ 



4.2 An algorithm to reduce Laurent differential monomials to T-shape 

In this section, an algorithm is given to reduce any set of Laurent differential monomials to 
a set of Laurent differential monomials in T-shape, which has the same differential transcen- 
dence degree with the original one. 

First, we will define the transformations that will be used to reduce any symbolic support 
matrix to a T-shape one. A Q-elementary transformation for a matrix M consists of two 
types of matrix row operations and one type of matrix column operations. To be more 
precise, Type 1 operations consist of interchanging two rows of M; Type 2 operations consist 
of adding a rational number multiple of one row to another; and Type 3 operations consist 
of interchanging two columns. 

Let Bi, . . . , Bm be Laurent differential monomials and M their symbolic support ma- 
trix. Then Q-elementary transformations of M correspond to certain transformations of the 
monomials. Indeed, interchanging the i-th and the j-th rows of M means interchanging Bi 
and Bj , and interchanging the i-th and the j-th columns of M means interchanging yi and yj 
in . . . , Bm{or in the variable order). Multiplying the i-th row of M by a rational number 
r and adding the result to the j-th row means changing Bj to B'-Bj. 

Lemma 4.9 Let Bi, . . . , Bm be Laurent differential monomials and Ci, . . . , Cm obtained by a 
series of Q-elementary transformations from Bi, . . . , Bm- Then d.tr.degQ(-Bi, . . . , Bm)/Q = 
d.tr.degQ(C7i, ...,C„)/Q. 

Proof: It is a direct consequence of Lemma 14.31 □ 

Now, an algorithm RDM(M) will be given to reduce a given symbolic support matrix 
to a T-shape matrix by a series of Q-elementary transformations. We sketch the algorithm 
below. Note that we still denote by M the matrix obtained by Q-elementary transformations 
from M. We assume that m < n and hence p = max(m, n) = n. The case m > n can be 
shown similarly. 

Let be a sub-matrix of M. Then the complementary matrix oi N in M is the sub- 
matrix of M from which all the rows and columns associated with N have been removed. 

The algorithm consists of three major steps. In the first step, a procedure similar to the 
Gauss elimination will be used to construct a reduced square sub-matrix R oi M such that 
the complementary matrix of ii in M is a zero matrix. Precisely, choose a column of M, say 
the first column, which contains at least one non-zero element. Then, choose an element. 
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say dii, of this column, which has the largest degree among all elements in the same column. 



ai 



If there exists a dii,i > 1 such that deg(dji) = deg(dii), then replace dij by di 
for j = l,...,n, where aj,ai are the leading coefficients of dii,dii respectively. This is a 
Q-elementary transformation of Type 2. Repeat the above procedure until the first column 
is in reduced form, that is deg(dji) < deg(dii) for i = 2, ... ,n. Consider the lower-right 
(m — 1) X (n — 1) sub-matrix of M and repeat the above procedure for N. In this way, 
we will obtain a reduced square matrix whose complementary matrix is a zero matrix Z in 
the lower-right corner of M. 

In the second step, a recursive procedure is used to construct a reduced form of M. Let 
the zero matrix Z obtained above be an i x j matrix. Denote r = i + j to he the 0-rank of 
it. If j = n, the last i rows of M are zero rows. Delete the last i rows from M, then we have 
a strictly smaller matrix, which can be treated recursively. 

If r > n + 1, M cannot be of full row rank, which will be considered in step three. 
Otherwise, let Mc be the lower-right (m+r— max(m, n)) x (n+r— max(m, n)) = (m+r—n)xr 
sub-matrix of M, Mci the lower-left i x [n + i — max(m, n)) = i x i sub-matrix of Mc, and 
M(72 the upper-right (m + j — max(m,n)) x j = [m + j — n) x j sub-matrix of Mc- In 
Figure [2l|a,b), Mc is represented by the pink area. Here, Mc is chosen to be the minimal 
(m — q) X (n — q) sub-matrix of M at the lower-right corner, which may have full rank. 
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(a) Sub-matrices for m < n (6) Sub-matrices for m > n 



(c) Two zero sub-matrices 
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(d) N3 and (e) Compute (/) Interchanging rows 

Figure 2: Matrix forms in Algorithm [H the blue parts are reduced ones 

Let Ni =RDM(Mci) and N2 =RDM(Mc2). Note that the Q-elementary transforma- 
tions of these sub-procedures are for the whole rows and columns of M. By doing so, the 
sub-matrix consisting of the first n — r columns of M remains to be a reduced one. 

If A^i and N2 are reduced matrices, we can obtain a reduced matrix for M by a suitable 
column interchanging. Otherwise, either A^i or is not of full rank. Assume A'^i is not of 
full rank. Then merging the zero sub-matrix of A^i and Z, we obtain a zero matrix with 
0-rank larger than that of Z (Figure [2]^c)). Repeat the second step for M with this new zero 
sub-matrix. 
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In the third step, M contams a "large" zero sub-matrix and a T-shape matrix of M can 
be constructed directly as follows. Let the zero matrix Z at the lower-right corner of M be 
an i X J matrix and r = i + j. Let Mc3 be the lower- left i x {n — j) sub-matrix of M and 
N2, = RDM(Afc3). In this case, Mc^ has more rows than columns. We can assume that 
A^3 is of full column rank. Otherwise, a sub-matrix of A''3 can be used as A''3. 

Let Mci be the upper-right (m — i) x j sub-matrix of M, A^4 = RDM(Mc4), and s = 
rk(A'^4) (see Figure [2jd)). If A'4 is of full row rank, then by suitable column interchangings, 
we can obtain a T-shape matrix. Otherwise, let the lower-left (m — s) x (n — j) sub-matrix 
of M be Mc5, and = RDM(Mc5), which is a reduced matrix with full column rank, see 
Figure [2)[e) . Now, by suitable column interchangings, we can obtain a T-shape matrix (see 
Figure EJf)). 

The idea of the algorithm is as follows. Try to use the first step to construct a reduced 
matrix. If the first step fails to do so, use the second step to change the matrix so that it 
contains a larger zero sub-matrix after each iteration. The procedure will end until either 
a T-shape matrix is obtained or the matrix has a zero sub-matrix with size larger than 
max(m,?i) + 1, in which case a T-shape matrix can be obtained directly. 

We now use the following example to illustrate the first two steps of the algorithm. 

Example 4.10 Let = yM'ysv's, = yliy'ify'^fvliy'^f , ^3 = yliy'if ^^(2/2")' 
2/3(7/3)'^ . Then, the symbolic support matrix is 



M = I 2x1+3 2x1 + 

3x^ + X2 






We will use this matrix to illustrate the algorithm. 

2 

X2 

The matrix after =^ is obtained with the first step of the algorithm. We first use du = xi + 1 
to reduce the degrees of 2xi + 2 and 3xi + 3 with Q-elementary transformations of Type 2. 
Since x\ is of greater degree than X2, nothing needs to do. Finally, we obtain a 1 x 1 zero 
matrix at the lower-right corner at the end of step 1. 

Now, goto the second step of the algorithm. We have r = 2 < max(m, n) + 1 = 4. Mc is 
the lower-right 2x2 sub-matrix of M , Mqi = (3^2); o,nd Mc2 = (!)■ 

Since both Mci and Mc2 o-'^^ reduced, we interchange the second and third columns of 

M to obtain the final matrix after =>, which is reduced. The corresponding monomials are 
Di = 2/12/1 y2"2/3y3' B>2 = 2/12/2 2/3' '^''^^ ^3 — v'llv^- T-shape under the variable order 

y\>yi> 2/2- 

We use the following example to illustrate the third step of the algorithm. 
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Example 4.11 Let Bi 

yl 



B2 



B3 



^ Then, the symbolic support matrix is M given below. 



y'my's' 



y'l, 



B5 = 



M = 



( 4 
xl 

X\ 

Xl 

2 



xl 
X2 






xl + X3 






2 \ 

2 



0/ 



X 



q + X 

Xl 
Xl 

2 



id) 



X 



V 



-xf + xf 

Xl 


2 











X3 

X3 + 1 
-X3 - 1 




2\ 





0/ 



v 



X-? 
' 




X3 



X3 + I 
-X3 - 1 
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X3 + I 
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xl 
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—xl + xl 








Xl 

















2 





0) 



For step 1 of the algorithm, we do nothing to M and the zero matrix Z obtained at the end 
of this step is a 2 X 2 zero sub-matrix at the lower-right corner of M. In step 2, Mc is set 

to be the lower-right 4x4 sub-matrix of M , Mqi = ( ^ n V '^'^^ 



In step 3, we have Mc3 



^ y \^ 

Merging Z and Mci, we obtain a 2 x 4 zero sub-matrix at the lower-right comer of M. 
Up to now, M is not changed. Then, step 3 of the algorithm is applied. 

Xl 

2 

1 2 



Since Mc3 is reduced and of full rank, we execute 



X3 



case 1 by setting Mqa = 



X2 X3 + X3 
X3 + I 




and N4 = IlDMfMci) which is a T- 



shape matrix with index (1,1) and is not of full rank. Now, M becomes the matrix after 
which contains N4. Since N4 is not of full rank, let Mc5 = {xi, xi, 2)-^ and compute = 

RDM(^Mc5j. Now M becomes the matrix after We interchange the first column and 

the 2, 3-th columns of M to obtain the final matrix which is in T-shape with index (1,2). 

The corresponding monomials are Di = y'ly'-^'ysyiy'i, D2 = y'ly3/yi, D3 = 2/12/32/3; 
-D4 = 1/(2/32/3), -D5 = 2/i- It is of T-shape under the variable order 2/2 > 2/3 > 2/i > 2/4 > 2/5- 
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Algorithm 1 — RDM(M) 

Input: Laurent differential monomials Bi, . . . , Bm in Y or their symbolic support matrix 

M = {dij)i<i<rn,l<j<n- 

Output: A T-shape matrix which is obtained from M by Q-elementary transformations. 
Initial: Let s = 1, p = max(m,n). 

1. While s < min(m, n) do 

1.1 If for any j, / > s, deg{dji) = — oo, letz = m — s + l,j = n — s + 1 and go to Step 2. 

1.2 Select j,l > s such that — oo ^ deg{dji) > deg{dii) for any i > s. Interchange the j-th 
row and the s-th row, the l-th column and the s-th column of M. Using dgs to do 
Q-elementary transformations such that deg{dss) > deg{dis) for i > s. 

1.3 s = s + l. 

2. Let r = i + j be the 0-rank of the i x j zero sub-matrix in the lower-right side of M. 

2.1 If M is already a T-shape matrix, return M. 

2.2 If j = n, delete the last i rows from M, and let A=RDM(M). 
Then add i rows of zeros at the bottom of N and return this matrix. 

2.3 If r >p-M, go to Step 3. 

2.4 Let Mc be the lower-right {m + r — p) x (n + r — p) sub-matrix of M. 
Let the lower-left i x {n + i — p) sub-matrix of Mc be Mci and the 
upper-right {m + j — p) x j sub-matrix of Mc be Mc2- (see (a, b) of Fig. [2]) 

2.5 Let iVi = RDM(Mci) and N2 = RDM(Mc2). 

2.6 If Ni,N2 are reduced matrices, interchange the p — r + 1 ton — j columns and the 
n — j + 1 to n + m — p columns of M, return the obtained T-shape (reduced) matrix. 

2.7 If the k x / zero sub-matrix Zi of Ni has 0-rank k + l > max(i, n + i — p) + 1 = i-|-l, 
combine Zi and the ? x j zero matrix to obtain a /c x {I + j) zero matrix with 0-rank 
k + 1 + j > i + j (see (c) of Fig. [2]). Let i = k,j = I + j, go to Step 2. 

2.8 Else, the kx I zero sub-matrix Z2 of N2 has 0-rank k + l > ma.x{m+j—p,j) + l = j + 1, 
combine Z2 and the i x j zero matrix to obtain a {k + i) x I zero matrix with 0-rank 
k + l + i>i+j. Let i = k + i, j = I, go to Step 2. 

3. Let Mc3 be the lower-left i x {n — j) sub-matrix of M and A3 = RDM(Mc3) with index 
{k,l). 

3.1 If / = 0, delete the last i - k rows from M, let A=RDM(Af), add i - A; zero 

rows at the bottom of N and return this matrix. 

3.2 If A3 is not of full rank, put the A: + 1, . . . , (A; + /)-th columns as the first / columns 
of M. Let i = i — k, j = n — I, A3 the lower-left i x {n — j) sub-matrix of M. 

3.3 Now A3 is of full column rank. Let the upper-right (m — i) x j sub-matrix of M 
be Mc4, A4 = RDM(Mc4), and s = rk(A4). 

3.4 Let the lower-left {m-s)x {n-j) sub-matrix of M be Mc5 and A5 = RDM(Mc5). 
Interchange the first n — j columns and the n — j + 1 to n — j + s columns of M, and 
return the obtained T-shape matrix. (See (d,e,f ) of Fig. [2j) 

In this algorithm, the Q-elementary transformations in RDM(Mcj) (i = 1, . . . , 5) are also 
for the whole m x n matrix. And in each step the new m x n matrix obtained after doing 
Q-elementary transformations is also denoted by M. 
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The following theorem proves what we claimed before. 

Theorem 4.12 The symbolic support matrix of any Laurent differential monomials Bi, . . . , 
Bffi can be reduced to a T-shape matrix by a finite number of Q^- elementary transformations. 

Proof: We assume that m < n and hence p = max(m, n) = n. The case m > n can be 
proved similarly. 

We prove the theorem by induction on the size of the matrix M, that is, m + n. One can 
easily verify that the claim is true when m + n = 2, 3, 4. Assume it holds for m + n < s — 1, 
we consider the case m + n = s. 

If a T-shape matrix is obtained in Step 1, then the theorem is proved. Otherwise, let Z 
be the i x j zero matrix obtained in this step. Since the complementary matrix of Z in M 
is a square matrix, the 0-rank of Z is larger than max(m, n) — min(m, n) + 1. 

In Step 2.2, M contains zero rows. By deleting these zero rows, the size of M is decreased 
by one at least. By induction, the algorithm is valid. 

In Step 2.3, from r > max(m, n) + 1, we have r = i + j >n + l and i > n — j. Then the 
i X (n — j) left-lower sub-matrix of M has more rows than columns. As a consequence, and 
M cannot be of full rank. 

In Step 2.4, Mq is chosen as the minimal sub-matrix of M such that it is of type (m — 
q) X {n — q) which may have full row rank. This implies that Mci must be an z x i square 
matrix, and hence q = n — r and Mc is an {m + r — n) xr matrix. Since the complementary 
matrix of Z in M is a square matrix, we have j>j— i = n — m. Hence m + r — n > i and 
Mc contains Z as a sub-matrix for the first loop, and this is always true since Z is from AIq 
and the size of Mq is increasing for each loop. 

In Step 2.5, by the induction hypothesis, A^i = RDM(Mci) and N2 = RDM(Afc'2) can 
be computed. Moreover, the lower-left mx {n — r) sub-matrix of M is always a reduced one 
although the Q-elementary transformations are for the whole rows of M. 

In Step 2.6, A'^i and are reduced with full rank. The algorithm terminates and returns 
a reduced matrix by suitable column interchanging given in the algorithm. 

In Step 2.7, A'^i is not of full rank. Then by Lemma 14.71 the k x / zero sub-matrix of A'^i 
has 0-rank k + 1 > max(i, n + i — max(m, n)) -|- 1 = i -|- 1. The i x j zero sub-matrix Z and 
this k X I zero sub-matrix form a k x {I + j) zero-matrix, with 0-rank k + j + l>i + j + l 
(Figure [2][|c)). Step 2.8 can be considered similarly. Since in each loop of Step 2, the 0-rank 
of the zero-matrix Z of M increases strictly, this loop will terminate. 

Step 3 treats the case when M is not of full rank. Note that Mas has more rows than 
columns. Step 3.1 is correct due to the induction hypothesis. 

For Step 3.2, since N3 is not of full rank and does not contain zero rows, we have / > 
and i > k. These conditions make the constructions given in the algorithm possible. 

In Step 3.3, A^3 is an i x {n — j) reduced matrix with full column rank and the lower-right 
ix j sub-matrix of M is a zero matrix. Due to this condition, the remaining steps are clearly 
valid. Also note that Mcs is obtained from by adding several more rows. Then Mcs is 
also of full column rank and hence is a reduced matrix of full column rank. □ 

Theorem 4.13 The differential transcendence degree of the Laurent differential monomials 
Bi, . . . , Bm over Q equals to the rank of their symbolic support matrix. 
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Proof: By Lemma 14.91 Q-elementary transformations keep the differential transcendence 
degree. The result follows from Theorems 14.81 and 14.121 □ 

Theorem 14.131 can be used to check whether the Laurent polynomial system ([2]) is differ- 
entially essential as shown by the following result. 

Corollary 4.14 The Laurent differential system (0j is Laurent differentially essential if and 
only if there exist Mij- (i = 0, . . . , n) with 1 < ji < h such that the symbolic support matrix 
of the Laurent differential monomials Moj^/Moq, . . . , Mnj„/Mno is of rank n. 

By Corollary 3.4 of [15j , the complexity to compute the determinant of a sub-matrix Mg 

2 

of M with size k x k is bounded by 0{k^^^L^''^ A), where L = log ||-/Vfs||, 7 denotes the 
number of arithmetic operations required for multiplying a scalar vector by the matrix Ms, 
and A is the degree bound of Ms- So, the complexity to compute the rank of M is single 
exponential at most. 

Remark 4.15 A practical way to check whether the Laurent differential system ^ is Lau- 
rent differentially essential is given below. 

• Choose n-\- 1 monomials Mij- (i = 0, . . . , n) with 1 < ji < h- 

• Use Algorithmic to reduce the symbolic support matrix of Moj^/Moq, . . . ,Mnj„/Mno to 
a T- shape matrix M . 

• Use Theorem \4.8\ to check whether the rank of AI is n. 

• If the rank of M is n, then the system is essential. Otherwise, we need to choose 
another set of n + 1 monomials and repeat the procedure. 

The number of possible choices for the n + 1 monomials is HILo which is very large. But, 
the procedure is efficient for two reasons. Firstly, Algorithm [1] is very efficient, since we are 
essentially doing numerical computation instead of symbolic ones. Secondly, the probability 
for n-\- 1 Laurent monomials to have differential transcendence degree n is very high. As a 
consequence, we do not need to repeat the procedure for many choices of n -|- 1 monomials. 

By Corollary 14.141 property 3) of Theorem 11.11 is proved. 

4.3 Differential transcendence degree of generic Laurent differential poly- 
nomials 

Algorithm [1] and Theorem 14.131 show how to reduce the computation of the differential 
transcendence degree of a set of Laurent differential monomials to the computation of the 
rank of their symbolic support matrix. In this section, this result will be extended to compute 
the differential transcendence degree of a set of generic Laurent differential polynomials. 
Consider m generic Laurent differential polynomials 

k 

Pi = UioMio + ^ UikMik (i = 1, . . . , m), (7) 

k=l 



28 



where all the Uik are differentially independent over Q. Let /Jjfc be the symbolic support 
vector of Mik/MiQ. Then the vector Wi = Yl\^QUik(iik is called the symbolic support vector 
of Pj, and the matrix Mp with wi, . . . , Wm as its rows is called the symbolic support matrix 
of Pi, . . . , Pm- Then, we have the following results. 

Lemma 4.16 Let Mki^...^km be the symbolic support matrix of the Laurent differential mono- 
mials (Mifcj/Mio, . . . ,Mmkm/Mmo)- Then rk(Mp) = maxi<fc^<i^rk(Mfc,,...,fc^). 

Proof: Let the rank of Mp be r. Without loss of generality, we assume that the r xr leading 
principal sub-matrix of Mp, say Mp^^, is of full rank. By the properties of determinant, 

h Ir 

det(Mp_r) = X] ■ ■ ■ S YVi=i^ikiT>iki, . . . , kr) where D{ki, . . . , kr) is the determinant of 

ki = l kr = l 

the r X r leading principal sub-matrix of Myfc^^,,, ^^. So det(Mp^r) 7^ if and only if there 
exist ki, . . . ,kr such that D{ki, . . . , kr) 7^ 0. Hence, the rank of is no less than the 

rank of Mp. On the other hand, let s = maxi<fc.<;.rk(Mfcj^ ,,^fc^). Without loss of generality, 
we assume D{ki, . . . , kg) 7^ 0, then, det(Mp_s) 7^ 0. Hence, s is no greater than the rank of 
Mp. □ 

The following result is interesting in that it reduces the computation of differential tran- 
scendence degree for a set of generic differential polynomials to the computation of the rank 
of a matrix, which is analogue to the similar result for linear equations. 

Theorem 4.17 d.tr.degQ(U^iUi)(Pi/Mio, . . . ,P.„/M^o>/Q(U[^iUi) = rk(Mp), where 

= {UiO, ■ ■ ■ ,Uil-). 

Proof: By Lemma l2.ll the differential transcendence degree of Pi/Mio, . . . ,Pm/Mmo is no 
less than the maximal differential transcendence degree of Mi/j^/Mio, . . . , Mmkm/^mO- 

On the other hand, the differential transcendence degree will not increase by linear com- 
binations since d.tr.degQ(A)(ai + Xdi,a2, ■ ■ ■ ,ak)/Q{X) < max(d.tr.degQ(oi, 02 . . . , afc)/Q, 
d.tr.degQ(ai, 02, . . . ,ak)/Q) for any differential polynomial ai{l < i < k) and di. So, the 
differential transcendence degree of Pi/Mio, . . . ,Pm/Mmo over Q{U'^^ni) is no greater than 
the maximal differential transcendence degree of Mifc^/Mio, . . . , M^fc^/Mmo- 

So, we have d.tr.degQ(U[^iUi)(Pi/Mio, . . . , P„/M^o)/Q(U^iUi) = maxfc,,„„fc„d.tr.deg 
Q(Mife^/Mio, . . . , Mmkm/Mmo)/Q- By Theorem gTS] and LemmaHUl the differential tran- 
scendence degree of Pi/Mio, . . . ,Pm/Mmo equals to the rank of Mp. □ 

By Theorem 14.171 we have the following criterion for system ^ to be differentially 
essential. 

Corollary 4.18 The Laurent differential system (0) is Laurent differentially essential if and 
only i/rk(Mp) = n. 

The difference between Corollary 14. 141 and Corollary l4.18l is that, in the later case we need 
only to compute the rank of a single matrix whose elements are multivariate polynomials 
in X^^o(^« + 1) + n variables, while in the former case we have to compute the ranks of 
nr=o ^« matrices whose elements are univariate polynomials in n separate variables. One 
also can replace Uj^- by v^' in Mp, where Vi is a new variable, then the elements of Mp will 
be multivariate polynomials in 2n -|- 1 variables. 
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In the rest of this section, properties for the ehmination ideal 



Tu = ([Pf , . . . , P^^^] : m) n Q{ui, . . . , u^} (8) 

will be studied, where Pj are defined in ([7|) and Uj = (ujo, • • • ,Uii-). These results will lead 
to deeper understandings for the sparse differential resultant. 

Theorem 4.19 Let Zu be defined in (B^. Then X„ is a differential prime ideal with codi- 
mension m — rk(Mp). 

Proof: Let rj = {rji, . . . ,r]n) be a generic point of [0] over Q(u), where u = {uik : i = 
1, . . . ,m;k = 1, . . . ,li} and 



k=l 



Similar to the proof of Theorem l3.9l we can show that 6 = (771, . . . , ry^; Ci^uu, . . . , uu-^ ; . . . ; (^m, 
Umi, ■ ■ ■ , Umim) is a generic point of [P^, . . . , P^] : m, which follows that [Pf^, . . . , P^] : m is 
a prime differential ideal in Q{Y, ui, . . . , Um}- As a consequence, X„ is a prime differential 
ideal. Since (^1, . . . , C,Yn are free of UiQ (i — 1, . . . , w-), by Theorem 14. 17^ 

d.tr.degQ(u)(Ci,...,Cm)/Q(u) 

= d.tr.degQ(ui, . . . , Um)(Ci, ■ ■ ■ , Cm)/Q(ui, . . . , u^) 

= d.tr.degQ(ui, . . . , u^)(^, . . . , |f^)/Q(ui, . . . , u^) 

= rk(Mp). 

Hence, the codimension of is m — rk(Mp). □ 

In the following, two applications of Theorem l4. 191 will be given. The first application is to 
identify certain Pj such that their coefficients will not occur in the sparse differential resultant. 
This will lead to simplification in the computation of the resultant. Let T C {0, 1, . . . 
We denote by Pt the Laurent differential polynomial set consisting of Pj (i € T), and Mp^ 
its symbolic support matrix. For a subset T C {0, 1, . . . ,n}, if card(T) = rk(Mp^), then Px, 
or {Ai : i G T}, is called a differentially independent set. 

Definition 4.20 LetT C {0, 1, . . . ,n}. Then we sayT or¥j is rank essential if the following 
conditions hold: (1) card(T) — rk(Mpj.) = 1 and (2) card(JI) = rk(Mpj) for each proper subset 
J o/T. 

Note that rank essential system is the differential analogue of essential system introduced 
in [20]. Using this definition, we have the following property, which is similar to Corollary 
1.1 in m. 



Theorem 4.21 If {Po,...,P„} is a Laurent differentially essential system, then for any 
T C {0, 1, . . . ,n}, card(T) — rk(Mp^) < 1 and there exists a unique T which is rank essential. 
In this case, the sparse differential resultant o/Po,...,Pn involves only the coefficients of 

P* (ieJ). 
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Proof: Since n = rk(Mp) < rk(Mp^)+card(P)-card(PT) = n+l+rk(Mp^)-card(T), we have 
card(T) -rk(MpTr) < 1. Since card(T) -rk(Mp^) > 0, for any T, either card(T) -rk(Mp^) = 
or card(T) — rk(Mp^) = 1. From this fact, it is easy to check the existence of a rank essential 
set T. For the uniqueness, we assume that there exist two subsets Ti, T2 C {!,••• , m} which 
are rank essential. Then, we have 

^WtiutJ < rk(Mp,J + rk(Mp,J - rk(Mp,^,,J 

= card(Ti) - 1 + card(T2) - 1 - card(Ti n T2) = card(Ti U T2) - 2, 

which means that Mp is not of full rank, a contradiction. 

Let T be a rank essential set. By Theorem 14.191 [PjJigT H Q{uj}jgT is of codimension 
one, which means that the sparse differential resultant of Pq, . . . ,Pn involves of coefficients 
of Fi {i G T) only □ 

Using this property, one can determine which polynomial is needed for computing the 
sparse differential resultant, which will eventually reduce the computation complexity. 

Example 4.22 Continue from Example \3.14\ {Po,Pi} is a rank essential sub-system since 
they involve yi only. 

A more interesting example is given below. 

Example 4.23 Let F he a Laurent differential polynomial system where 

IPo = uooyiy2 + -uoij/s 
IPi = uioyiy2 + -uiiysyg 

P2 = U20yiy2 + ^21^3 

Tro (o) I (o) I (o) 

IP3 = -"soyi + ^*3iy2 + ""322/3 

where o is a very large positive integer. It is easy to show that P is Laurent differentially 
essential and P = {Po,Pi,P2} is the rank-essential sub-system. Note that all yi,y2,y3 o'^e in 
P. P is rank essential because yiy2 can be treated as one variable. 

The second application is to prove the dimension conjecture for a class of generic differ- 
ential polynomials. The differential dimension conjecture proposed by Ritt [45, p. 178] claims 
that the dimension of any component of m differential polynomial equations in n > m vari- 
ables is no less than n — m. In [16], the dimension conjecture is proved for quasi-generic 
differential polynomials. The following theorem proves the conjecture for a larger class of 
differential polynomials. 

k 

Theorem 4.24 Let Pj = UiQ + ^ UikM,ik{i = l,...,m; m < n) be generic differential 

k=l 

polynomials in n differential indeterminates Y and Uj = {uio, . . . , uu.). Then [Pi, ... , P^,] C 
Q(ui, . . . ,Um){Y} is either the unit ideal or a prime differential ideal of dimension n — m. 

Proof: Use the notations introduced in the proof of Theorem 14.191 with Mjo = 1. Let 
Zo = [Pi,...,P„] C Q{ui,...,u^,Y} andXi = [Pi,...,P™] C Q(ui, . . . , u^){Y}. Since P, 
contains a non- vanishing degree zero term Uio, it is clear that [Pi, . . . ,Pm] : m = Zq. 
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From the proof of Theorem 14.191 is a prime differential ideal with = (r/i, . . . , ??n; Ci; 
Till, . . . , uii^ ; • • • ; Cm; Umi-, • • • , Umi^) as a generic point. Note that rk(Mp) < m and two cases 
win be considered. If rk(Mp) < m, by Theorem 2:„ = [Fi, . . . , P^] n Q{ui, . . . , u^} is 
of codimension m — rk(M]p) > 0, which means that Xi is the unit ideal in Q(ui, . . . , Um){Y}. 
If rk(Mp) = m, by the proof of Theorem gTll d.tr.degQ(u)(Ci, . . . ,Cm)/Q(u) = m and 
Tu = [0] follows. Since Xq = Xi n Q{ui, . . . , u^, Y} and Zq is prime, it is easy to see that Xi 
is also a differential prime ideal in Q(ui, . . . , Um){Y}. Moreover, we have 

n = d.tr.degQ(u)(r/i,...,77„,Ci,...,Cm)/Q(u) 
= d.tr.deg Q(u) (r/i Ci Cm)/Q(u, Ci , 

+d.tr.degQ(u,Ci,...,Cm)/Q(u) 
= d.tr.deg Q(u) (r/i Ci Cm)/Q(u, Ci , 

Hence, d.tr.deg Q(u, Ci, . . . ,Cm)(m> • • • , f?n)/Q(u, Ci, • • • ,Cm) = n-m. Without loss of gen- 
erality, suppose r]i, . . . ,7]n^m are differentially independent over Q{uXi, ■ ■ ■ Xm)- Since 
lo = ^1 n Q{ui, . . . ,Um,,Y}, {yi, . . . ,yn-m} is a parametric set of Xi. Thus, [Pi, . . . ,Pm] C 
Q(ui, . . . , Um){Y} is of dimension n — m. □ 

By Theorem 14.171 Theorem I4.19| and Corollary 14.181 properties 1) and 2) of Theorem 
11.11 are proved. 

5 Basic properties of sparse differential resultant 

In this section, we will prove some basic properties for the sparse differential resultant 
R(uo, . . . ,u„). 

5.1 Sparse differential resultant is differentially homogeneous 

Following Kolchin [32j, we now introduce the concept of differentially homogenous polyno- 
mials. 

Definition 5.1 A differential polynomial p € J~{yo, ■ ■ ■ ,yn} is called differentially homoge- 
nous of degree m if for a new differential indeterminate X, we have p{XyQ,Xyi . . . ,Xyn) = 

A"'p(yo,yi,---,yn). 

The differential analogue of Euler's theorem related to homogenous polynomials is valid. 



• • • ) Cm) 

...,Cm)+rn. 



Tfieorem 5.2 13^ f € J-'{yQ,yi, . . . ,yn} is differentially homogenous of degree m if and 
only if 



'k + r\ {k) df{yo,...,yn) ^ [ mf r = 
r ) y^ Qy(k+^) ~ 1 r / 

Sparse differential resultants have the following property. 



EE 

j=o km 
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Theorem 5.3 The sparse differential resultant is differentially homogenous in each Uj which 
is the coefficient set o/Pj. 

Proof: Suppose ord(R, Uj) = hi > 0. Follow the notations used in Theorem l3.91 By Corollary 
13.121 R(u; (q, . . . , Cn) = 0. Differentiating this identity w.r.t. nj^'* (j = 1, . . . , k) respectively, 
we have 

~§R~ I ~m~ ( Mij(n) s ~aR~ / , Mij{ri) ,is ~aR~ / r Mij('7) i//\ , BR I /?t,:'i r ^ij (')) 1 ffa, ) _ n /r|*~l 

5+ +^( - 11^) +^( - (i)[ifeMio + ■ ■ ■ + - Ci = ° 



+ --- + r^(-(;::)[t^i'°') = o e^-*) 



In the above equations, -^Wr (A: = 0, . . . , /ij; j = 0, . . . , /j) are obtained by replacing UiQ by 



Ci = 0, 1, . . . , n) in each respectively. 
Now, let us consider Y^-^q Y^k>o ( fc '^) 

— j^fP) ■ Of course, it needs only to consider 
r < hi. For each r < hi and each j £ {1, . . . , /j}, 

= ir*)xOu,, + (r + U)x{^+')u[^ + ■■■ + jh,.) x (^Q^g'"^^ 

S^'r'* V / *J^Mio(>?)J V r J^ijMioin)) ~^ 

I 3R / {hi\„, rMiill!l]ih,~r) _ (r+l\( h, \ , ! Mjjjr,) ^(h,-r -1) _ (hA (hA (lu-r) M^jjjfl^ 

9«<'''' V yrJ^-VlMioiv)^ \ r )\r + ll"-ij\-M^oiv)' "' \ r ) \hj ""ij Mio{v) J 

- ('A,, ~gR~ (''^+1^,,' I AA„,(^i-'-) 3R I <'r\ ~gR~ / _ ...MiiSsl\ 

I (r+l\ —dR~ f _ „. M.j(r,) Y QR f _ „ M.j (»?) V^''^''^ 

It follows that E^i C:)-..^ + E^i ('■r)-^.^ + • • • + E^L, a^^'-"^ + 

By Corollary Em G = T.k>oT!^=o )4f TTR:^ S sat(R). Since ord(G) < ord(R), 

G can be divisible by R. In the case r = 0, ^ ^ u-J^^ = ^n- ■ -R for some m G Z. While 

j=0/fc=0 '5">V 

in the case r > 0, if G 7^ 0, it can not be divisible by R. Thus, in this case, G must be 
identically zero. From the above, we conclude that 

'^ + ^\ (fc) 9R _ / r^O 



EE 



r 



*J I mi? r = 



j=0 k>0 ^ ' ^"-ij 

By Theorem 15.21 R(uo, . . . ,u„) is differentially homogenous in each Uj and the theorem is 
obtained. □ 

With Theorem 15. 3 1 property 1) of Theorem 11.21 is proved. 
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5.2 Order bound in terms of Jacobi number 



In this section, we will give an order bound for the sparse differential resultant in terms of 
the Jacobi number of the given system. 

Consider a generic Laurent differentially essential system {Po,...,Pn} defined in ([2|) 
with Uj = {uio,Uii, . . . , UiiJ being the coefficient vector of Pj (i = 0, . . . , n). Suppose R is the 
sparse differential resultant oi Fq, . . . ,¥n- Denote ord(R, Uj) to be the maximal order of R 
in Uik {k = 0,. . . that is, ord(R, Uj) = maxfcord(R, ttjfc). If Uj does not occur in R, as 
shown in Example 13.141 then set ord(R, Uj) = — oo. Firstly, we have the following result. 

Lemma 5.4 For each i, z/ord(R, Uj) = /ij > 0, then ord(R, njfc) = hi {k = 0, . . . 

Proof: Firstly, we claim that ord(R, Ujo) = h^. For if not, suppose ord(R, Ujfc) = hi > 
for some k ^ 0. Then by differentiating R(u; (^O) • • • ) Cn) = w.r.t. we have 

^(^) (u;Co, . . . ,Cn) = 0, where Ci are defined in dH). By Corollarv 13.121 we have -^Sh- G 

sat(R), a contradiction. Thus, ord(R, Ujo) = hi. For each A; 7^ 0, ord(R, Uj^) < hi. If 
ord(R, Uik) < hi, differentiate R(u; (q, . . . Xn) = w.r.t. u[l'\ we have -^^(u; Co, • • • , Cn) • 

(-Sm) = 0- So -^(u;Co, . . . ,Cn) = and G sat(R), a contradiction. Thus, for 

each k = 0, . . . Ji, ord(R, Uj/c) = hi. □ 

Let A = (fflij) be an n X n matrix where a^j is an integer or —00. A diagonal sum of 
A is any sum oij^ + 02^2 + • • • + ani„ with ii, . . . , i„ a permutation of 1, . . . , n. If ^4 is an 
m X n matrix with M = min{m,n}, then a diagonal sum of ^ is a diagonal sum of any 
M X M submatrix of A. The Jacobi number of a matrix A is the maximal diagonal sum of 
A, denoted by Ja,c{A). 

Let ord(Pf,yj) = eij {i = 0,...,n;j = l,...,n) and Eord(Pi) = ord(Pf) = a. We 
call the (n + 1) x n matrix A = (cij) the order matrix of Pq, . . . ,Pn- By A:-, we mean the 
submatrix of A obtained by deleting the {i + l)-th row from A. We use P to denote the 
set {P^,...,P^} and by P^, we mean the set P\{Pf}. We caU Ji = 3ac{A.) the Jacobi 
number of the system P^, also denoted by Jac(Pj). Before giving an order bound for sparse 
differential resultant in terms of the Jacobi numbers, we first give several lemmas. 

Given a vector ^ = {ko, ki, . . . ,kn) € ^>o^, we can obtain a prolongation of P: 

p[i?] _ I \(W>N\[k 



[J{¥ff< (10) 



Let tj = maxjeoj + ko, eij + ki, . . . , Cnj + kn). Then pl^l is contained in Q[u[^], Y[^^]], where 
ul^l = Uf^ouf'l and YI^] = ^J=iyf'^. 

Denote u(¥^~^^) to be the number of Y and their derivatives appearing effectively in 
p[-^] . In order to derive a differential relation among Uj (i = 0, . . . , n) from pl-^1 , a sufficient 
condition is _^ 

|p[^l[ > zy(p[^l) + l. (11) 
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Note that iy(p[^]) < |Y[^] | = E"=i(*i + 1) = Ej=i niax(eoj + A^o, + A;i, . . . , Cnj + K) + n. 
Thus, if |P[^]| > Y[^] + 1, or equivalently, 

n 

ko + ki-\ h /cn > ^ max(eoj + /cq, eij + ki,..., e„j + A:„) (12) 

is satisfied, then so is the inequahty (fTT]l . 

Lemma 5.5 Let ¥ be a Laurent differentially essential system and ^ = {kQ,ki, . . . ,kn) G 
Z"q^ a vector satisfying H^) . Then ord(R, Uj) < ki for each i = 0, . . . ,n. 

Proof: Denote ml-^l to be the set of ah monomials in variables Y^'^\ Let I = (P'-^l) : m'-^l 
be an ideal in the polynomial ring Q[Y[-^1, u^-^l]. Denote U = uf^^ U"^q u^^'K Assume 
= Ek=o '^^kNik = 0, . . . , n). Let Qu = -(ELi UikNik/Nio)^^) for i = 0, 1, . . . , n; / = 
0, 1, . . . , fcj. Denote C = {U, Cofco, • • • , Coo, • • • , Cnfc„, • • • , Cno)- It is easy to show that (YI^I, C) is 
a generic point of X. Indeed, on the one hand, each polynomial in X vanishes at 

(Y[^1,C)- On 

the other hand, if / is an arbitrary polynomial in Q[Y[ 1, u^-^l] such that /(Y[^],C) = 0, sub- 
stitute uS = ((Ff -Efc=i UikNik)/N^Qf^ into /, then we have 11^=0 ^fo / = mod (p[^]), 
where /i G Q[Y[^], [/]. Clearly, /i = and / G Z follows. 

Let Xi = X r\ Q[u["^l]. Then Xi is a prime ideal with C, as its generic point. Since 
Q(C) C Q(Y[^1,C/), Codim(Xi) = \U\ + T.l=o{ki + 1) - tr.degQ(C)/Q > \U\ + lp[^l| - 
tr.degQ(Y[^],[/)/Q = [p[^| - [Y[^]| > 1. Thus, Xi / (0). Suppose / is any nonzero 
polynomial in Xi. Clearly, ord(/, Uj) < k^. Since Xi C [P^, . . . ,P^] : m fl Q{uo, . . . , u„} = 
sat(R), / G sat(R). Note that R is a characteristic set of sat(R) w.r.t. any ranking by 
Lemma 12.31 Thus, ord(R, Uj) < ord(/, Uj) < k^. □ 



Lemma 5.6 Let F be a Laurent differentially essential system and Ji > for each i = 
0, . . . , n. Then ki = Ji {i = 0, . . . ,n) satisfy [W\) in the equality case. 

Proof: Let A = {eij) be the (n + 1) x n order matrix of P, where eij = ovd{T f ,yj). Without 
loss of generality, suppose Jq = en + 622 + • • • + enn- 

Firstly, we will show that for each k ^ 1, en + Ji > eui + Jk- Since Jk is the Jacobi 
number of P^ and k ^ 1, Jk has a summand of the form eip^. Consider the longest sequence 
of summands in Jk in the following form: 

To = eip^ + ep^p2 + ■ ■ ■ + (^Pm-lPm 

and suppose Jk = Tq + Ti. Since Jk is a diagonal sum, pi ^ pj for 1 < i < j. For otherwise, 
Jk contains ep._^p. and ep^_^p^ as summands, a contradiction. Also note that pi ^ Q for 
1 < i < m. Now we claim that pm is either equal to 1 or equal to k. Indeed, if p„i = 1 
or Pm = k, Tq cannot be any longer and these two cases may happen. But if pm / 1 and 
Pm 7^ k, then we can add another summand Cp^p^^-^ to Tq, which contradicts to the fact 
that To is the longest one. Now three cases are considered. 
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Case 1) If pi = 1, Jfc = eii + Ti and eki + Jk = en + e^i + Ti . Since eki + Ti is a diagonal 
sum of P|, eki + Ti < Ji. Thus, en + Ji > e/ci + Jfc. 

Case 2) If = 1 for m > 1, Tq = eip^ + Cp-^p^ H h ep^_-^i. Since Jo = en H h e^n, 

To < en + epipi + - • • + ep„_ip„_i. For otherwise, since / 0, Tb + X;fcg{2,...,n}\{pi,...,p„_i} efcfc 
is a diagonal sum of Pg which is greater than Jq. Then e^i + Jfc = e^i + Tq + Ti < e^i + 
en + Spipi + ■ ■ ■ + ep^_ipm-i + 7i < en + Ji, where the last inequality follows from the fact 
that eki + epipi H + ep^_-^p^_^ + Ti is a diagonal sum of P|. 

Case 3) If pm = k, Tq = eip^ + ep^p^ + • • • + ep^_-^k- Then, similar to case 2), we can 
show that efei + eip^ + ep^p^ H h ep^^^k < en + e^fc + ep^p^ H h ep^_-^p^_^. Thus, 

efci + Jfc = e/ci + eipi + Cp^p^ H + ep^_-^fc + Ti 

< efcfc + en + e^^p^ H h ep^_-^p^_^ + Ti 

< en + Ji- 

Similarly, we can prove that for each j, Cjj + Jj > Ckj + Jk with < k < n. So 

n 

max(eoj + Jq, ■ ■ ■ , e^j + J„) = en + Ji + e22 + ^2 H h e^n + Jn 

= Jo + Jl + . . . + Jn- 

□ 

Corollary 5.7 Lei P 6e a Laurent differentially essential system and Jj > /or each i = 
0, . . . ,n. Then ord(R, Uj) < Jj (z = 0, . . . , n). 

Proof: It is a direct consequence of Lemma 15.51 and Lemma 15.61 □ 

The above theorem shows that when all the Jacobi numbers are not less that 0, then 
Jacobi numbers are order bounds for the sparse differential resultant. In the following, we 
deal with the remaining case when some Jj = — oo. To this end, two more lemmas are 
needed. 

Lemma 5.8 |3, \34f Let A be an m x n matrix whose entries are 's and 1 's. Let Jac(^) = 
J < min{m, n}. Then A contains an a x b zero sub-matrix with a + b = m + n — J . 

Lemma 5.9 Let P &e a Laurent differentially essential system with the following (n+ 1) x n 
order matrix 

A = ( "^^^ {-oo)rxt \ 

~ V ^21 A22 ) ' 

where r + t > n + \. Then r + t = n + 1 and Jac(yl22) ^ 0. Moreover, when regarded as 
differential polynomials in yi, . . . ,yr-i, {Po) • • • )Pr~i} is 0, Laurent differentially essential 
system. 

Proof: From the structure of A, it follows that the symbolic support matrix of P has the 
following form: 

V B21 B22 J 
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Since P is Laurent differentially essential, by Corollary 14.181 rk(Mp) = n. As rk(Mp) < 
rk(5ii)+rk((52i B22)), n < {n-t) + {n+l-r) = 2n + l- (r + t). Thus, r + t < n + 1, and 
r + t = n + 1 follows. Since the above inequality becomes equality, Bn has full column rank. 
As a consequence, rk(Mp) = rk(i?ii) + rk(i?22)- Hence, B22 is a t x t nonsingular matrix. 
Regarding Pq, . . . ,Pr-i as differential polynomials in yi, . . . ,yr-i, then Bu is the symbolic 
support matrix of {Pq, . . . ,Pr-i} which is of full rank. Thus, {Pq, . . . ,Pr--i} is a Laurent 
differentially essential system. 

It remains to show that Jac(^22) > 0. Suppose the contrary, i.e. Jac(yl22) = —00. Let 
^22 be a t X t matrix obtained from A22 by replacing — cxo by and replacing all other 
elements in ^22 by I's. Then Jac(A22) < t, and by Lemma |5.8| A12 contains an a x 6 zero 
submatrix with a + b = 2t — Jac(A22) > t + 1. By interchanging rows and interchanging 
columns when necessary, suppose such a zero submatrix is in the upper-right corner of ^22- 
Then 

Cll (-Oo)axfe 
C2I C22 



A 



22 



where a + b > t + 1. Thus, 



B' 



22 



Dn Oaxb 

D21 D22 



which is singular for a + b > t + 1, a contradiction. Thus, Jac(A22) > 0. □ 
Now, we are ready to prove the main result of this section. 

Theorem 5.10 Let ¥ be a Laurent differentially essential system and R the sparse differ- 
ential resultant o/P. Then 



ord(R, Ui 



—00 if Ji = —00, 

hi < Ji if Ji > 0. 



Proof: Corollarv 15.71 proves the case when Jj > for each i. Now suppose there exists at 
least one i such that Ji = —00. Without loss of generality, we assume = — 00 and let 
An = (eij)o<j<n-i;i<j<n be the order matrix of Pf^. By Lemma 15.81 and similarly as the 
procedures in the proof of Lemma 15. 9^ we can assume that An is of the following form 

. _ f All (-oo)rxt 
" " V ^1 ^22 

where r + t > n + 1. Then the order matrix of P is equal to 

A = ( (-oo)rxt 
V A21 A22 

Since P is Laurent differentially essential, by Lemma 15.9^ r + t = n + 1 and Jac(^22) ^ 
0. Moreover, considered as differential polynomials in yi, . . . ,yr-i, P = {po, . . . ,Pr~i} is 
Laurent differentially essential and An is its order matrix. Let J, = Jac{{Aii)~-). By applying 
the above procedure when necessary, we can suppose that Jj > for each i = 0, . . . , r — 1. 
Since [P] n Q{uo, . . . , u„} = [P] n Q{uo, . . . , u^-i} = sat(R), R is also the sparse differential 
resultant of the system P and u^, . . . , u„ will not occur in R. By Corollary 15.71 ord(R, Uj) < 
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Jj. Since Ji = Jac(A22) + Ji > Jj for < i < r — 1, ord(R, Uj) < Jj for < i < r — 1 and 
ord(R, Uj) = —oo for i = r, . . . , n. □ 

Corollary 5.11 Let P he rank essential. Then Ji > for i = 0, . . . ,n and ord(R, Uj) < Jj. 

Proof: From the proof of Theorem 15.101 if Ji = —oo for some i, then P contains a proper 
differentially essential subsystem, which contradicts to Theorem 14.211 Therefore, Jj > for 
i = 0, . . . ,n. □ 

The following example shows that in spite of Jj > 0, ord(R, Uj) = — oo may happen. 

Example 5.12 Let P = {Pq, Pi, P2, P3} be a Laurent differential polynomial system where 

Po = uoo + uoiyiy[y2y2 
Pi = uio + unyiy'iy2y2 

IP2 = U20 + U2iyi + 1*222/2 

Then, the corresponding order matrix is 



A 



Lt is easy to show that P is Laurent differentially essential and {Po,Pi} is the rank- essential 
sub-system. HereH = uqqUh—uqiUiq. Clearly, ord(R, uq) = ord(R, ui) = and ord(R, U2) = 
ord(R, U3) = —00, but Jq = 2, Ji = 2, J3 = 3, J4 = —00. Also note that in this example, the 

1 2 

sub-matrix An in the proof of Theorem 15.101 corresponds to | 1 2 



We conclude this section by giving two improved order bounds based on the Jacobi bound 
given in Theorem 15.101 

For each j € {1,. . . ,n}, let Oj = m.m{k G N| Bis.t. deg(Pf^, yj'^^) > 0}. In other words, 

Oj is the smallest number such that y^~^^ occurs in {Pg^, . . . ,P^}. Let B = [cij — Oj) be 
an (n + 1) X n matrix. We call Ji = 3ac{B~-) the modified Jacobi number of the system P^. 
Denote 7 = Yl^=iQ-j- Clearly, J, = Ji — j. Then we have the following result. 

Theorem 5.13 Let ¥ be a Laurent differentially essential system and R the sparse differ- 
ential resultant o/P. Then 



( 1 


2 


—00 


\ 


1 


2 


—00 










—00 




V 1 


—00 





/ 



ord(R, Uj 



—00 if Ji = —00, 

hi < Ji-^ if Ji> 0. 
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(k) 

Proof: Follow the notations given above. Let Fj be obtained from Pj by replacing yj by 



{j = 1, . . . ,n;k > Oj) in Pj (i = 0, . . . , n) and denote 



^o,...,Pn}. Since 



Ma 



v 











it follows that rk(M]p) = rk(Mp) = n. Thus, I = [P] PI Q{uo, . . . , u„} is a prime differential 
ideal of codimension 1. We claim that Z = sat(R). Suppose Pj = moMio + Tj and Pj = 
UioMio + T;. Let Ci = -Ti/Mio and 0^ = -f./Mi^. Denote u = Uf^ouA{itio}- Then 
C = (u, ^Oi • • • ) Cn) is a generic point of sat(R) and 6 = (u,9q, . . . , On) is a generic point of 
1. For any differential polynomial G £ sat(R), G(C) = = (E '/'(^)-^<^(u))/(nr=i ^fo ) 
where </>(¥) are distinct differential monomials in Y. Then F^{u) = for each (p. Thus, 
G{9) = (E'^W-^</<(u))/(nr=i^fo) = and G e Z follows. So sat(R) C I. In the similar 
way, we can show that I C sat(R). Hence, R is the sparse differential resultant of P. Since 
Jac(Pj) = Jac(Pj) — 7, by Theorem 15.101 the theorem is proved. □ 



Remark 5.14 Let 1^ = (e — eo, e — ei, . . . , e — e„,) where e = Yl7=o Clearly, ipl-^l | = 
ne + n + 1 = |yH| + 1 > |Y[^1| + 1. Then by Lemma\5^ deg(R, Uj) < e - Cj < s - Sj. Here 
Si is the the order o/ Pj = 0, . . . , n). If Li = e — ei — 7(P) where 7(F) = J2j=ii2j + ^j) 
and ej = minjjej — ord(Pf', ?/j)|ord(Pf', yj) 7^ —00}. By 147}, (Lq, . . . ,L„) also consists of 
a solution to Then deg(R, Uj) < Lj. One can easily check that Ji < Li < e — Ci for 

each i, and the modified Jacobi bound is better than the other two bounds as shown by the 
following example. 

Example 5.15 Let A = (ejj)o<i<n,i<j<n be the order matrix of a system P.- 





—00 





\ 


5 





—00 







3 


5 




V 5 


2 


—00 


/ 



Then {Jo, Ji, J2, J3} = {12,12,7,10}, {Lq, L^, L2, L3} = {13,13,13,13}, {e - eo, e - ei, e - 
62, e — 63} = {15, 15, 15, 15}. This shows that the modified Jacobi bound could be strictly less 
than the other two bounds. 

Now, we assume that P is a Laurent differentially essential system which is not rank 
essential. Let R be the sparse differential resultant of P. We will give a better order bound 
for R. By Theorem 14.211 P contains a unique rank essential sub-system Fj. Without loss 
of generality, suppose T = {0, . . . , r} with r < n. Let be the order matrix of ¥j and for 
z = 0, . . . , r, let Aj'^ be the matrix obtained from Aj by deleting the (i + l)-th row. Note 
that Aj-- is an r X n matrix. Then we have the following result. 
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Theorem 5.16 With the above notations, we have 



ord(R, Uj) 



— oo 



r 

r + 1, . . . ,n. 



Proof: It suffices to show that ord(R, Uj) < Jac(j4^.) for i = 0, ...,r. Let L 



M,, 



1, . . . , r) such 



Z]j=i '^ijVj for i = r + 1, . . . , n. Since Pt is rank essential, there exist -j^^^ 
that their symbolic support matrix B is of full rank. Without loss of generality, we assume 
that the r-th principal submatrix of B is of full rank. Consider a new Laurent differential 
polynomial system P = U {L^+i, . . . , L„}. This system is also Laurent differentially 
essential since the symbolic support matrix of -j^^, • • • , , Vr+i, ■ ■ ■ ,yn is of full rank. 

And R is also the sparse differential resultant of P, for Pj is the rank-essential subsystem 
of P. The order vector of Lj is (0, . . . , 0) for i = r + 1, . . . , n. So Jac(Pj) = Jac{Aj:j^) for 
i = 0, . . . , r. By Theorem 15.101 ord(R, Uj) < Jac{{Aj'-) for i = 0, . . . , r. □ 



Example 5.17 Continue from Example \4.23 The corresponding order matrix is 



A 





V o 







o 



1 
1 

o ) 



Here R = uoiUio{{u2iUio)' U20U11 - 'U2iiiio(u20^tii)') - uqiUiquIquI^. Clearly, ord(R, uq) = 
0,ord(R, ui) = ord(R, U2) = 1, and ord(R, U3) = —00. But Jq = Ji = J2 = o + 1, J3 = 1, 
and ord(R, Uj) <C Ji for i = 0, 1, 2. // using Theorem \5.1(A then Af consists of the first three 
rows of A and Jacobi numbers for Aj are 1, 1, 1 respectively, which give much better bounds 
for the sparse differential resultant. 



With Theorem 15.101 property 2) of Theorem 11.21 is proved. 



5.3 Differential toric variety and sparse differential resultant 

In this section, we will introduce the concept of differential toric variety and establish its 
relation with the sparse differential resultant. 

We will deal with the special case when all the Ai coincide with each other, i.e., Aq = ■ ■ ■ = 
An = A. In this case, A is said to be Laurent differentially essential when Aq, . . . , An form 
a Laurent differentially essential system. Let A = {Mq = (yW)"o,Mi = (Y[°1)"i, . . . , Mi = 
(yl"!)'^'} be Laurent differentially essential where G Then by Definition 13.61 

I > n and there exist indices ki,...,kn € {1,...,/} such that '■^^J^^^^ , . . . , ''^ [ Ij-jao" ^^e 
differentially independent over Q. Let 

Pi = UioMo + UiiMi + ■■■ + uuMi (i = 0, . . . , n) 

be n + 1 generic Laurent differential polynomials w.r.t A. 
Consider the following map 
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defined by 

...,u) = {{&Y\i&'^r\- ■ ■ , i&^r) (13) 

where P(/) is the /-dimensional differential projective space over £ and ^ = (^i,...,^n) 
G (f^)". Note that ((^[°1)"°, (?[°1)"S . . . , (C'°^)"0 is never the zero vector since E £^ for all 
i. Thus is well defined on all of (£^^)", though the image of (/>^ need not be a differential 
projective variety of P(/). Now we give the definition of differential toric variety. 

Definition 5.18 The Kolchin projective differential closure of the image of (pA ^■^ defined to 
be the differential toric variety w.r.t. A, denoted by Xj^. That is, Xj^ = 4>a{{£^)^) ■ 

Then we have the following theorem. 

Theorem 5.19 X^ is an irreducible projective differential variety over Q of dimension n. 

Proof Denote Pf = Yl[=o'^ik^k (i = 0, ...,n) where G m. Clearly, Mfc/Mo = 
Nk/No {k = l,...,l). Let J = [Nozi - NiZo, ...,NoZi- NiZq] : m G Q{Y; zq, ^i, • • • , zi}, 
where m denotes the set of all differential monomials in Y. Let rj = {rji, . . . , r/„) be a generic 
point of [0] over Q and v a differential indeterminate over Q(r?). Let 6 = {v, ^^v, . . . , /yp^^j ^) • 
We claim that (ry; 0) is a generic point of which follows that ^7 is a prime differential ideal. 
Indeed, on the one hand, since each Nqz^ — N^zo {i = !,...,/) vanishes at {rj; 6) and rj an- 
nuls none of the elements of m, (ry; 9) is a common zero of J . On the other hand, for any 
/ G Q{Y; zo, 2^1, . . . , which vanishes at {rj;6), let /i be the differential remainder of / 
w.r.t. NoZi — NiZQ (i = 1, . . . ,1) under the elimination ranking zi y . . . y zi y zq y Y . Then 
/i G Q{Y;zo} satisfies that Ngf = /i mod [iVo^i - NiZo,NoZ2 - N2Z0, . . . , NqZi - Nizq]. 
Since f{r];0) = 0, . . . ,i]n,v) = 0, and /i = follows. Thus, f € and the claim is 

proved. 

Let J\ = J T\ Ql^o, zi, . . . , z;}. Then Ji is a prime differential ideal with a generic 
point B. Denote z = (zq, -Zi, • • • , z;)- Fo'^ ^-^y / ^ ^1 ^ z, since z^f G Ji, zq/ vanishes 
at B and f((f) = follows. So / G Ji, and it follows that J\ : "l = J\. And for any 
f ^ Ji <Z J and any differential indeterminate A over Q{r],v), let /(Az) = X] </'(A)/(^(z) 
where 0(A) are distinct differential monomials in A and /^(z) G Q{z}. Then f{X9) = 
= J2H^)Ui^)- So each f^{e) = and ^ G Ji follows. Thus, /(Az) G Q{A}Ji. By 
Definition 12. 2^ Ji is a differentially homogenous differential ideal. Then V = V(j7i) is an 
irreducible projective differential variety in P(/). Since is a generic point of V, dim(y) = 
d.tr.degQ(^^, . . . , ^^)/Q = n. If we can show X^, = V, then it follows that X^ is an 
irreducible projective differential variety of dimension n. 

For any point ^ G (^^)"', it is clear that (C; -^o(O) ^i(O) • • • ^^liO) is a differential zero 
of J and consequently {NojC), N^jQ, . . . , NijC)) eYjJi) = V . So 0(0 = {No{C), Ni{C), . . . , 
Ni{^)) G V. Thus 0^((^^)") C V and Xa = 0^((£'^)") C y follows. Conversely, since 
^i^l) = (1' W)^ • • • ' ^ generic point ofV,VQ X^. Thus, = X^. □ 

Now, suppose zq,...,zi are the homogenous coordinates of P(/). Let 

Lj = UioZQ + -Uiizi H h lijiZ/ (i = 0, . . . , n) 
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be generic differential hyperplanes in P(/)- Then, clearly, Pj = Lj o (p^. In the following, we 
will explore the close relation between Res^ and the differential Chow form of Xj[. Before 
doing so, we first recall the concept of projective differential Chow form (|36j). 

Let V be an irreducible projective differential variety of dimension d over Q with a 
generic point ^ = (,^0, Suppose ^0 0- Let Lj = - Yl[=o UikZk {i = 0,...,d) 

be d + 1 generic projective differential hyperplanes. Denote Ci = ~ Ylk=i "^ikCo^Ck {i = 
0, . . . , d) and Uj = (ujO) • • • > 'Uu)- Then we showed in [36j that the prime ideal I(Co) • • • > Cd) 
over Q(UiUj\{uio}) is of codimension one. That is, there exists an irreducible differential 
polynomial F G Q{uo,...,Urf} such that 1{Co,uqi, . . . ,uqi; . . .]Cd,Udi, ■ ■ ■ ,Udi) = sat(F). 
This F is defined to be the differential Chow form of V(X) or X. We list one of its properties 
which will be used in this section. 

Theorem 5.20 136^ Theorem 4- 7] Let -F(uo, ui, . . . , u^) be the differential Chow form of V 
with ord(-F) = h and Sp = . Suppose that Uj are differentially specialized over Q to 

sets Vi C £ and Pj are obtained by substituting Uj by Vj in¥i{i = 0, . . . ,d). // Pj = = 
0, . . . ,d) meet V, then sat(F) vanishes at (vq, . . . , v^^). Furthermore, if F(vo, . . . , v^) = 
and Sp{vQ, . . . , v^;) ^ 0, then the d+1 differential hyperplanes Pj = {i = 0, . . . ,d) meet V . 

With the above preparations, we now proceed to show that the sparse differential resul- 
tant is just the differential Chow form of X_^. 

Theorem 5.21 Let Res^ be the sparse differential resultant 0/ Pq, . . . , P^. Then Res^ is 
the differential Chow form of Xj[ with respect to the generic hyperplanes Lj (i = 0, . . . , n) . 

Proof: By the proof of Theorem 15.191 is an irreducible projective differential variety 
of dimension n with a generic point (1, ^j^)- Let d = - Ylk=i '^ik^^ = 

0, . . . ,n). Then sat(Chow(X^)) = ]I((Co,^toi, • • • ,""0;; • • • ;Cn,Uni, . . . ,Uni)). And by the def- 
inition of sparse differential resultant, sat(Res_4) = I((Co; ^oii ■ ■ ■ jUqi; . . . ; CmUni, ■ ■ ■ , Uni))- 
By Lemma 12. 3| Chow(X_4) and Res_4 can only differ at most by a nonzero element in Q. 
Thus, Res^ is just the differential Chow form of X^. □ 

Based on Theorem 15.201 we give another characterization of the vanishing of sparse 
differential resultants below, where the zeros are taken from 8 instead of 8^. 

Corollary 5.22 Let L, = ViQZQ + ViiZi + - ■ ■ + Viizi = 0{i = 0, . . . , n) be projective differential 
hyperplanes over 8 with Vj = (^,0, . . . ,Vii). Denote ord(Res^) = h and Sr, = ^^^r^- If X_a 

meets Lj = (i = 0, . . . , n), then Res^(vo, . . . , v„) = 0. And if Res^(vo, . . . , v„,) = and 
'S'r(vo, . . . , v„) ^ 0, then meets Lj = (i = 0, . . . , n). 

Proof: It follows directly from Theorems 15.211 and 15.201 □ 

Example 5.23 Continue from Example \'j. 1 7\ Following the proof of Theorem \5.1iA consider 
J = \y\Z\ — yiZo,yiZ2 — yfzo] : m. It is easy to show that Xj^ is the general component of 
Z1Z2 — {zqZ2 — ZqZ2), that is, Xj[ = V(sat(ziZ2 — {^0^2 — z'qZ2))). And Res_A is equal to the 
differential Chow form of X_4 . 

By Theorems 15.191 and 15. 2H property 4) of Theorem 11.21 is proved. 
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5.4 Poisson-type product formulas 

In this section, we prove formulas for sparse differential resultants, which are similar to the 
Poisson-type product formulas for multivariate resultants [30] , 

Denote ord(R, Uj) hy hi {i = 0, . . . , n), and suppose ho > 0. Let u = U"^QUj \ {uqo} and 
'^oo" ^^)- Consider R as an irreducible algebraic polynomial r(n£")) in 
Qo[nQQ''^]. In a suitable algebraic extension field of Qo, r{u^°^) = has to = deg(r, Uqq"^) = 
deg(R,M[)o''^) roots ji, . . . ,'yto- Thus 

to 

R(uo, ui, . . . , u„) = ^ J] - 7r) (14) 

r=l 

where A G Q(ui, . . . , u„)[u{('°'\nQg'''']. For each r such that I < t <to, let 

Qr = Qo(7r) = Q(U) ' • • - 45""'^ ' 7r) (15) 

be an algebraic extension field of Qo defined by r{u^Q^) = 0. We will define a derivation 
operator 6r on so that Qr becomes a 5,-- field. This can be done in a very natural way. 
For e € Q(u), define Jt-c = 6e = e'. Define (5*noo = ^^00 for « = 0, . . . , /iq — 1 and 



5^ 



Since R, regarded as an algebraic polynomial r in u^^\ is a minimal polynomial of 7,-, Sr = 
— does not vanish at iiQQ°^ = jr- Now, we define the derivatives of S^-Uoq for i > /iq by 



induction. Firstly, since r(7T-) = 0, 5T-(r(7T-)) = Sr,| (hg)_ (5t-(7t)+^| {ho)_ =0, where T = 

R' — SrUqq"^"'^^ We define d^^^^uoo to be (5r(7r) = (h > • Supposing the derivatives 

of 6t'^~^'' uqq with order less than j < i have been defined, we now define S!^°~^''uqq. Since 

R^') = Sr45°^'^ + ?i is linear in n^''''^ we define <5>+^noo to be , . , 

In this way, {Qr, ^r) is a differential field which can be considered as a finitely generated 
differential extension of Q(u). Recall that Q(u) is a finitely generated differential extension 
field of Q contained in S. By the definition of universal differential extension field, there 
exists a differential extension field Q <Z 8 Q(u) and a differential isomorphism ip-r over 
Q(u) from {Qti^t) to {Q-,5). Summing up the above results, we have 

Lemma 5.24 (Qt-,<^t) defined above is a finitely generated differential extension field of 
Q(u), which is differentially Q{u) -isomorphic to a subfield of £. 



the symbol p 



Let p be a differential polynomial in -F{uo, ui, . . . , u„} = J-{\i, uqq}. For convenience, by 
,^ . , we mean substituting u^qq^^'^ by (5^7,- (i > 0) in p. Similarly, by saying 

p vanishes at u^°^ = 7,-, we mean p ,^ y = 0. It is easy to prove the following lemma. 

Lemma 5.25 Let p be a differential polynomial in J-'{u, uqq}. Then p £ sat(R) if and only 
if p vanishes at Uqq"^ = 7,-. 
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When a differential polynomial p € Q(u){Y} vanishes at a point t] € Q!^, it is easy to 
see that p vanishes at ^riij) ^ convenience, by saying is in a differential variety V 

over Q(u), we mean (Prir]) G V- 

With these preparations, we now give the following theorem. 

Theorem 5.26 Let R(uo, ui, . . . , u„) be the sparse differential resultant o/Pq; • • • jl^n with 
ord(R, uq) = /lo > 0. Let deg(R, Uqq°^) = to- Then there exist in an extension field 
(Qri ^r) of (Q(u), 5) for r = 1, . . . , to O'^d k = 1, . . . Jq such that 

to lo 

K = All{uoo + ^uokCrk)^^°\ (16) 

T = l k = l 

where A is a polynomial in Q(ui, . . . , u„)[uQ^°'\itQQ°^]. Note that equation il6\) is formal and 
should be understood in the following precise meaning: (uoo + X]p=i ^OpCrp)^'^''^ = ^'^'^uqo + 

'5^(E;=l"ope.p)• 

Proof: We will follow the notations introduced in the proof of Lemma 15.241 Since R is 
irreducible, we have R^o = ^rWKT (h,-,) ^ 0- ^et = R^p/R^o (p = 1, • • • ,^o), where 



R 



00 "OO =Tt 

Note that R^p and ^rp are in J>. We will prove 



,{^o)_^ 
-'00 ~ '■r 



7r = -^T°(^OlCrl + 'U02Cr2 H h UoIo(,tIq[ 

''Op 



Differentiating the equality R(u; Co> Ci) • • • j Cn) = w.r.t. Mq^''^ we have 



9R , OR ^_MoM^^Q^ (17) 



^45°) ^oo(»?) 



where ,^ . are obtained by substituting Q to Ujo (z = 0, 1, . . . , n) in ?• ^ . Multiplying uop 

•S^op 5«0p 
to the above equation and for p from 1 to Iq, adding them together, we have 



Thus, g = uop ^(^^) + Boo ^(ftp) £ sat(R). Since g is of order not greater than R, it 

must be divisible by R. Since q and R have the same degree, there exists an a € Q such that 
q = aR. Setting Uq^"^ = in both sides of g = aR, we have X]p=i ''^OpRrp + ^^ooRrO = 0- 
Hence, as an algebraic equation, we have 

^00 + ^ UQpirp = (18) 
p=l 
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under the constraint u^qq^ = 7r- Equivalently, the above equation is vahd in {Tr,6r)- As a 
consequence, 7,- = — "Uop^rp)- Substituting them into equation (^^, the theorem 
is proved. □ 

Note that the quantities ^^-p are not expressions in terms of yi. In the following theorem, 
we will show that if {Ai{i = 0, ...,n)} satisfies certain conditions. Theorem 15.261 can be 
strengthened to make ^rp as productions of certain values of yi and its derivatives. Following 
the notations introduced before Lemma I5.3H we have 

Theorem 5.27 Assume that 1 ) any n of the Ai {i = 0, . . . ,n) form a differentially indepen- 
dent set and 2) for each j = 1, . . . , n, ej E Span^jajfc — Oio : k = 1, . . . , /j; i = 0, . . . , n}. 
Then there exist r]-rk {t = 1, . . . ,to] k = 1, . . . ,n) such that 

R = ^n(-oo + E-o.^^^-^J (19) 



r=l ^ k=l 

to 



r=l 



{ho) 



where rjr = (rjri, ■ ■ .,rjrn)- 



Moreover, rjr {t = 1, . . . , to) lies on Pi, . . . 



Proof: Follow the notations in this section and those introduced before Lemma 15.311 By 
condition 1), each hi > 0. Denote Pf = MjPj (i = 0, ...,n) where Mi are Laurent dif- 
ferential polynomials. Then (r/; Co; ^oii ■ ■ ■ , uqi^; . . . 'Xn, Uni, ■ ■ ■ , Uni„) is a generic point of 
[P^,...,P^] : m where r] = (771, . . . , r/„) is a generic point of [0] over Q(u) and Q = 
~ Sfc=i ^tfc M-o(r;j • proof of Lemma [5. 3 H there exist Sj and Tj which are products of 

nonnegative powers of -^^-j such that Sjyj — Tj G [Pq^, . . . ,P^] : m. That is, r]j = Tj/ Sj 

for j = l,...,n, where Tj^Sj are obtained by substituting (liQo, • • . , ttno) = (Coi-'-iCn) 
in Tj,Sj respectively. Since R is an irreducible polynomial, every does not van- 

ishes at 45°^ = 7t- Let r]rj = Jlj^o)^^^ and rjr = (ry^i, . . . , 77™). By ([IT]), = 

n So /, X ; n so r /TTr\ (k) 

n n {vf^)(^ok-aoohk = / A. So n n ' 

S be the differential polynomial set consisting of -^^-j and (Sj)™'"'"^ (^) for all i = 

'^'^ik ^ 

0, . . . , n; A; = 0, . . . , Zj; j = 1, . . . , n and m € N. By Corollary 13.121 there exists a finite set Si 
of5andaGNsuchthat/?=( Syifi ft [{T, / S.^'Y'^^''^''' - ^ / ^ 



(aok-aoohk _ -gp- / ^p- 
""ofc ' ""00 



sat(R). By Lemma 15.251 H vanishes at Uq^°^ = 7^. And by the proof of Theorem I5.29| 

S n sat(R) = 0. So ^rk = By Theorem S R = ^nt°=i(^ioo + E ^^0^?^^)^''"). 

Thus, (UnD follows. 

To prove the second part of this theorem, we need first to show that S^rjrj ^ for 
each A; > 0. Suppose the contrary, that is, there exists some k such that d^rjrj = 0. From 
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Vrj = ^1 (^o) > SrVri = i^Y^^l (hn) = 0. Thus, S^+^i^Y^^ € sat(R). It follows that 



ry- = i=) = 0, a contradiction to the fact that rjj is a differential indeterminate. 

■' '-'j ^ ^ ^ 

Follow the above procedure, we can show that = -^Sh- / -^jxH- where -^S^ 

ik iV) ik 



ik 



,^ > . From ([20]) . it is easy to show that Ylk=n "^ik-^wi ~ some b in Q. So, 

'00 =7r _____ 9" -ft' 

for each i / 0, Ylk=o ^ik-^^ = 0. It follows that for each i / 0, Pi(r/r) = I]fc=o '^ikMik{r]r) = 
^^(Eto^^fc^) =0- Sory.liesonPi,...,P„. □ 

Under the conditions of Theorem 15.271 we further have the following theorem. 



Theorem 5.28 The elements rjr {t = I, . . . ,to) defined in Theorem 5.27 are generic points 
of the prime ideal [Pf , . . . : m C Q(u){Y}, where u = Uf^^Uj. 

Proof: Follow the notations in this section. Let JTq = [P]^,...,P^] : m C Q{Y, u} and 
J = [Pf , . . . ,P^] : m C Q(u){Y}. Similar to the proof of Theorem E31 it is easy to show 
that Jo is a prime differential ideal. And by condition 1), H Q{u} = [0]. Thus, J = [Jq] 
is a prime differential ideal and J fl Q{Y, u} = JTq. Let ^ = (^i, . . . ,6.n) be a generic point 
of J. Then (C;u) is a generic point of Jq. Let /3 = - Efc=i ^ofeAfofc(0/^oo(0- Then 

/3, uoi, . . . , uqzo; u) is a generic point of I = [P^,P^, . . . ,P^] : m C Q{Y; uq, u}. Since 
SjVj -Tj = l,...,n), = 3^(/3,noi,...,uo«o;u). 

By Theorem 15.271 r]r is a common non-polynomial solution of = (i = 1, . . . , n), thus 
also a differential zero of J. Recall rjrj = ^1 ('io)_ • If / is any differential polynomial in 
Q(u){Y} such that /(r/^) = 0, then f(^, ^)|^(ho)_^ = 0. There exist aj G N such that 
9 = U,S"'f{^,...,^) G Q{Y;uo,u}. Then g\Z)_^ = 0. By LemmaE^Sl g G sat(R) 

while Sj ^ sat(R). As a consequence, g{l3, uqi, . . . , uqi^; u) = and Sj{l3, uqi, . . . , uqiq] u) 7^ 
0. It follows that /(^i, . . . , ^n) = and f £ J. Thus, rj^ is a generic point of J. □ 

With Theorems 15.261 IB. 271 and 15.281 property 4) of Theorem 11.21 is proved. 



5.5 Structures of non-polynomial solutions 

In this section, we will analyze the structures of the non-polynomial solutions. Firstly, we 
will give the following theorem which shows the relation between the original differential 
system and their sparse differential resultant. 

Let Ai {i = 0, . . . ,n) be a Laurent differentially essential system of monomial sets. Then 
by Theorem I4.2H Ai{i = 0, . . . , n) can be divided into two disjoint sets {Ai : i G T} and 
{Ai : i G {0, 1, . . . , n}\T}, where T C {0,1,..., n} is rank essential. In this section, we 
will assume that {0, 1, . . . , n} is rank essential, that is, any n of the Ai{i = 0, . . . , n) form 
a differentially independent set, which is equivalent to the fact that each Uj occurs in R 
effectively. 
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Theorem 5.29 Let Pj = X^^'^q '^ik^ik (i = 0, . . . , n) he a system of differential polynomials 
such that any n of the Pj form a differentially independent set. Let R(uo,...,u„) he the 



sparse differential resultant ofFi with ord(R, Uj) = hi. Denote Qik = -^^Mj^ ^^Mjo 



and S to he the set consisting of {i = 0, . . . ,n) and {y^^^)i<i<nk>'d- Then 

[Po, . . . ,P„] : m = [R, {Qik)o<i<n;i<k<k] ■■ 5°° 

in Q{Y, uo, . . . ,u„}. 

Proof: Let I = [Pq, . . . , Pn] : m- Following the notations in the proof of Theorem 13.91 2^ is a 
prime differential ideal with a generic point {rj; Co, uoii • • • i uqi^; . . . ; Cn, Uni, ■ ■ ■ , Uni„) where 
Ci = — J2k=i ^ifc M-pfa) • definition of sparse differential resultant, XnQ{uo, . . . , u„} = 

sat(R). Differentiating the equality R(u; (qXi, ■ ■ ■ , Cn) = w.r.t. ttj-^'\ we have 



+ :r77rT(-ir^) = (20) 



SnJJ') MM 



where are obtained by substituting C^i to Ui^ii = 0,1,..., n) in -^S^. Let Qik 

9Kk 9u\^ 



-^Mik - Trh^io. Clearly, Q,k G I. 

Since any n of the Pj form a differentially independent set, ord(R, Uj) > 0. Substitut- 
ing Mik by [Qik + Mio-^)/-^ in each P^, we have P^ = ^1=0^^^^^^ = UioMio + 

EtiMQ'^k + M,o^)/^. So -^P, = J2'^^^u,kQ^k + {Eto^^k^)M^o■ Since 

Qik € X, X^fc^o ^ifc '^(^.) G Thus, there exists some a G Q such that X]fc=o ""tfc ^(^ ) = '^I^- 

It follows that Pj € [R, ((3iA:)o<j<n;i<fc<«J ^ 5*°°- For any differential polynomial / € I, there 
exists a differential monomial M G m such that Mf G [Pq, . . . ,Pn] C ([R, {Qik)o<i<n;i<k<ii] '■ 
S°°). Thus, / G [R, (Q^fc)o<i<n;i<fc<d : andX C [K, {Q^k)o<i<n;i<k<k] follows. 
Conversely, for any differential polynomial g G [R, {Qik)o<i<n;i<k<k] ■ S°°, there exist some 
differential monomial M and some 6 G N such that M(f3- ^fi , G [R, Qjfc] C X. Since I 

is a prime differential ideal, 5 G Z. Hence, Z = [R, {Qik)o<i<n;i<k<k] '■ S°°. □ 

Theorem 15.291 shows that under the condition R = 0, the non-polynomial solutions of 
Pj = 0{i = 0, . . . , n) are the solutions of some differential polynomials of two terms. 

Corollary 5.30 Let Pq, . . . ,P„ he a Laurent differentially essential system of the form ^ 

and R(uo, . . . , u„) its sparse differential resultant. Suppose ord(R, uq) = ho > and denote 

Sn = ) ■ Suppose that when Uj (i = 0, . . . , n) are specialized to sets Vj which are elements 
a«(V _ 

in an extension field of T, Pj are specialized to¥i{i = 0, . . . ,n). Suppose 5'r(vo, . . . , v„) 7^ 0. 

//Pj = 0(i = 0, . . . , n) have a common non-polynomial differential solution then for each 

k, we have 

MokiO _ dn 



(vo, . . . ,v„)/5r(vo, . . . , v„). (21) 
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Proof: Denote ¥^ = MjPj (i = 0, . . . ,n) where Mj are Laurent differential monomials.. By 
the proof of Theorem l5.29l for each k = 1, ... Jo, the polynomial SuMoMok — ^^^q) MqAI^q € 



[Pq^, . . . ,P^] : m. Thus, if ^ is a common non-polynomial differential solution of Pj = 0, 
then 5r(vo, . . . , v„) •Mofc(0 - t^(vo, • • . , v„)Moo(0 = 0. Since 5r(vo, . . . , v„) / 0, (EU) 

follows. □ 



We conclude this section by giving a sufficient condition for a differentially essential 
system to have a unique non-polynomial solution. 

Follow the notations in section [3^21 that is, Ai = {Mio, Mn, . . . ,Mj;^} are finite sets of 

h 

Laurent differential monomials where Mik = (Y^'^'l)"''', and Pj = ^ UikM-ik {i = 0, . . . 

fc=o 

aik G Z'^^^'"'"^) is an exponent vector written in terms of the degrees of yi, . . . , y„, y'l, . . . , y^, 
. . . , yi^\ . . . , Vn''^ . Let o = maxj{si}. Of course, every vector in Z"'^**^^) can be embedded 
in 'L^'^°^^\ Let ej be the exponent vector for yi in whose i-th coordinate is 1 and 

other coordinates are equal to zero. 

Lemma 5.31 Assume that 1 ) any n of the Ai {i = 0, . . . ,n) form a differentially independent 
set and 2) for each j = 1, . . . ,n, ej € Spangjajfc — Ojo : k = 1, . . . ,li;i = 0, . . . , n}. Denote 
ord(R, Uj) = hi where R(uo,...,u„) = ReSytQ,...,^,^ . Let Pj be a specialization of Fi with 
coefficient vector Vj (z = 0, . . . , n). // R(vo, . . . , v„) = and , (vq, . . . , v„) ^ for each 

i and k, then Pj = (i = 0, . . . , n) have at most one common non-polynomial solution. 

Proof: By hypothesis 1), each hi > 0. Denote P^ = MjPj where Mi are Laurent differential 
monomials. Similar to procedures to derive equation ([20]) . we have = ~^^/~^^) 



where are obtained by substituting d to Ujo {i = 0,1, ... ,n) in J^,^ . By hypoth- 

esis 2), there exist tjik € Z such that '^i k^jikictik — ckjo) = for j = 1, . . . ,n. So 

n,A0'" = w- Thus, u.,{^f = 1, = n..{^/^)'-*. It follows 

that Sjyj — Tj G [Pq',...,P^] : m where Sj and Tj are products of nonnegative pow- 
ers of -^^-j obtained from the above identity. Since — ^^(vq, . . . , v„) 7^ for each i 

and k, T,(vo, . . . , v„) • 5j(vo, . . . , v^) / 0. Let yj = Tj(vq^. . . , v„)/S'j(vo, . . . , v„). If 
^ = {£,1, . . . ,£,n) is any common non-polynomial solution of P^ = (i = 0, then 
^3 = Tj{vo, ... ,Vri)/ Sj{\ro, ... ,\rn) = ijj. Thus ^ = and P^ = (z = 0, . . . , n) 

have at most one common non-polynomial solution. □ 



Theorem 5.32 Assume that 1 ) any n of the Ai {i = 0, . . . ,n) form a differentially indepen- 
dent set and 2) for each j = 1, . . . , n, ej E Span^jajfc — Oio : k = 1, . . . ,li;i = 0, . . . , n}. 
Denote ord(R, Uj) = hi where R(uo, . . . ,u„) = ReSy^(,^...^^„. Let Pj be a specialization o/P, 
with coefficient vector ^i{i = 0, . . . Then there exists a differential polynomial set S C 

Q{uo, . . . , u„} such that V(R)\ U V{S) ^ and whenever (vq, . . . , v„) G V(R)\ U ^(5), 
_ s&s ses 

Pj = (i = 0, . . . , n) have a unique common non-polynomial solution. 
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Proof: Follow the notations in the proof of Lemma 15.311 Denote ¥^ = MjPj {i = 0, . . . ,n) 
where Mi are Laurent differential monomials. Then by the proof of Lemma 15. 3H there 
exist Sj and Tj which are products of nonnegative powers of -^^-j such that Sjuj — Tj € 

[Fq, . . . , P^] : m. Thus, R, Siyi — Ti, . . . , SnVn — Tn is a characteristic set of [P^, . . . , P^] : m 
w.r.t. any elimination ranking -< Vi < ■ ■ ■ ~< Un- 

Let S be the differential polynomial set consisting of J^Jh-) and (S'j)™'"*"^ (^)^'"'* for all 

i = 0, . . . , n; = 0, . . . , /i; j = 1, . . . , n and m G N. Firstly, we show that V(R)/ IJ Y{S) ^ 

Ses 

0. Suppose the contrary, viz. V(R) C |J Y{S). In particular, there exists one S & S 

Ses 

such that S vanishes at the generic point C of sat(R). It is obvious that -^^j does not 

vanish at C- If (5'j)'"+^ (^) vanishes at C for some m, (5^)"+^ (^)^"^ G sat(R). Since 
Srn+iyM _ g [pAf^ . . . ,P^] : m, it follows that 5f +^ e [P^, . . . ,Pi^] : m, a 

contradiction. 

Suppose (vq, . . . , v„) G V(R)/ |J Y{S). Since -^^(vq, . . . , v„) ^ for each i and k, 

s&s 9u\^ 

Tj{-VQ, . . . , v„) • 5j(vo, . . . , v„) / 0. Let = -^{^o, • • • , v„) and denote y = (yi, . . . ,i/n)- 
And for each m G N, yj"^ = (^)(")(vo, . . . , v„) / 0. Thus, y G (f^)'^. Since R,5iyi - 

Ti, . . . , /Snl/n — T„ is a characteristic set of [P^, . . . , P^] : m, if • MjPj = 0, mod [R, Siyi — 
Ti, ■ ■ ■ , Sny-n ~ Tri\ where ii is a product of powers of Hence, Mi{y) ■ Pi(y) = 0, which 

follows that Pi(i/) = 0. Thus, y is a non-polynomial common solution of Pj. On the other 
hand, for every non-polynomial common solution of Pj, Sj{vo, . . . ,Vn)yj — Tj{vQ, . . . , v„) 
vanishes at it. Thus, it must be equal to the point y. As a consequence, we have proved that 
Pj = have a unique common non-polynomial solution. □ 

Theorem 15.321 can be rephrased as the following geometric form. 

Corollary 5.33 Let Zi{Ao, . . . , An) be a subset of S^^^^ x ••• x consisting of points 

(vo, . . . , v„) for which the corresponding Laurent differential polynomials Fi = (i = 0, . . . ,n) 
have a unique non-polynomial common solution and Zi{Ao, . . . ,An) the Kolchin closure of 
Zi{Ao, . . . , An)- Then under the condition of Theorem \5.3^ we have Zi{Ao, . . . , An) = 
Y{sat{ResAo,...,An))- 

Example 5.34 Continue from Example \3.14\ Ln this example, the sparse differential resul- 
tant R o/Po,Pi,P2 is free from the coefficients 0/P2. The system can be solved as follows: 
yi can be solved from Pq = Pi = and ¥2 = uiq -\- uiiy'2 is of order one in y2 which will 
lead to an infinite number of solutions. Thus the system can not have a unique solution This 
shows the importance of the first condition in Theorem \5.3SX 

Example 5.35 Continue from Example \3.15[ Ln this example, the characteristic set of 
[Po,Pi,P2] w.r.t. the elimination ranking -< y2 ^ Vi is R, uiiUooy2 ~ ^01^102/2) ""012/22/1 + 
uqq. Here Ao,Ai,A2 do not satisfy condition 2) and the system {Po,Pi,P2} does not have 
a unique solution under the condition R = 0. 
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5.6 Sparse differential resultant for differential polynomials with non- 
vanishing degree zero terms 

As pointed out in the preceding sections, for Laurent differential polynomials, non-polynomial 
zeros are considered. But, for certain differential polynomials, this condition seems to be 
too demanding. In this section, we restrict to consider the sparse differential resultant for 
differential polynomials with non-vanishing degree zero terms. To be more precise, consider 
n + 1 differential polynomials of the form 

h 

Pi = Uio + Y^ UikMik (i = 0, . . . , n) (22) 

k=l 

where Mit = (y['*'1)"''= is a monomial in {yi, . . . , y„, . . . , y^''\ . . . , whose exponent 

vector aik G lH^^^'^'^ with \aik\ > 1, and all the Uik are differentially independent over Q. 
Denote Bi = {1, Mn , . . . , Mu^ } (i = 0, . . . , n). The set of exponent vectors = {0, aik ■ k = 
1, . . . is called the support of Pj, where is the exponent vector for the constant term. 
Denote Uj = {mo, . . . , uuj (i = 0, . . . , n) and u = UjUj\{njo}. 

Definition 5.36 Let Pj(i = 0, . . . ,n) be a differential polynomial system of the form i22\) . 
{Po, . . . ,Pra} is called a differentially essential system if they form a Laurent differentially 
essential system when considered as Laurent differential polynomials. Ln this case, we also 
call Bq, . . . ,Bn ffl differentially essential system. 

All results for sparse differential resultants proved in the previous sections can be natu- 
rally rephrased in this case by just setting M^o in ([2]) to 1. The main difference is that we 
do not need to consider non-polynomial solutions and hence results in section 15.51 could be 
modified. 

First, we show that Theorem 13.91 can be modified as follows. 

Theorem 5.37 Let Pq, . . . ,P„ be differential polynomials as defined in \2S\) . Then [Pq, . . . , 
Pn] is a prime differential ideal in Q{Y, Uq, . . . ,u„}. And ([Po,Pi, . . . ,Pn]) nQ{uo, . . . ,u„} 
is of codimension 1 if and only if {Pj, i = 0, . . . , n} is a differentially essential system. 

Proof: Let r] = {rji, . . . , rjn) be a generic point of [0] over Q(u). Denote Ci = ~ Ylk=i '^ikMik{ri). 
It is easy to show that (ry; C) is a generic point of [Pq, Pi, ... , Pn] C Q{Y, uq, . . . , u„} where 
C = (Co,^ioi, • ■ ■ , uqi^{, Cn,Uni, ■ ■ ■,UniJ, and it follows that [Po,Pi, . . . ,Pn] is a prime 
differential ideal. So [Pq, Pi, . . . , P„] : m = [Pq, Pi, • • • , P„,]. By Theorem [3Jl the second part 
follows. □ 

Then, for a differentially essential system {Pq, . . . ,P„} of form (I22p . its sparse differential 
resultant R can be defined as 

[Po, . . . , P„] n Q{uo, . . . , u„} = sat(R). (23) 

This equation is different from equation ^ in that the differential ideal [Pq, . . . ,Pn] here 
is a differential ideal in Q{Y, uq, . . . ,u„} while the other one is generated in the Laurent 
differential polynomial ring. 



50 



Now we first introduce some notations similar to Section 15.51 Let Bq, . . . , Bn be a dif- 
ferentially essential system of monomial sets. For every specific differential polynomial set 
(Fo, ...,Fn) with Fi = Vio + Ylk=i ^ik^ik ^ we also represent it by (vq, . . . , v„) G 

£■'0+1 X • • • X where Vj = {vio, vn,..., vuj. Let 

Zo{Bo, ...,Bn) = {(vo, . . . , v„) G f'o+i X ■ ■ ■ X ^'"+1 : Fo = • • • = F„ = have 

a common solution in if"} (24) 



Let Z{Bo, Bn) be the Kolchin differential closure of Zo{Bo, . . . , S„) in £^o+i x---x £^"+\ 
Note that zeros from £ are considered, instead of £^ as in ([6]). 

The following result shows that the vanishing of sparse differential resultant gives a 
necessary condition for the existence of solutions, and as well as gives a sufficient condition 
in some sense. 

Theorem 5.38 Suppose Bi{i = 0,...,n) form a differentially essential system. Then we 
have Z{Bo, . . . = V(sat(ResBo,...,i3„)) . 

Proof: Since sat(Reseo,...^B„) C [Pq > • • • > I^n] C Q{Y; Uq, . . . , u„}, it follows directly that 
Zo{Bo, ...,Bn)'^ V(sat(Reseo,...,e„)). Consequently, Z{Bo, ■■■ ,Bn) ^ V(sat (Resg o,...,g„)) . 

For the other direction, follow the notations in the proof of Theorem 15.371 By The- 
orem (5371 [Po,...,Pn] is a prime differential ideal with a generic point (r/, C) where rj = 
{rji,..., r]n) is a generic point of [0] over Q(u) and C = (Co, ""oi, • • • , uqio'^ . . . ; Cn, . . . , UniJ. 
Let {Fq, . . . ,Fn) be a set of differential polynomials represented by C- Clearly, r/ is a solu- 
tion of Fi = 0. Thus, C G Zq{Bo, . . . , Bn) C Z{Bo, . . . , Bn)- Since C is a generic point 
of sat(Rese()^...^g,J, it follows that V(sat(ResBQ_,,,^e^)) C Z{Bq, . . . ,Bn)- As a consequence, 
V(sat(ResBo,'...,Bj) = ^(^o, • • • , Bn). ' ' □ 

In the following, we will analyze the properties of the solutions as we did in section 15.51 
The following lemma shows the relation between the original differential system and their 
sparse differential resultant, which is a direct consequence of Theorem 15.291 

Lemma 5.39 Let Pj = t*io + X]A;=i ^ik^hk (i = 0, . . . , n) he a system of differential polynomi- 
als satisfying that any n of the Pj form a differentially independent set. Let R(uo, . . . , u„) he 
the sparse differential resultant of¥i with ord(R, Uj) = hi. Denote Qik = — ^S^Mj^ ^Sh- 



and S to he the set consisting of — (i = 0, . . . , n). Then 

9Kq 



[Po, . . . , Pn] — [R, {Qik)o<i<n;l<k<li\ ■ 

in Q{Y, uq, . . . ,u„}. 

Proof: It is a direct consequence of Theorem 15.291 by setting Mjo = 1 and from the fact that 
[Po, . . . , P„] : m = [Po, . . . , P„] as differential ideals in Q{Y, uq, . . . , u„}. □ 

Lemma 15.391 shows that under the condition R = 0, the solutions of Pj = (i = 0, . . . , n) 
generally are the solutions of some differential polynomials of two terms. 
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Corollary 5.40 Let Po)---iIPn be a differentially essential system of the form ^22\) and 
R(uo, . . . , u„) their sparse differential resultant. Suppose ord(R, uq) = ho > and denote 
Sn = ) • Suppose that when Uj (i = 0, . . . , n) are specialized to sets Vj over Q which are 

elements in an extension field of J-', Pj are specialized to Pj (i = 0, . . . , n). //S'r(vo, . . . , v„) ^ 
0, in the case that Pj = 0(i = 0, . . . , n) have a common differential solution then for each 
k, we have 

MokiO = — (7-t(vo,---,v„)/S'r(vo,...,v„). (25) 
Proof: By Lemma [5. 39 1 for each k = 1, . . . Jq, the polynomial SuMok ^ [1^0) • • • 7 ^n]- 



Thus, if ^ is a common differential solution of Pj = 0, then S'r(vo, . . . , v^) • MokiO 



an 



^(vo, . . . , v„) = 0. Since 5r(vo, ...,Vn)^0,m follows. □ 



In the rest of this section, we will gives a sufficient condition for a differentially essential 
system to have a unique solution. 

Theorem 5.41 Assume that 1) any n of the Bi {i = 0, . . . ,n) form a differentially indepen- 
dent set and 2) for each j = 1, . . . ,n, ej S Span^jaj^ : /c = 1, . . . , /j; i = 0, . . . , n}. Denote 
ord(R, Uj) = hi where R(uo, . . . , u„) = ReSy^Q^...^^^ . Let Pj be a specialization of Pj over Q 



with coefficient vector Vj (z = 0, . . . , n). // R(vo, . . . , v„) = and vj; , (vq, . . . , v„) 7^ /( 
each i and k, then Pj = = 0, . . . , n) have a unique common solution. 
Proof By hypothesis 1), each hi > 0. By LemmalLMl Qik = -^Mik ^ G [Pq, • • 



or 



C Q{Y; uo, . . . , u„}. Since (r/; C) is a generic point of [Pq, . . . ,P„], Mik{7]) = -^/-^, 
where are obtained by substituting Q to Uio (i = 0, 1, . . . , n) in -^m-j- By hypothesis 2), 

du,' du,' 

ik rk 



there exist tjik € Z such that J2i k tjikOtik = for j = 1, . . . , n. So Hj k i^ik) = Uj- Thus, 
U^,k {M^k{v)Y'" = Vj = U^,k /Jhf" ' I* foUows that S,y,-T, G [Pq, . . . ,P„] where 



ik iO 

Sj and Tj are products of nonnegative powers of — obtained from the above identity. 

Thus, R, Siyi — Ti, . . . , Snyn — T„ is a characteristic set of [Pq, . . . , Pn] w.r.t. any elimination 
ranking Uik -< yi -< • • • -< Un- For each Pj, there exists a product Ai of nonnegative powers 
of ) such that AjPj G [R, ^lyi — Ti, . . . , 5„|/„ — T„]. Now specialize Uj to Vj over Q{i = 

0, ...,n). Since -^(vq, . . . , v„) / for each i and k, Tj(vo, . . . , v„) • 5j(vo, . . . , v„) / 0. 

Let yj = Tj{vo, . . . , v^)/5j(vo, . . . , v„) and y = (y^ . . . Clearly, ¥i{y) = 0. Thus, y is a 
common solution of Pj (i = 0, . . . , n). 

On the other hand, for any solution ^ of Pj, 5j(vo, . . . , v„)yj — Tj(vo, . . . , v„) (j = 
l,...,n) vanishes at it. So = y. As a consequence, we have proved that in this case, 
Pj = (i = 0, . . . , n) have a unique solution. □ 
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6 A single exponential algorithm to compute the sparse dif- 
ferential resultant 



In this section, we give an algorithm to compute the sparse differential resultant for a Laurent 
differentially essential system with single exponential complexity. The idea is to estimate 
the degree bounds for the resultant and then to use linear algebra to find the coefficients of 
the resultant. 

6.1 Degree of algebraic elimination ideal 

In this section, we will prove several properties about the degrees of elimination ideals in 
the algebraic case, which will be used to estimate the degree bound for sparse differential 
resultants. 

Let P be a polynomial in K\X\ where K is an algebraic field and X = {xi, . . . , x„,} a set 
of algebraic indeterminates. We use deg(P) to denote the total degree of P. Let Z be a prime 
algebraic ideal in K\X\ with dimension d. We use deg(X) to denote the degree of I, which 
is defined to be the number of solutions of the zero dimensional prime ideal (X, Li, . . . ,Lrf) 
in K{U)[X], where Lj = Ujo + Yl^=i '^ij^j = 1, . . . , d) are d generic hyperplanes [23] and 
U = {uij (z = 1, . . . , d, j = 0, . . . , n)}. That is, 



Clearly, deg(X) = deg(X, Li, . . . ,Lj) for i = 1, . . . , d. deg(X) is also equal to the maximal 
number of intersection points of V(X) with d hyperplanes under the condition that the 
number of these points is finite |25) . That is, 



We investigate the relation between the degree of an ideal and that of its elimination ideal 
by proving Theorem 16.21 

Lemma 6.1 Let I be a prime ideal of dimension zero in K\X\ and = Z n K[xi, . . . , x^] 
the elimination ideal of X with respect to xi, . . . ,Xk- Then deg(Xfc) < deg(X). 

Proof: Since both I and Ik are prime ideals of dimension zero, deg(Xfc) = \Y{Ik)\ and 
deg(X) = |V(X)|. To show deg(Zfc) < deg(X), it suffices to prove that every point of V(Xfc) can 
be extended to a point of V(X). Let (^i, . . . , ^fc) € Y{Ik)- For any point (771, ... , rjn) € Y(I), 
(771, ... , rjk) is a zero point of Ik- So we have -fC(^i, . . . , Ck) — ^ivi^ • • • > By [52l Propo- 
sition 9, Chapter 1, §3], there exist ^k+i, ■ ■ ■ such that . . . , = K{r]i, . . . ,rjn)- 
Thus, (Ci) ■ ■ ■ ) Cn) is a zero of Z, which completes the proof. □ 



Theorem 6.2 Let I be a prime ideal in and Ik — IC^I^[xi , . . . , Xk ] for any 1 < k < n. 
Then deg(Xfc) < deg(X). 



deg(X) = |V(X,Li,...,Ld)|. 



(26) 



deg(X) = max{|V(X) Ci Hi Ci ■ ■ ■ D Hd] : Hi areaffine hyperplanes 

such that |V(X) n i?i n • • • n Fd| < 00} 



(27) 
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Proof: Suppose dim(X) = d and dim(Xfc) = di. Two cases are considered: 

Case (a): di = d. Let Pj = Uio + uaxi + • • • + nj^x^ {i = 1, . . . ,d). Denote u = {uij : i = 
1, . . . ,d;j = 0, . . . ,k}. Then by [231 Theorem 1, p. 54], J' = (X^, Pi, . . . , P^) is a prime ideal 
of dimension zero in K{u)[xi, . . . , x^] and has the same degree as Ik- We claim that 

i) {I,¥i,...,¥a)nK{u)[xi,...,Xk] = J. 

ii) {I, Pi, ... , Pcf) is a 0-dimensional prime ideal over K{u). 

To prove i), it suffices to show that whenever / is in the left ideal, / belongs to J'. 
Without loss of generality, suppose / G -ftr[u][xi, . . . Then there exist hi,qi G JC[u][X] 
and gi ^ I such that / = Yli ^i9i + X^f=i Qi^i- Substituting Uio = —unxi — ■ ■ ■ — Ui^Xk into 
the above equality, we get / = Ylii f^i9i S I- Thus, / G X^. But f = f mod(Pi, . . . ,Pd), so 
f £ J', and i) follows. 

To prove ii), let (,^i, . . . , be a generic point of X. Denote Uq = {uio, . . . , Udo}- Then 
J^o = {I,Fi, . . . ,Prf) C i^(u\[/o)[X, C/q] is a prime ideal of dimension d with a generic point 
(6, • • • , Cn, Ci, • • • , Cd), where Ci = - I]j=i '^ijCj = 1, • • • , d)- Since di = d, there exist d 
elements in {^i, . . . algebraically independent over K. So by [22, p. 168-169], Cit ■ ■ Xd 
are algebraically independent over K{u\Uo). Thus, j7o H K{u\Uo)[Uo] = (0) and ii) follows. 

By Lemma EU deg(J) < deg(X,Pi, . . . ,Pd). So by ([271), deg(X) > |V(X,Pi, . . . ,Pd)| > 
deg(J) = deg(Xfc). 

Case (b): di < d. Let Lj = Ujo + "UjiXi + • • • + UinXn {i = 1, . . . , d — di). By [23^, Theorem 
1, p. 54], J' = (X,Li, . . . , hd-di) ^ ^(u)[X] is a prime ideal of dimension di and deg(^) = 
deg(X), where u = {uij : i = 1, . . . , d — di; j = 0, . . . , n}. Let Jk = J ^ K{vl)\x\, . . . , x^]. We 
claim that J^, = (Ik) in K{u)[xi, . . . ,Xk]- Of course, J'k 2 (Xfc). Since both J'k and (X^) are 
prime ideals and dim((Xfc)) = di, it suffices to prove that dim(j7fc) = di. 

Suppose (^1, . . . ,^„) is a generic point of X, then (^i, . . . ,^k) is that of X^. Let j7o = 
(X,Li, . . . ,Lrf_rfJ C K{u\Uo)[^,Uo], then (Ci,..., Cn, - E"=i ^ijO' • • • ' " Ej 
is a generic point of it, where Uq = {uio, • • • ,Ud-difi}- Since dim(Xfc) = di, without loss of 
generality, suppose ^i, . . . ,^di is a transcendence basis of K{^i, . . . ,S^jS)/K and ^i, . . . -.idii 
ik+ii ■ ■ ■ , S,k+{d-di) is that of K{^i, . . . , ^„)/i^. Then by [22l p. 168-169], it is easy to show that 
JonK{u\Uo)[xi, . . . ,Xd^,Uo] = (0), and JknK{u)[xi, . . . ,XdJ = follows. So dim(Jfc) = di 
and Jk = (Xfc). 

Since dim(^fc) = dim(J'), by case (a), we have deg(j7fc) < deg(j7') = deg(X). Due to the 
fact that deg(JA;) = deg((Xfc)) = deg(Xfc), deg(Xfc) < deg(X) follows. □ 

Li this article, we will use the following two results. 

Lemma 6.3 /5i, Corollary 2.28] Let Vi, . . . ,Vr C P" (r > 2) be pure dimensional projective 
varieties in . Then 

r 

l[degiV,)>Y,deg{C) 

i=l C 

where C runs through all irreducible components o/ Vi fl • • • H l^- 

Lemma 6.4 /55, Proposition 1, p. 151] Let Fi, . . . ,Fm G K\X] be polynomials generating an 
ideal I of dimension r. Suppose deg(Fi) > • • • > deg(Fm) and let D := ]Xi=i deg(Fj). Then 
deg(X) < D. 
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6.2 Degree bound for sparse differential resultant 

In this section, we give an upper bound for the degree of the sparse differential resultant, 
which will be crucial to our algorithm to compute the sparse resultant. 

Theorem 6.5 Let Po,...,Pn be a Laurent differentially essential system of form ^ with 
Eord(Pi) = Si and deg(Pf',Y) = m^. Suppose ¥f = X]jc=o '^"'^ ^'^ modified 

Jacobi number of {P^, . . . ,P^}\{Pf }. Denote m = maxj{mj}. Let R(uo, . . . , u„) be the 
sparse differential resultant ofFi{i = 0, . . . ,n). Suppose ord(R, Uj) = hi for each i. Then 
the following assertions hold: 

1) deg(R) < lYl=oimi + < (m + l)I^r=o(-^»+i), where m = maxi{mi}. 

2) R has a representation 

n n hi 

n iv(J^+^)'i^*5(^) • R = ^ ^ G,, (Pf ) (28) 

j=0 j=0 j=0 

where dj G Q[u[f . . . , ut'\Y^''^ and h = max{hi + ej such that deg(Gij(Pf < 
[m + l + Er=o(^* + l)deg(iV,o)]deg(R). 

Proof: 1) By the definition of sparse differential resultant, [Pq , . . . , P^] : mPl Q{uo, . . . , u„} = 
sat(R). Let rj = (771, ...,??„) be a generic point of [0]. Denote Q = -Y^k=i^ik ff^J^] (i = 

0, . . . Then (t/;Co,'Uoi, • • • j-uoio; • • • ;Cn,Uni, . . .,UniJ is a generic point of [P^, . . . ,Pn] : 
m. Clearly, Pq',...,P^ is a characteristic set of [Pq',...,P^] : m w.r.t. the elimination 
ranking n„,o y ■ ■ ■ >- uio y uqq ;^ u ^ Y. Taking the differential remainder of R w.r.t. this 
characteristic set, we have 

n hi 
i=0 k=0 

for Qi € N. Denote h = maxj{/ij + e^} and by m'^l we mean the set of all monomials in 
YN. Let J= ((P^)['^o],...,(PJJ)['^"]) iml'^l be an ideal in 7^ = Q[yW, uf''^, . . . , u[f"l]. Then 
R G jr. Furthermore, it is easy to show that is a prime ideal in TZ with a generic point 
(7?W; S, . . . , Ct^) and J D Qrf . . . , u^^^] = (R), where S = U,uf Let H.^ 

be the homogenous polynomial corresponding to (P^^) ^'^^ with xq the variable of homogeneity. 
Then = ((i^jfc)i<i<n;0<fc<hj : m is a prime ideal in Q[xo, YI'*], u|f°', . . . , uif'"'] with a 
generic point {v,vr]^'^^;vu,vCll^''\ . . . ,fCra'"b where m is the whole set of monomials in Y''^! 
and XQ. Then deg(J^°) = deg(J). 

Since V((i7ifc)i<j<„;0<fc</ii 
ducible component oiY{{Hik)i<i<n-o<k<h, 

). By LemmaESl deg(jO) < niLollfcLoK+l) = 
nr=o("^i + Thus, deg( J) < lYl=oi^^ + Since J n Q[nf'\. . . , u^^] = (R), 

by Theorem [O deg(R) < deg(J") < lYi=o{'mi + 1)^*+^ < (m + l)Er=o(-^i+i) follows. The 
last inequality holds because hi < Ji by Theorem 15.131 
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2) To obtain the degree bounds for this representation for R, we first substitute Uio by 
(P^ — Ylk=i '^ik^ik) l^io i^to R and then expand it. To be more precise, we take one mono- 
mial M(u; ttoo, • • • , Uno) in R(uo, . . . , u„) for an example. Denote M = -uC HiLo IlfcLoC^io^)'^'*' 
with I7I + Y1^=Q X^fcLo ^ifc ~ deg(R), where u''' represents a monomial in u and their deriva- 
tives with exponent vector 7. Substitute UiQ by (P^ — X]fc=i '^ik^ik) /^io into M, we have 

M(u; uoo, . . . , "«o) = n n ( (('f"^'^ - E "'fe^«0/^»o^ 



i=Ofc=0 ^ k=l 

When expanded, the denominator is of the form 11^=0 -^10"'°'^'^^^^'^''° ^'^'^ every term of the 
numerator has total degree I7I + Y17=oYlk=oi(^ + l)deg(iVio) + ("T-i + 1 - deg{Nio))]dik in 
Ug'*"', . . . ,u|f"' and Y'''] with h = max{/ij -|- Cj}. For every monomial M in R, Ylk=o(^ 

< {h, + l)deg(R). Thus, nr=o<"'''^'''^''^ ■ R = Er=oE'=oG.,(Pf)^'Vr where 
Gij,T G Q[uq^"', . . . , u[f"^, Y['*1] and T is free from Uio for i = 0, . . . , n. It is easy to see that 
T G [P^, . . . , P^] : m and T = follows. By the above substitution for every monomial in 
R, we can see that 



n hi 

< max(^,rf^^){|7| + X^Z^P + l)deg(iV*o) + K + 1 - deg(7Vio))]difc 

n n 

+ Y^i^i + l)deg(R)deg(A^,o) - ^^(^(^ + '^)d^k)deg{Nio)] 

i=0 i=0 k 

n 

< deg(R) + ^(mi - deg(7Vio))difc + deg(R) ^^^(/li + l)deg(iV, 

i=( 

n 

< (m + l)deg(R) + deg(R) ^ihi + l)deg(iV,o) 

i=0 

n 

+ 1 + + l)deg(iV,o)]deg(R). 



Ho) 

i=0 



\m 

i=0 



□ 



Example 6.6 Continue from Example \3.15\ In this example, Jq = 2, Ji = J2 = 1 and 
tuq = mi = m2 = 2. The expression of R shows that ho = ord(R, uq) = 1 < Jq, hi = 
ord(R, Ui) = < Ji [i = 1, 2) and deg(R) = 5 << S*^ = ULoi'^i + l)'''^^- 

For a differentially essential system of form (|22p , the second part of Theorem 16.51 can be 
improved as follows. 

Theorem 6.7 Let ¥q, . . . ,¥n be a differentially essential system of form \2S\) with m = 
maxi{deg(Pf , Y)} and Ji the modified Jacobi number o/{P^, . . . ,P.^}\{Pf }. Let R(uo, . . . , 
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u„) be the sparse differential resultant o/Pj (i = 0, . . . , n). Suppose ord(R, Uj) = hi for each 
i and h = max{/ij + Sj}. Then R has a representation 

n hi 

R(uo,...,u„) = ^X^G,,-pP 

n 

where dj G Q[uf . . . , u[f"^, yI'^I] such that deg(GijPp^) < (m+l)deg(R) < (m+l)»5'^'^'^^^^ 

Proof: Regarding Pj as Laurent differential polynomials, P^ = Pj and NiQ = 1. By setting 
NiQ = 1 in Theorem 16.51 these three assertions directly follow. □ 

The following result gives an effective differential Nullstellensatz under certain conditions. 

Corollary 6.8 Let fo, . . . , fn G J^{yi, • • • , Un} have no common solutions with deg(/i) < m. 
Let Jac({/o, . . . , fn}\{fi}) = Ji- If the sparse differential resultant of fo, . . . , fn is nonzero, 
then there exist Hij G T{yi, ...,yn} s.t. E"=o E/=o ^u/i^^ = 1 "'^'^ deg{Hijfj;^^) < (m + 
l)EILoW+i)+i. 

Proof: The hypothesis implies that P(/i) form a differentially essential system. Clearly, 
R(uo, . . . ,u„) has the property stated in Theorem 16.71 where Uj are coefficients of P(/i). 
The result follows directly from Theorem 16.71 bv specializing Uj to the coefficients of /j. □ 

With Theorem 16. 5| properties 6) and 7) of Theorem 1 1 . 2 1 are proved. 



6.3 A single exponential algorithm to compute sparse differential resultant 

If a polynomial R is the linear combination of some known polynomials Fi{i = 1, . . . , s), that 
is R = X]i=i -f^i^ij ^iid we know the upper bounds of the degrees of R and HiFi, then a 
general idea to estimate the computational complexity of R is to use linear algebra to find 
the coefficients of R. 

For sparse differential resultant, we already gave its degree bound and the degrees of the 
expressions in the linear combination in Theorem 16.51 

Now, we give the algorithm SDResultant to compute sparse differential resultants based 
on the linear algebra techniques. The algorithm works adaptively by searching for R with 
an order vector [ho, . . . ,hn) G N"+^ with hi < Ji by Theorem 16.51 Denote o = J2^=o^i- 
We start with o = 0. And for this o, choose one vector {ho, . . . ,hn) at a time. For this 
(/iQ, . . . , hn), we search for R from degree d = 1. If we cannot find an R with such a degree, 
then we repeat the procedure with degree d+1 until d > nr=o(™'« + 1)'*'"''^. In that case, we 
choose another (ho, . . . , hn) with Y17=o ^« ~ °- (^O' • • • ' with < Ji and 

X^ILo hi = a, H cannot be found, then we repeat the procedure with o + 1. In this way, we 
will find an R with the smallest order satisfying equation (|28p . which is the sparse resultant. 

Theorem 6.9 Let Pq, . . . ,P„ be a Laurent differentially essential system of form De- 
note P = {P^, . . . ,P^}, Ji = Jac(Pj), J = X;r=o Ji (^^d m = max^^odeg(Pj, Y). Algorithm 
SDResultant computes sparse differential resultant R o/Pq, . . . ,P„ with the following com- 
plexities: 
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Algorithm 2 — SDResultant(Po, . . . ,P„) 



Input: A generic Laurent differentially essential system ¥q, . . . ,¥n- 
Output: The sparse differential resultant R(uo, . . . , u„) of Pq, . . . , Pn• 
l. For i = 0, . . . , n, set = J2k=o u^kNik with deg(Aio) < deg{Nik). 

Set eij = ord(Pf rrii = deg(Pf ), rriio = deg(Aio), Uj = coeff(Pi) and \ui\ = /j + 1. 

Set A = (eij) and compute J, = Jac{A:-). 

2. Set R = 0, o = 0, m = maxjjmj}. 

3. While R = do 

3.1. For each vector (ho, . . . , hn) G N"^-'^ with ^"^q ^ ^^"^ ^* — '^i 

3.1.1. f7 = U^=ouf ^ = maxi{/i, + e^}, d = 1. 

3.1.2. While R = and d < IHLoC^-i + 1)'''^^ ^o 

3.1.2.1. Set Ro to be a homogenous GPol of degree d in U . 

3.1.2.2. Set Co = coeff(Ro, U). 

3.1.2.3. Set Hij{i = 0, . . . , n; j = 0, . . . , /ij) to be GPols of degree 
[m + l + Er=o(^» + ^)m^o]d - - 1 in Y^^\U . 

3.1.2.4. Set Cij = coeS{Hij,Y^^^ U [/). 

3.1.2.5. Set V to be the set of coefficients of nr=o ^l?'^^^'^R-o(uo, • • • , u„)- 
Er=oEi=o-^ij(IPf )^^^ as a polynomial in Y['^1,C/. 

3.1.2.6. Solve the linear equation P = in variables cq and c^. 

3.1.2.7. If Co has a nonzero solution, then substitute it into Rq to get R and go to 
Step 4, else R = 0. 

3.1.2.8. d:=d+l. 

3.2. o:=o+l. 

4. Return R. 

/*/ GPol stands for generic algebraic polynomial. 

/*/ coeff(P, V) returns the set of coefficients of P as an ordinary polynomial in variables V. 



1) In terms of the degree bound D o/R, the algorithm needs at most ("^^(■^+"'+^)) — 
<Q-arithmetic operations, where I = Y27=o(^i + -^^^^ ^/ 

2) The algorithm needs at most 0(^- —^i j Q- arithmetic opera- 
tions. 

Proof: The algorithm finds a differential polynomial P in Q{uo, . . . , u„} satisfying equation 
(|28p. which has the smallest order and the smallest degree in those with the same order. 
Existence for such a differential polynomial is guaranteed by Theorem 16.51 Such a P must 
be in sat(R) by equation Since each differential polynomial in sat(R) not equal to R 
either has greater order than R or has the same order but greater degree than R, P must 
be R. 

We will estimate the complexity of the algorithm below. Denote D to be the degree 
bound of R. By Theorem[631 D < [m + l)E?=o(-^>+i) = (m + l)-^+"+i, where J = Y.1=o -h- 
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In each loop of Step 3, the complexity of the algorithm is clearly dominated by Step 3.1.2, 
where we need to solve a system of linear equations "P = over Q in cq and Cjj. It is easy 
to show that |co| = C^+^Y^) and |c,,| = ^^ere L = ^^.{h, + m + 1) 

and di = [m + 1 + Yl^^^ihi + l)mio]d. Then P = is a linear equation system with 

N = ('^t^7VEr=o(^^ + l)r"T+i'^if ^'^) variables and M = (^^iit'iT^) 

To solve it, we need at most (max{M, N})^ arithmetic operations over Q, where w is the 

matrix multiplication exponent and the currently best known lo is 2.376. 

The iteration in Step 3.1.2 may go through 1 to nr=o("^* + 1)'''^^ < (m + l)Sr=o('^«+i)^ 
and the iteration in Step 3.1 at most will repeat nr=o('^« ~^ ^) t™6s. And by Theorem f 
Step 3 may loop from o = to X]r=o('^« ~^ '^^^ whole algorithm needs at most 

EE E Kx{M,iv})-- 

0=0 hi<.Ji d=l 



n 

< (j+i)(n('^^+i)) 



i=0 



{J 



EtoiJ'^ + ^)m^o]D + L + nih + l] 
L + n{h + 1) 



-I 2.376 



< iJ + n + 2)3-376( ELo(-/. + !) )»+! . ^ . r. ^ „ ^ . 2.376(L+„(.+1)) 

n + 1 

< {J + n + 2f-^7,iJ + n+ir+' ^ ^^^j ^^^^ 2^^^2.376((.+n)(.+ l)+n) 



arithmetic operations over Q. The above inequalities follow from the fact that h < J, 
L = Y17=oif^i + + 1) < IJ + I, L + n{h + 1) < {I + n)J + I + n = {I + n){J + I) and the 
assumption [m + 1 + J2i=oi'^i + ^)m]D > L + n{h + 1). 

Since / > 2(n + 1), the complexity bound is q^j- complexity 

assumes an 0(l)-complexity cost for all field operations over Q. Thus, the complexity 
follows. Now 1) is proved. To prove 2), we just need to replace D by the degree bound for 
R in Theorem 16.51 in the complexity bound in 1). □ 



Remark 6.10 As we indicated at the end of Section 3.3, if we first use Algorithm [1] to 
compute the rank-essential set T, then the algorithm can be improved by only considering the 
Laurent differential polynomials Pj (i G T) in the linear combination of the sparse resultant. 

Remark 6.11 Algorithm SDResultant can be improved by using a better search strategy. 
If d is not big enough, instead of checking d + 1, we can check 2d. Repeating this procedure, 
we may find a k such that 2^ < deg(R) < 2^+^. We then bisecting the interval [2^,2^+^] 
again to find the proper degree for R. This will lead to a better complexity, which is still 
single exponential. 



Remark 6.12 // the given system is algebraic, that is J = 0, then the complexity bound 
given in 1) of Theorem \6.9\ is essentially the same as that given in I49j[p. 288] since D ^ m 
and D ^ n. 
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For differential polynomials with non- vanishing degree terms given in (1220 . a better de- 
gree bound is given in Theorem 16. 7i Based on this bound, we can simplify the Algorithm 
SDResultant to compute the sparse differential resultant by removing the computation for 
and Nio in the first step where Nio is exactly equal to 1. 

Theorem 6.13 Algorithm SDResultant computes sparse differential resultants for a dif- 

ferentially essential system of form ^2^1 with at most Oy- j Q- 

arithmetic operations. 

Proof: Follow the proof process of Theorem 16.91 it can be shown that the complexity is one 
mentioned in the theorem. □ 

With Theorem 16. 9^ Theorem 11.41 is proved. 

6.4 Degree bound for differential resultant in terms of mixed volumes 

The degree bound given in Theorem l6.5l is essentially a Bezout type bound. In this section, a 
BKK style degree bound for differential resultant will be given, which is the sum of the mixed 
volumes of certain polytopes generated by the supports of certain differential polynomials 
and their derivatives. 

We first recall results about the degree of algebraic sparse resultant given by Sturmfels 
(|i50j). Let /C[X] = /C[a;i, . . . , be the polynomial ring defined over a field /C. For any vector 
a = (ai, . . . , a„) E Z*^, denote the Laurent monomial x\^X2^ ■ ■ ■ x^" by X°. Let Bq, . . . ,Bn C 
Z" be subsets which jointly span the affine lattice Z". Suppose = {0, . . . ,0) £ Bi for each 
i and \Bi\ = + 1 > 2. Let 

Fi(xi,...,x„) = c,o+ Yl c,,„X'^(i = 0,l,...,n) (29) 

aeBA{0} 

be generic sparse Laurent polynomials defined w.r.t Bi {i = 0,1, . . . ,n). Bi or {X" : a € Bi\ 
are called the support of Fj. Denote Cj = [cia)a&Bi ^ — Ui(cj\{cjo})- Let Qj be the 
convex hull of Bi in M", which is the smallest convex set containing Bi. Qi is also called the 
Newton poly tope of Fj, denoted by NP(Fj). In [50j, Sturmfels gave the definition of algebraic 
essential set and proved that a necessary and sufficient condition for the existence of sparse 
resultant is that there exists a unique subset {Bi}i^ i which is essential. Now, we restate the 
definition of essential sets in our words for the sake of later use. 

Definition 6.14 Follow the notations introduced above. 

• A collection of {Bi}i£ j, or {Fjjjgj of the form i29\) . is said to be algebraically inde- 
pendent i/ tr.degQ(c)(Fj — CiQ : i € J)/Q(c) = Otherwise, they are said to be 
algebraically dependent. 

• A collection of {Bi}i^ j is said to be essential if {Bi}i^j is algebraically dependent and 
for each proper subset J of I, {Bi}i^ j are algebraically independent. 

In the case that {Bq, . . . ,Bn} is essential, the degree of the sparse resultant can be 
described by mixed volumes. 
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Theorem 6.15 ([50j) Suppose that {Bq, . . . ,Bn} is essential. For each i E {0, 1, . . . ,n}, 
the degree of the sparse resultant in Cj is a positive integer, equal to the mixed volume 

JC{0,...,i-l,i+l,...,n} ieJ 

where vol(Q) means the n- dimensional volume of Q C M" and Q1 + Q2 means the Minkowski 
sum of Qi and Q2- 

The mixed volume of the Newton polytopes of a polynomial system is important in that 
it relates to the number of solutions of these polynomial equations contained in (C*)", which 
is the famous BKK bound. The following theorem explains it. 

Theorem 6.16 (Bernstein's Theorem) (f^) Given polynomials /i, •••,/« over C with 
finitely many common zeroes in (C)", let Qi be the Newton polytope of fi in M". Then the 
number of common zeroes of the fi in (C*)" is bounded by the mixed volume A4{Qi, . . . , Qn)- 
Moreover, for generic choices of the coefficients in the fi, the number of common solutions 
in (C*)" is exactly M{Qi, . . . , Qn)- 

It is well known that for a given polynomial system over C, the Bezout bound gives a 
bound for the number of isolated solutions in (C)". Comparing the BKK bound with Bezout 
bound, we have the following lemma. 

Lemma 6.17 Follow the notations in Theorem \6.16l Then A4{Qi, . . . , Qn) < Yli=i deg(/i). 

Proof: Suppose fi{i = l,...,m) are a system with generic coefficients. Then by Theo- 
rem 16.161 the number of common zeroes of the fi in (C*)" is equal to the mixed volume 
A4{Qi, ■ ■ ■ , Qn)- And by Lemma 16.41 the number of common zeroes of the fi in (C)" is 
bounded by Uti deg(/i). Thus, M{Qi,. . . , Q„) < EULi deg(/^) follows. □ 

In the rest of this section, the degree of sparse resultant will be used to give a degree bound 
for differential resultant in terms of mixed volumes. A system of n + 1 generic differential 
polynomials with degrees mo, . . . , m„ and orders sq, . . . ,Sn respectively of the form 

^^ = Uio+ Yl ni„(YW)"(i = 0,...,n), (30) 
1 < laj < rui 

of course forms a differentially essential system and their sparse differential resultant is 
exactly equal to their differential resultant defined in |16j . So Theorem 16.71 also gives a 
degree bound for differential resultant. But when we use Theorem 16. 71 to estimate the degree 
of R, not only Bezout bound is used, but also the degrees of Pj in both Y and Uj are 
considered. 

The following theorem gives a better upper bound for degrees of differential resultants, 
the proof of which is not valid for sparse differential resultants. Precisely, in the following 
result, when estimate the degree of R, the BKK bound is used rather than the Bezout bound 
as did in Theorem 16.71 
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Theorem 6.18 Let Pj (i = 0, . . . , n) he generic differential polynomials inY = {yi, . . . , 
with order Si, degree rui, and coefficients Uj respectively. Let R(uo, . . . , u„) he the differential 
resultant of Fq, ... ,Fn- Denote s = Ylii=QSi. Then for each i G {0,1, ... ,n} , 

S-Si 

deg(R, Uj) < ^ M{{Qji)jj^ifl<i<s-sj, Qio, ■ ■ ■ , Qi,k-i, Qi,k+i, ■ ■ ■ , Qi,s-s,) (31) 

k=0 



dIs-sq] 
' 



where Qji is the Newton polytope of as a polynomial in y^^^ , ■ ■ ■ , yn^ • 

Proof: By [16l Theorem 6.8], ord(R, Uj) = s — Si{i = 0,...,n) and (R) = 
IPn n Qlujf"*"', . . . , u^n~'"^]- Regard pf ^ {i = 0, ... ,n,k = 0, ... ,s - Si) as polynomials 
in the n{s + 1) variables = {y^, . . . ^ y^^ y[, . . . , y'^, . . . , y[^\ . . . , yn^}, and we denote its 
support by Bik. Let Fj^ be the generic sparse polynomial with support Bik. Denote Vj^ to be 
the set of coefficients of Fj^ and in particular, suppose Viko is the coefficient of the monomial 
1 in Fjfc. Now we claim that 

CI) B = {Bik : < i < n; < k < s — Si} is an essential set. 

C2) B = {Bik :0<i<n;0<k<s-Si} jointly span the affine lattice Z^^'^+i). 

Note that \B\ = n{s + 1) + 1. To prove CI), it suffices to show that any n(s + 1) of 
distinct Fj^ are algebraically independent. Without loss of generality, we prove that for a 
fixed k & {0, . . . , s — So}, 

Sk = {(Fj7)l<j<n;0<i<s-Sj)Foo, • • • ,Fo,fc-l,Fo,fc+li • • • 1^0, s-so} 

is an algebraically independent set. Remember that {yi, . . . , yn,y'i, ■ ■ ■ , y'm ■ ■ ■ , y[^''~^'^\ ■ ■ ■ , 
^(s,+0| -g ^ subset of the support of Fj/. Now we choose a monomial from each F^; and 

denote it by m{¥ji). For each j € {1, . . . , n} and I € {0, . . . , s — sj}, let m{¥ji) = yj'^^^'^ 
which belongs to the support of ¥ji. For the fixed k, there exists a r € {0, 1, . . . , n — 1} such 
that either Yll=i ^ ^ ^ Yli=i some r G {0, 1, . . . , n — 2} or ^J^i Si < k < J2l=i ^« 

for T = n — 1. Here when r = 0, it means < A; < si — 1. Then for / ^ k, let 

< / < si - 1 
y^2~''^ Sl<l<Si + S2-l 



01) 



. yn 22i=l Si + l<l< }2i=l Si = S-So 

So m{Sk) is equal to {y^^'^ : 1 < J < n}, which are algebraically independent over 
Thus, the n(s + 1) members of Sk are algebraically independent over Q. For if not. 
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¥ji — Vjio are algebraically dependent over Q(v) where v = ^f=oYlk=o '^ik\{'>^iko}- Now 
specialize the coefficient of m(¥ji) in F^; to 1, and all the other coefficients of ¥ji — vjio to 

0, by the algebraic version of Lemma [XTl {m(¥ji) : Fji £ Sk} are algebraically dependent, 
which is a contradiction. Thus, claim CI) is proved. Claim C2) follows from the fact that 
{1, yj'^' : 1 < j < n} is contained in the support of Fo,s_so. 

Prom the claims CI) and C2), the sparse resultant of {¥ik)o<i<n;0<k<s-Si exists and 
we denot e it by G. Then (G) = {{¥ik)o<i<n;0<k<s-s,) f] Q[{'^ik)o<i<n;0<k<s-s,], and by 
TheoremEjg deg(G,Vjfc) = M{{Qji)j^ifl<i<s-sp Qio, ■ ■ ■ , Qi,k-i, Qi,k+i, ■ ■ ■ , Qi,s~s,) ■ 

Now suppose is a generic point of the zero ideal (0) in Q(v)[Y[''1]. Let Qk = —¥ikiO~^''^iko 
and Cik — ~'^i^\0~^'^io^ {i = 0, . . . ,n; k = 0, . . . , s — Si). Clearly, Q and Ci are free of Viko and 
u['^^ respectively. It is easy to see that v. Coo, • • • , Co,s-so, • • • , CnO, • • • , Cn,s-sJ is a generic 
point of the algebraic prime ideal ((Fj^) 

0<i<n;0<k<s-Si) C Q[YW , (vjfc)o<i<n;0<fc<s-sJ 5 while 

(C; ^7=oi^i\{uio})^'~''h Coo> • • • , Co,s-.o' • • • ' C„0' • • • , C„,s-s„) is a generic point of the alge- 
braic prime ideal ((Pf ^)o<j<n;0<fc<s-sJ C Q[YH, u[f"''°', . . . , uL''"''"']. If we regard G as a 
polynomial in Viko over Q(v), then G is the vanishing polynomial of (Coo, • • • , Co,s-so, • • • , CnO, 
. . . ,Cn,s-s„) over Q(v). Now specialize the coefficients Vj^ of Fj^ to the corresponding coef- 
ficients of Pj-'^^ Then Ci are specialized to Ci- In particular, Viko are specialized to m-q^ which 
are algebraically independent over the field Q(C,U^^ouj'' ''\u% ^''). We claim that there 

- J 1-1 TT/, ,-n \s — Si]\ \s — Sj] (s — sn) (s — Sn)\ ^ 

exists a nonzero polynomial H[[J!f^QUi \uIq , uqq, . . . , Mqq , . . . , Uno, ■ ■ ■ , w„o ) ^ 
Qfujf'''^, . . . , uL'~'"'] such that 

C3) H{U^^qu[ '\u\q Coo, • • • , Co,s-so, ■ ■ ■ , CnO> • • • , Cn,s-Sn) ~ ^ 

C4) degiH,uf-''^) < deg(G,U^ro Vife). 

We obtain H by specializing v one by one in G. For each t; G v, denote u to be its 

(k) 

corresponding coefficient in P^- . Now we first specialize u to u and suppose Cifc is spe- 
cjalized to Cifc correspondingly. Clearly, G{v\{v}, n; Coo, • • • , Co,s-so, CnO, • • • , Cn,s-sJ = 0. If 
G = G{v\{v},u;vooo,voio, ■ ■ ■ ,vo,s-sofi, ■ ■ ■ ,VnOO,Vnio . . . ,v„^s-s„fi) 0, denote G by Hi. 
Otherwise, there exists some a G N such that G = {v — u)'''G\ with G\\i,=u 

/ 0. But 

G{-v\{v},U] Coo, • • • , Co,s-so, CnOj • • • , Cn,s-Sn) = = (^^ " uYGi{-v\{v}, U; Coo, • • • , Co,s-so, CnO, 

■ ■ ■ ,Cn,s-s„), SO Gi(v\{7;},u;Coo,---,Co,s-so,CnO,---,Cn,s-s„) = 0. Denote Gi|„=„ by Hi. 
Clearly, deg(i?i,u[* '''' IJUfcVifc) < deg(G, UfcVjfe) for each i. Continuing this process for |v| 
times till each f € v is specialized to its corresponding element u, we will obtain a nonzero 
polynomial i7|v|(U^=o("A{^^jo})'''"'''^;wooo,i'oio, • • • , ^'o.s-scO, • • • , fnoo, f^nio, • • ■ ,Vn,s-sr^fi) sat- 
isfying H\^\{\J'^^Q{Mi\{uio})^'-''^\CQQ, . . . , Co,s-.o' CnO, • • • , Cn,s-s„) = and moreover, for each 

1, deg(ff|v|, U Ufe{wifco}) < deg(G, UfcVjfc). Since are algebraically independent 
over the field Q(C, ur=o(ui\{n,o})[^-^'l), H = ii'|,|(Uto(uA^*o})[^-^'Uoo, • • • , 4r'°^ • • • > 
Uno, ■ ■ ■ , u"nO ^ '^[^0 ■ ■ ■ : ^"') is a polynomial satisfying C3) and C4). 

From C3), H G (p[f-^°], . . . , Plf"^"'). Since {¥^^-''\. . . ,p|f-^"])nQ[u[,^-^°], . . . , ulf"^"'] = 
(R) and R is irreducible, R divides H. It follows that deg(R, uj** '^'') < deg(i7, *'^) 

< deg(G,UfcVifc) = deg(G,Vifc) = -A^((Qii)jVi,o<z<s-s,, Qio, • • • , Qi,fc-i, Qi,fc+i, • • • , 
fc=0 fc=o 
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As a corollary, we give another degree bound for differential resultant by using Bezout 
bound, which is better than the bound given in Theorem 16.71 in that only the degrees of Pj 
in Y are considered in the bound. 



Corollary 6.19 Let Pj (i = 0, . . . , n) be generic differential polynomials inY = {yi, . . . , y„} 
with order Si, degree mj and coefficients Uj respectively. Let R(uo, . . . , u„) be the differential 
resultant o/ Pq, . . . , P„. Denote s = Y17=o ^i- Then for each i G {0, 1, . . . , n}, deg(R, Uj) < 

s-s,+l T-rn ^^-^J+^ 
m, 1 li=0 

Proof: Follow the notations in the proof of Theorem l6.181 Since {Bik : < i < n;0 < k < s — 
Si} is an essential set, for any fixed A; € {0, ... , s — Si}, the polynomials in Sk together generate 
an ideal of dimension zero in YW. By lemma[6J/71 M{{Qji)j^i^o<i<s-sp Qio, ■■■ , Qi,k-i, Qi,k+i, 
Qi,s-s,) < Wj=Q^'r'^^- Hence, by TheoremEIEl 

deg(R, Uj) < ^ A^((Qj7)jVi,o<«<s-s,, Qio, • • • , Qi,fc-i, Qj,fc+i, • • • , Qi,s~sJ 

A:=0 

s—Si n I 1 " 

^ sr-^ i TT + l S — Sj + i T-r s-s,+l 

< > — \ \m- ' = 

'-^ rrii ■' ruj -l-l J 

A:=0 j=0 j=0 

□ 



Example 6.20 Consider two generic differential polynomials of order one and degree two 
in one indeterminate y: 

Po = uoo + uoiy + uo2y' + uosy"^ + uoivy' + uo5{y')^ , 
Pi = uio + uny + uuy' + uisy'^ + uuyy' + ui5{y'f . 

Then the degree bound given by Theorem \6.5\ is deg(R) < (2 + 1)^ = 81. The degree bound 
given by Corollary \ 6. 19\ is deg(R, uq) < 2^ = 16 and hence deg(R) < 32. The degree bound 
deg(R,uo) given by TheoremWdE is M(Qio, Qii, Qoo) + M(Qio, Qii, Qoi) = 4 + 6 = 10 
and consequently deg(R) < 20, where Qoi = Qio = conv{(0, 0, 0), (2, 0, 0), (0, 2, 0)}, Qoi = 
Qii = conv{(0,0,0),(2,0,0),(0,2,0), (0, 0, 1), (1, 0, 1), (0,1,1)}, and conv(-) means taking 
the convex hull in R^. 

We will end this section by giving Algorithm DResultant to compute differential resul- 
tant based on the degree bound given in Theorem 16.181 

Theorem 6.21 Let Pq, . . . ,Pn be a generic differential polynomial system of the form (30\). 
Denote s = ^"^q ord(Pj, Y) and m = max"^Qdeg(Pj, Y). Algorithm DResultant computes 
the differential resultant R o/Pq, . . . ,Pn with the following complexities: 

1) In terms o/deg(R), the algorithm needs at most 0[{ns + n) 2-376 [^deg(R)]^('"("+i))) 
Q-arithmetic operations where I = max"^Q (™'^"^^'-|t^^) is the size of system Pj. 
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Algorithm 3 — DResultant(Po, . . . ,P„) 

Input: A generic differential polynomial system Pq) • • • i^^jt- 
Output: The differential resultant R(uo, . . . , u„) of Pq, . . . , Pn- 

1. For i = 0, ... ,n, set Si = ord(Pi), rrij = deg(Pj, Y) and Uj = coeff(Pj). 

2. Set R = 0, s = J27=o Si, m = maxilmj, d = n + l,U = U^^ouj'"''l 

3. While R = do 

3.1. Set Ro to be a homogenous GPol of degree d in U. 

3.2. Set Co = coeff(Ro,?7). 

3.3. Set Gik{i = 0, . . . ,n; k = 0, . . . , s — Si) to be GPols of degree (m + l)d — — 1 
in YH,[/. 

3.4. Set Cik = coefr(Ga., yW U U). 

3.5. Set V to be the set of coefficients of Ro(uo, . . . , u„) — Yli=o Ylk=o ^ik^i''' 
as a polynomial in YW, U. 

3.6. Solve the linear equation P = in variables cq and Cj^. 

3.7. If Co has a nonzero solution, then substitute it into Rq to get R and go to 
Step 4, else R = 0. 

3.8. d:=d+l. 

4. Return R. 

/*/ GPol stands for generic algebraic polynomial. 

/*/ coeff(P, U) returns the set of coefficients of P as an ordinary polynomial in variables U. 



2) The algorithm needs at most 0(^{ns + n)'^'^'^^[mD]'~^^''"'^^~^^^^') Q-arithmetic operations, 
where D is the degree bound of R given by Theorem \6.18[ 



Proof: The algorithm terminates by Theorem 16.181 and returns a differential polynomial P 
in (P[f"''°^, . . . ,p[f"''"^)nQ[u!f"''°', . . . , ui~'"^] with the smallest degree, which is exactly the 
differential resultant. 



{i = 0, ...,n), and I = max^^p/j. So \U\ = Y17=o^iis — Si + 1). In each loop of Step 
3, the complexity of the algorithm is clearly dominated by Step 3.5., where we need to 
solve a system of linear equations P = over Q in cq and Cj^. It is easy to show that 



|co| = ('Ifl^lT') and |c,fc| = (("^+'^"7;;;v;(t-E5''"^'+'^)- Then P = is a linear equation 
system with N = |co| + X^^Lo Sfc=o' l^*'^! variables and M = equations. 



To solve it, we need at most (max{M, N})'^ arithmetic operations over Q, where uj is the 
matrix multiplication exponent and the currently best known to is 2.376. 

Step 3 may loop from d = n + 1 to deg(R) < D, where D is the degree bound of deg(R) 
given by Theorem 16.181 The whole algorithm needs at most 



We will estimate the complexity of the algorithm below. Denote = |uj 




deg(R) 



(max{M,A}) 



2.376 



< 0{{ns + n)2-376[^deg(R)]°('"(^+i))) 



d=n+l 



< 0{{ns + n)^-^'^^[mD] 



0{ln(s+l)) 
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arithmetic operations over Q. In the above inequahties, we assume that (m + l)deg(R) > 
l{ns + n + 1) + n{s + 1), which is generally true. Otherwise, mdeg(R) need to be replaced 
by Ins to give a single exponential complexity bound. Our complexity also assumes an 
0(l)-complexity cost for all field operations over Q. Thus, the complexity follows. □ 

Remark 6.22 One might suggest to use an approach similar to Algorithm DResultant to 

compute the sparse differential resultant and to obtain a complexity bound similar to that 
given in Theorem \6.21[ In this way, a differential polynomial P € sat(R) will be obtained, 
which cannot be proved to be R due to the reason that P might have a higher order and a 
lower degree than that of R. 

With Theorem 16.151 Theorem 11.31 is proved. 

7 Conclusion 

In this paper, we first introduce the concepts of Laurent differential polynomials and Lau- 
rent differentially essential systems, and give a criterion for Laurent differentially essential 
systems in terms of their supports. Then the sparse differential resultant for Laurent differ- 
entially essential system is defined and its basic properties are proved, such as the differential 
homogeneity, necessary and sufficient conditions for the existence of solutions, differential 
toric variety, and the Poisson-type product formulas. Furthermore, order and degree bounds 
for the sparse differential resultant are given. Based on these bounds, an algorithm to com- 
pute the sparse differential resultant is proposed, which is single exponential in terms of the 
order, the number of variables, and the size of the Laurent differentially essential system. 

In the rest of this section, we propose several questions for further study. 

It is useful to represent the sparse differential resultant as the quotient of two determi- 
nants, as done in [111 [T4] in the algebraic case. In the differential case, we do not have such 
formulas, even in the simplest case of the resultant for two generic differential polynomials 
in one variable. The treatment in [6] is not complete. For instance, let f,g be two generic 
differential polynomials in one variable y with order one and degree two. Then, the differ- 
ential resultant for /, g defined in [6] is zero, because all elements in the first column of the 
matrix M{6,n,m) in [6, p. 543] are zero. Although using the idea of Dixon resultants, the 
algorithm in [5l] does not give a matrix representation for the differential resultant. 

From ()28|) , a natural idea to find a matrix representation is trying to define the sparse dif- 
ferential resultant as the algebraic sparse resultant of P = {pf ^ (i = 0, . . . , n, fc = 0, . . . , hi)} 

(i) 

considered as Laurent polynomials in , which will lead to a matrix representation for the 
sparse differential resultant. As far as we know, this is actually very difficult even in the case 
of the resultant for two generic differential polynomials in one variable. 

The degree of the algebraic sparse resultant is equal to the mixed volume of certain 
polytopes generated by the supports of the polynomials [40j or [19, p. 255]. A similar degree 
bound is given in Theorem 11.31 for the differential resultant. We conjecture that the bound 
given in Theorem 11.31 is also a degree bound for the sparse differential resultant. 

There exist very efficient algorithms to compute algebraic sparse resultants [13^ [T^ [T5] , 
which are based on matrix representations for the resultant. How to apply the principles 
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behind these algorithms to compute sparse differential resultants is an important problem. 
A reasonable goal is to find an algorithm whose complexity depends on deg(R), but not on 
the worst case bound of deg(R) which is the case in Algorithm SDResultant. 

In the algebraic case, it is shown that the sparse polynomials Pi (i = 0, . . . , n) can be 
re-parameterized to a new system = 0, . . . , n) with the help of the Newton polygon 
associated with Pj such that vanishing of the sparse resultant gives a sufficient and necessary 
condition for §>i{i = 0, ...,n) to have solutions in C^, where C is the field of complex 
numbers \10\ page 312]. It is interesting to extend this result to the differential case. To do 
that we need a deeper study of differential toric variety introduced in Section 15.31 

As a less important problem, we guess that assuming the first condition in Theorem l5.32| 
the second condition € Span^jajfc — aio : k = 1, ... = 0, ... ,n} is also a necessary 
condition for the system Pj to have a unique solution under the condition of R = in the 
generic case. 
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